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CHAPTER  I 


INTRODUCTION 

Background  for  Investigation 

Methods  of  excavating  soft  and  medium-hard  rocks  have  advanced  in 
the  last  decade  both  in  performance  and  cost  advantages  much  more  than 
methods  of  excavating  hard  rock.  Difficulty  in  penetrating  hard  rocks 
with  mechanical  cutters  has  been  the  primary  factor  delaying  the  advance¬ 
ment  of  more  nearly  continuous  methods  of  hard  rock  excavation.  The 
problems  of  excavating  hard  rock  by  boring  machines  (bearing  and  cutter 
wear  due  to  high  thrust)  for  all  types  of  openings  are  well-known,  and 
currently  drill  and  blast  methods  are  usually  required.  New  methods  of 
disintegration  are  desired  which  may  approach  being  continuous,  and  pos¬ 
sess  other  favorable  characteristics  which  will  make  them  technically 
feasible  and  economically  applicable. 

Thus,  rapid  excavation  of  hard  rock  by  mechanical  cutters  alone  is 
not  presently  feasible,  and  currently  used  drill  and  blast  methods  of 
excavation  are  cyclic  in  nature  and  lack  the  advantages  of  continuous 
operation  and  smooth  boring  ability.  When  heat  sources  are  used  to  ex¬ 
cavate  hard  rock  by  melting,  very  large  quantities  of  energy  are  required, 
and  for  removal  of  rock  by  spallation  the  process  is  usually  limited  to 
types  of  spallable  rocks  composed  of  a  minimum  of  ten  percent  quartz. 

The  combination  of  thermal  weakening/spalling  and  mechanical  disintegra¬ 
tion  (thermal -mechanical  fragmentation)  appears  to  offer  a  promising  poten¬ 
tial  for  hard  rock  excavation  because  of  the  non-cyclic  nature,  methods 
of  application,  potential  economic  feasibility,  and  possible  applicability 
to  hard  rocks  not  containing  quartz. 
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This  research  project  to  study  the  excavation  of  hard  rock  using 
combined  mechanical  and  thermal  methods  of  rock  removal  is  concerned 
with  theoretical  studies  of  the  process  of  heat  weakening  and  spalla¬ 
tion,  including  pertinent  aspects  of  thermoelastic  stresses,  heat  energy 
coupling  efficiency,  rock  surface  temperature,  and  temperature  distribu¬ 
tion  with  depth.  Theoretical  studies  of  disintegration  include  surface 
chipping  for  particular  indentors,  kerf  ridge  removal,  energy  require¬ 
ments,  and  estimation  of  efficiency  for  heat  assisted  mechanical  chipping. 
Types  of  rocks  with  dissimilar  composition  and  spalling  characteristics 
were  selected  for  experimental  examination  to  determine  the  applicability 
of  the  theoretical  studies.  Several  important  technical  parameters  of 
a  technique  utilizing  pneumat’c  chippers  and  rock  weakening  by  heat  were 
investigated  experimentally,  and  pertinent  theories  of  thermoelasticity, 
heat  transfer,  and  impact  by  two  types  of  bits  were  studied  in  relation 
to  fragmentation. 

The  experimental  difficulties  anticipated  in  the  research  were  test 
equipment  design,  and  devising  techniques  for  collection  of  large  fragments 
as  well  as  small  particles  broken  from  the  large  blocks  of  rock.  Evalua¬ 
tion  of  specific  energy,  surface  energy,  and  efficiency  of  energy  trans¬ 
fer  were  among  the  major  technical  problems  to  be  solved. 

The  data  and  theory  developed  has  added  vital  knowledge  of  the  response 
of  rock  to  combined  high  temperature  and  impact.  The  experiments  were 
designed  to  lead  to  the  development  of  the  prototype  of  an  excavation 
machine. 

The  disintegration  of  soft  rock  by  tunnel  bor'ing  machines  exemplified 
many  of  the  characteristics  of  excavation  methods  which  are  desired  for 
shafts  and  tunnels  in  more  effective  excavation  of  hard  abrasive  rock. 


In  general,  some  of  the  more  novel  types  of  rccK  breakage  that  are  being 
currently  studied  are  aimed  at  the  development  of  rapid  methods  of  tunnel 
boring  several  years  in  the  future. 

Rock  Properties 

Considerable  laboratory  research  has  been  performed  on  the  properties 
of  rock  which  affect  their  resistance  to  disintegration.  Thermal  and 
chemical  means  of  weakening  rock  have  been  devised  on  small  laboratory 
specimens,  but  none  of  these  appear  to  have  been  applied  to  large  blocks 
or  to  rock  in  situ.  The  science  and  technology  of  mechanical  and  thermal 
fragmentation,  however,  has  advanced  far  enough  that  it  has  proven  profit¬ 
able  to  apply  recent  findings  to  disintegration  of  rock  on  a  larger  scale. 
Review  and  analysis  of  existing  theory  and  data  (Refs.  1-3)  indicate  that 
mechanical  properties  of  rock  such  as  strength,  moduli,  hardness,  etc. 
measured  by  conventional  methods,  are  often  only  indirectly  related  to 
rock  resistance  to  disintegration  processes.  Thermal  properties,  however, 
can  be  utilized  in  extending  the  applications  of  the  theory  of  thermal 
stresses,  heat  transfer,  and  related  phenomena. 

While  extensive  studies  have  been  made  of  the  fracture  of  metals  and 
rock  under  static  loads,  only  a  relatively  small  amount  of  research  has 
dealt  with  rock  fracture  under  dynamic  loads.  Effective  dynamic  tensile 
strengths  have  been  found  to  be  5  to  13  times  as  great  as  static  (Ref.  4). 
However,  most  fracture  studies  were  performed  by  means  of  the  Hopkinson 
bar  which  breaks  the  rock  by  means  of  a  pulse  reflected  from  a  free  surface. 

The  ambient  temperature  static  data  available  for  the  rocks  tested 
in  this  project  are  listed  in  Table  (A-l)  of  Appendix  A. 


Excavation  System  Characteristics 


The  critical  characteristics  of  an  effective  method  of  rock  disintegration 
excavation  system  are: 

1.  Capable  of  operation  in  all  types  of  rock  and  rock  conditions 

2.  Minimum  shutdown  time 

3.  Continuous  operation  to  approach  simultaneous  breaking,  loading, 
and  haulage 

4.  Economical  operation  costs 

5.  Favorable  or  tolerable  environmental  factors 

6.  Reasonable  capital  costs  and  depreciation 

7.  Reasonable  overhead  costs 

The  first  step  in  the  evaluation  of  a  prospective  method  of  excava¬ 
tion  is  the  verification  of  its  technical  feasibility.  If  the  investiga¬ 
tion  of  the  fragmentation  basic  process  or  processes  establishes  that 
critical  parameters  are  favorable,  then  the  feasibility  of  incorporating 
them  into  an  operating  system  and  the  economics  of  the  system  must  be 
established.  Important  specific  factors  which  influence  the  overall 
evaluation  are  as  follows: 

Rock  Disintegration  and  Excavation  Parameters 
A.  Characteristics  of  technique 

1.  State  of  research  and  development 

a.  Theory 

b.  Research  methods 

c.  Research  data 

d.  Development 


2.  Energy  requirements 

a.  Type,  electrical,  mechanical,  chemical,  etc. 

b.  Source,  fuel,  explosives,  etc. 

c.  Degree  of  fragmentation 

d.  Efficiency  of  energy  transfer 

e.  Efficiency  of  energy  usage  -  specific  energy  of  disintegration 

f.  Overall  efficiency  of  process 

g.  Possible  use  of  waste  energy 

3.  Fracturing  process 

a.  Geometry  of  <_.;ergy  application,  external,  internal,  etc. 

b.  Temperature  levels  and  effect 

c.  Character  of  energy  application  -  impact,  thermal  expansion,  etc. 

d.  Relationships  to  other  techniques  and  possible  complementary  uses 

e.  Degree  of  disintegration  (2c  above) 

Application  of  Technique 

1.  Ability  to  control  separation  of  rock  masses 

2.  Ability  to  preserve  strength  of  peripheral  rock  in  opening 

3.  Ability  to  make  smooth  opening 

4.  Power  requirements  and  efficiency 

a.  At  the  rock  interface 

b.  Overall 

5.  Environmental  factors,  noise,  fumes,  etc. 

6.  Adaptability  to  changes  in  properties  and  structure  of  rock 

7.  Equipment  requirements  for  application 

a.  At  face  -  state  of  development 

b.  Back  up  -  state  of  development 

c.  Limiting  requirements 
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Applicability 

a.  To  horizontal  openings 

b.  To  vertical  or  inclined  openings 

9.  Character  of  method  using  technique 

a.  Cyclic 

b.  Semi -cyclic 
i_ .  Continuous 

d.  Down  time  of  continuous  method 

10.  Thrust  requirements 

11.  Metallurgical  requirements  of  machines 

12.  Research  and  development  requirements 

a.  Type  of  additional  research  required 

b.  Estimate  of  time  for  research  and  development 

c.  Estimate  of  relative  costs 

13.  Possible  use  of  complementary  methods 

14.  Relationships  of  disintegration  to  mucking,  haulage  and  hoisting 

15.  Estimated  costs  -  unit  and  overall 

a.  Research  and  development 

b.  Capital  investment  in  operating  machine 

c.  Operating  costs,  energy,  labor,  etc. 

d.  Support  costs 

e.  Overhead 

f.  Other 

16.  Constraints 
Excavation  Techniques 

The  most  critical  factor  in  the  functioning  of  a  disintegration 
technique  is  almost  always  the  effective  transfer  of  breakage  energy 
into  the  rock.  Of  almost  equal  importance  is  the  efficiency  and  cost 
of  producing  the  energy. 
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Processes  of  rock  disintegration  by  drilling  are  similar  to  those 
which  are  used  or  proposed  for  use  in  full  face  disintegration  of  tunnels, 
shafts,  etc.  Characteristics  of  novel  methods  have  been  summarized  (Ref.  4} 
according  to  status,  rock  removal  mechanism,  specific  energy,  maximum  power 
to  rock  and  maximum  potential  drilling  rate,  as  well  as  those  for  convention¬ 
al  drills. 

For  full  face  disintegration  versus  drill  hole  cutting  there  are  fewer 
restrictions  on  the  power  which  can  be  delivered  at  the  rock.  Field  tests 
have  been  made  of  several  novel  drills  but  the  only  ones  in  operation  are 
conventional  rotary  and  jet  piercing,  although  the  Russians  claim  to  have 
used  explosive  drilling  successfully  to  limited  depths  in  some  types  of 
rock.  In  addition,  shaped  explosive  charges  may  be  employed  to  produce 
high  velocity  metallic  jets  which  are  effective  in  penetration.  Some 
specific  energy  values  for  tunnel  boring  machines  are  available  and  much 
of  the  fragmentation  energy  is  utilized  in  producing  very  fine  particles. 

The  values  of  specific  energy  for  drilling  by  Maurer  (Ref.  4)  give 
only  the  energy  to  fracture  the  rock,  and  not  the  total  energy  require¬ 
ments  or  the  efficiency  of  the  system.  Approximations  for  total  energy 
requirements  are  available  for  some  of  the  preferred  methods  (Table  1). 

Table  1 

Estimates  of  Total  Energy  Consumption* 
joules/cc 

Chemical  explosive  10-30  total  available  energy 

Pneumatic  moil  30-110  total  required  energy 

Flame  jet  15,000+  spalling 

Electron  beam  12,000+  melting 

Laser  12,000+  heat  weakening  and  melting 

Radiant  heat  6,000  heat  weakening  and  spallation 

*  From  estimates  given  at  Rapid  Excavation  Conference,  Deerfield,  Mass., 

July  20-25,  1969. 
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A  study  of  the  feasibility  of  flame  jet  tunneling  by  the  United  Air¬ 
craft  Research  Laboratories  (UARL)  (Ref.  5)  concluded  that  it  may  be  feas¬ 
ible  in  rocks  with  good  spal labi 1 ity.  One  geometry  favorable  for  suggested 
breakage  was  by  means  of  annular  kerfing  and  fracture  of  large  blocks  of 
the  kerf  with  a  breaker  arm.  Thus,  the  overall  specific  energy  for  com¬ 
bined  thermal  spall inq  and  breaking  was  believed  to  be  economical  even 
though  the  specific  energy  of  flame  jet  cutting  may  vary  from  1500  to 
6000  joules/cc.  The  predicted  rate  of  advance  for  a  30  foot  tunnel,  however, 
was  comparable  only  to  rates  achievable  by  conventional  drilling  and  blast¬ 
ing. 

For  complete  small  particle  fragmentation  by  flame  jets  major  problems 
exist  in  environmental  control  in  that  noxious  fumes  and  high  temperatures 
are  continuously  created  by  the  flames  which  cause  the  spallation.  The 
attendant  engineering  problems  and  cost  of  creating  acceptable  environ¬ 
mental  conditions  increase  rapidly  with  the  increase  in  size  of  tunnel  so 
that  the  rate  of  advance  makes  tunnels  larger  than  30  foot  diameter  im¬ 
practical  (Fig.  7,  p.  13,  Ref.  5).  Mechanical  breakage  of  rock  between 
kerfs  was  not  experimentally  evaluated. 

The  investigators  at  UARL  also  evaluated  the  spallability  of  some  rock 
types  in  the  Northeast  Corridor,  i.e.,  between  Washington,  New  York,  and 
Boston,  About  one-third  of  them  had  high  spallability  indices.  Thus, 
th’s  method  apparently  has  somewhat  limited  applicability  for  tunneling 
in  this  geographical  area.  The  presence  of  mud  seams  or  water  in  the  face 
of  a  tunnel  could  stop  spallation  altogether. 

Infrared  heat  units  have  been  recently  adapted  for  use  at  a  copper 
mine  for  heat  weakening  and  spallation  of  a  shale  ore  bed.  Their  effici¬ 
ency  of  heat  generation  and  transfer  into  rock  are  both  reported  to  be 
good.  Temperatures  are  high  and  the  energy  has  excellent  penetrability. 
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i.e.,  usual  heat  diffusion  processes  are  not  involved  to  the  same  degree 
at  the  rock-air  interface.  The  high  intensity  radiant  heat  units  employed 
for  this  type  of  heating  are  stock  items,  but  depth  of  penetration  of  infra¬ 
red  heat  into  rock  has  not  been  evaluated. 

To  the  present,  laser  beams  have  been  employed  only  in  the  heat  weaken¬ 
ing  of  small  beams  and  cutting  of  slabs  of  rock  in  the  laboratory.  At 
their  current  energy  levels  their  overall  efficiency  of  transfer  of  energy 
to  rock  ranges  from  3  to  6  percent  for  heating,  and  as  high  as  30  percent 
for  melting.  Their  dispersed  thermal  effects  are  much  the  same  as  radiant 
heating,  although  more  concentrated. 

Thermal -mechanical  Fragmentation 

Thermal  spalling  may  be  caused  either  by  (1)  differential  thermal  ex¬ 
pansion  of  mineral  constituents,  (2)  thermal  decrepitation,  (3)  expansion 
of  entrapped  water,  or  (4)  thermal  destruction  of  the  binding  material  be¬ 
tween  grains. 

The  properties  of  rock  significant  in  thermal  weakening  or  disintegra¬ 
tion  (Ref.  6)  are  thermal  expansion,  temperature  phase  changes,  diffusivity 
or  conductivity,  heat  capacity,  heat  of  fusion,  heat  of  vaporization,  de¬ 
crepitation,  and  change  in  chemical  form  (Ref.  2).  Some  rocks  with  above 
ten  percent  quartz  spall  because  of  thermal  expansion  and  phase  changes 
(Ref.  6).  Dolomites  may  spall  because  of  a  two-stage  decomposition. 

Shales  containing  cognate  water  have  been  found  to  spall  with  considerable 
energy  (Ref.  7),  and  some  sandstones  experience  grain-by-grain  disintegra¬ 
tion  because  of  weakening  of  the  binding  material  by  heat  (Ref.  2). 

The  response  of  rocks  to  heat  varies  with  the  intensity  of  heat  (tem¬ 
perature),  the  length  of  time  of  application,  the  nature  of  the  energy 
which  heats  the  rock,  and  the  composition  and  structure  of  the  rock. 

The  same  factors  which  cause  spallation,  which  may  be  due  to  thermal  ex¬ 
pansion,  decrepitation,  expansion  of  entrapped  water  or  destruction  of 
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mineral  grain  binders,  may  also  cause  weakening  of  the  rock  structure. 

Such  weakening  changes  the  resistance  of  the  rock  to  disintegration  by 
various  mechanical  processes  such  as  chipping,  impact  cutting  or  abrasion. 

Fragmentation  of  solid  rock  caused  by  thermal  expansion  is  produced 
by  a  set  of  stresses  imposed  on  the  rock  material  in  a  given  set  of  con¬ 
ditions.  Modification  of  the  response  of  the  rock  may  result  if  the  im¬ 
posed  stress  state  and/or  the  properties  of  the  rock  material  are  altered. 
Thermal -mechanical  fragmentation  alters  the  rock  material  as  well  as  in¬ 
creasing  the  local  stress  state.  The  material  properties  are  changed  dur¬ 
ing  heating  because  of  microcracks  which  form,  and  because  the  basic  me¬ 
chanical  properties  are  different  at  elevated  temperatures.  The  stress 
state  is  changed  by  the  superposition  of  thermal  and  mechanical  stresses 
and  by  the  fact  that  the  original  continuous  medium  becomes  discontinuous 
after  thermal  treatment. 

Mechanical  stresses  are  produced  when  indentors  are  forced  or  impacted 
against  a  rock  surface.  Wedge-  and  spherical -shaped  tungsten  carbide  in¬ 
dentors  are  commonly  used  for  hard  rock  and  the  force  is  usually  supplied 
by  pneumatic  hammers  or  hydraulic  loading  of  rolling  cutters.  A  pneumatic 
hammer  was  used  in  this  research.  Fracturing  occurs  when  induced  stresses 
exceed  the  tensile  or  shear  strength  of  the  rock.  Beneath  the  area  of 
indentor  contact  a  zone  of  finely  crushed  rock  is  created  from  which  sig¬ 
nificant  fractures  are  initiated.  These  fractures  propagate  along  curved 
paths  to  the  surface,  breaking  off  chips  of  rock.  The  fine  crushing  re¬ 
sults  from  intense  stresses  created  by  the  initial  approximate  line  or 
point  contact  of  the  indentor.  It  should  be  noted  that  the  indentor  im¬ 
pacts  on  softened  or  fractured  material  if  thermal  treatment  precedes  the 
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mechanical  device.  However,  fruitful  theoretical  studies  of  randomly  dis¬ 
continuous  material  are  not  possible  with  the  present  state-of-the-art. 


The  assumption  that  the  material  is  continuous  permits  a  reasonable  first 
approximation  that  allows  important  process  parameters  to  be  investigated. 

Ncnuniform  heating  of  rock  produces  thermal  stresses  by  the  interaction 
of  forces  from  different  local  sections  of  the  rock  that  expand  different 
amounts,  either  due  to  rock  geometry  or  the  following  factors.  Because  of 
the  low  thermal  conductivity  of  rock  (Table  I,  Appendix  A)  the  surface  may 
be  at  the  spalling  temperature  while  only  1/8  inch  or  less  below  the  rock 
may  be  at  its  ambient  temperature.  This  sharp  thermal  gradient  is  a  primary 
factor  in  creating  the  thermal  stresses  that  cause  fracture.  The  nonhomo- 
geneous  or  orthotropic  composition  of  rock  assists  the  formation  of  local¬ 
ized  thermal  stresses  because  individual  grains  and  crystals  of  the  rock 
have  different  coefficients  of  thermal  expansion,  and  orthotropic  rocks  ex¬ 
pand  different  amounts  along  twj  perpendicular  axes.  Other  important  pro¬ 
cesses  in  creating  thermal  stresses  are  the  phase  changes  that  occur  in 
some  minerals  with  increase  in  temperatures.  Theii  importance  is  charac¬ 
terized  by  quartz,  which  undergoes  0.82  percent  volumetric  expansion  during 
the  alpha -to-beta  phase  transition  at  573  degrees  C.  Minerals  surrounding 
the  quartz  crystals  constrain  the  expansion  and  thus  produce  high  thermal 
stresses.  High  thermal  stresses  are  also  created  in  relatively  uniform 
rock  by  the  geometry  of  the  heat  source  and  its  application  on  flat  surfaces, 
in  corners  or  holes  in  the  rock. 

Rocks  of  many  types  permanently  lose  strength  after  being  heated,  es¬ 
pecially  if  the  rock  is  heated  to  temperatures  above  600  degrees  C.  Reasons 
for  this  permanent  reduction  in  strength  are  chemical  changes,  spalling 
and/or  anisotropic  thermal  expansion,  and  the  formation  of  microcracks 
along  grain  boundaries.  Microcracks  form  as  the  rock  crystals  and  grains 
expand  and  contract  due  to  heating  and  cooling.  In  this  weakened  condition 
a  rock  may  be  more  susceptible  to  removal  by  some  mechanical  processes. 
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If  mechanical  chipping  is  applied  to  the  rock  while  thermal  stresses 
are  still  present,  crack  initiation  and  propagation  may  be  enhanced  due 
to  the  superposition  of  thermal  and  mechanical  stresses.  The  direct  com¬ 
bination  of  these  two  processes  has  been  referred  to  as  thermal -mechanical 
fragmentation. 

Plan  and  Objectives  of  Research 

The  research  program  was  divided  into  three  general  areas:  (1)  theo¬ 
retical  studies  of  heat  energy  coupling  and  the  mechanics  of  heat  weaken¬ 
ing  of  rock  by  surface  heating;  (2)  theoretical  studies  of  mechanical  dis¬ 
integration  (chipping)  and  the  effects  of  heat  thereon;  and  (3)  laboratory 
investigations  of  the  important  parameters  of  heat  assisted  near  surface 
chipping.  Significant  deep  kerfing  tests  were  performed  for  three  rocks. 

The  heat  energy  coupling  studies  were  concerned  primarily  with  the 
efficiency  of  heat  transfer  from  high  temperature  moving  flame  sources 
(such  as  flame  jets  arid  acetylene  torches)  into  the  surface  of  the  rock. 

The  classical  equations  for  static,  quasi-static  and  moving  sources  have 
been  employed  to  describe  the  transfer  phenomena  into  homogeneous  solids 
and  the  equation  for  the  quasi -static  case  was  analyzed  employing  the  para¬ 
meters  of  the  flame  jet-rock  system.  For  simple  conditions,  if  the  tempera¬ 
ture  distribution  is  known  the  thermoelastic  stress  field  in  the  rock  may 
then  be  calculated  and  possible  modes  of  failure  predicted,  high  stress 
concentrations  being  created  by  large  thermal  gradients. 

The  existing  theories  of  penetration  of  indentors  into  rock  required 
extension  to  apply  theories  of  plasticity  to  the  reaction  rf  rock  to  pres¬ 
sure  from  wedge  indentors  of  various  angles  and  friction  at  the  rock  inter¬ 
face,  for  flat-edged  indentors  and  for  hemispherical  indentors.  Preheat¬ 
ing  of  the  rock  and  microfracturing  decreased  the  effectiveness  of  the  in¬ 
dentors  in  removing  chips  from  some  types  of  rock  but  also  decreased  the 
friction  and  wear  in  the  indentors  in  some  cases. 
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Experimental  work  was  conducted  to  evaluate  the  basic  parameters  in¬ 
volved  in  the  above  theories  and  the  extent  to  which  the  theories  can  be 
applied  to  six  types  of  rock  (Appendix  A).  This  included  (1)  the  evalua¬ 
tion  of  energy  requirements  for  fragmentation  with  and  without  thermal 
weakening;  (2)  investigation  of  important  parameters  of  fragmentation  in¬ 
cluding  width  of  kerf,  optimum  indexing,  and  size  and  shape  of  chippers; 

(3)  temperature  measurements  of  rock  during  thermal  weakening;  and  (4) 
the  development  of  tests  to  measure  resistance  to  disintegration,  and  a 
comparison  of  results  between  mechanical  and  thermal -mechanical  methods  of 
disintegration. 

Tests  were  conducted  with  a  commercially  available  pneumatic  chipper 
(jackhammer)  equipped  with  cross-wedge  and  button  bits.  Heat  tests  were 
conducted  largely  with  a  standard  flame  jet,  and  a  few  tests  were  run  with 
acetylene  welding  torches.  A  drop  crushing  test  similar  to  that  employed 
by  Gross  and  Zimmerly  (Ref.  8)  offers  promise  as  a  means  of  evaluating  rock 
resistance  to  impact  and  thermal  breakage  as  well  as  for  calculating  values 
of  breakage  (surface  energy)  by  means  of  Rittinger's  law  of  crushing. 

Variation  in  some  parameters  was  somewhat  restricted  by  the  character¬ 
istics  of  the  heat  source  (flame  jet)  and  pneumatic  chipper  (jackhammer) 
chosen.  That  is,  the  flame  jet  operates  essentially  at  one  level  of  heat 
output  which  is  required  to  maintain  flame  stability,  and  heat  flux  may 
be  varied  by  changing  the  standoff  distance.  Likewise,  the  energy  output 
of  the  jackhammer  can  be  varied  only  over  a  small  rang:,  because  of  its 
optimum  operating  characteristics,  and  energy  transmitted  to  the  rock  is 
difficult  to  measure  accurately. 
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CHAPTER  II 


HEAT  ENERGY  COUPLING 

Heat  Coupling  -  Moving  Flame  Source 

For  one  to  evaluate  the  process  of  thermal  fragmentation  of  rock 
materials  resulting  from  the  application  of  a  heat  source  the  temperature 
distribution  inside  the  rock  material  must  be  determined.  The  heat  energy 
coupled  into  a  rock  mass  from  an  energy  source  must  be  known  to  determine 
whether  the  source  can  be  economically  used  for  thermal  weakening,  spall¬ 
ing,  or  splitting  of  rock  materials.  A  combined  experimental  and  numeri¬ 
cal  method  was  used  to  estimate  the  thermal  energy  crossing  a  rock  surface 
from  a  flame. 

Measurements  of  the  surface  and  subsurface  temperatures  were  made  to 
define  the  heat  penetration  from  a  heat  source  into  a  rock  material.  From 
the  temperature  measurements  the  heat  energy  coupling  efficiency  was  esti¬ 
mated. 

General  Theory 

Heat  Penetration.  The  temperature  distribution  within  the  rock  material 
can  be  determined  by  either  solution  r  at  transfer  equations  directly 

or  experimental  measurement  of  the  tt  s.  Rocks  are  complex  materials 

with  inherent  inhomogeneity  and  anisoti  .  well  as  temperature-dependent 
properties,  and  because  of  their  character  the  heat  conduction  equation 
becomes  nonlinear.  The  heat  transfer  equation  for  a  moving  heat  source 
over  a  rock  surface  is  not  only  nonlinear  but  transient  and  three-dimensional . 
A  solution  for  a  problem  of  this  type  is  not  available.  An  experimental 
approach  was  taken  to  measure  the  surface  temperature  and  selected  sub¬ 
surface  temperatures  of  the  rock.  The  experimental  temperature  measure¬ 
ment  then  defines  the  deplh  of  penetration  of  the  heat. 
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V 

Several  rock  types  were  tested  to  determine  if  measurable  differences 
in  thermal  response  can  be  observed  (as  a  function  of  mineral  content). 

The  surface  temperatures  were  correlated  with  measurements  made  at  various 
depths  in  order  to  better  define  the  heat  penetration  characteristics. 

Similar  tests  were  run  by  Browning,  et  al.,  (Ref.  9)  using  embedded  thermo¬ 
couples.  Results  of  measurements  made  for  a  flame  jet  heat  source  will 
later  be  compared  with  similar  measurements  for  other  heat  sources  during 
the  second  phase  of  this  research  project. 

Heat  Energy  Coupling 

The  calculation  of  the  heat  energy  coupling  efficiency  requires  the 
calculation  of  the  heat  energy  entering  the  rock  surface.  The  heat  trans¬ 
ferred  to  a  plane  surface  normal  to  a  gas  jet  was  studied  by  J.E.  Anderson 
and  E.F.  Stresino  (Ref.  10)  for  various  torches*.  Their  experimental  work 
was  performed  using  a  copper  plate  which  was  cooled  by  water  flowing  past 
one  side  of  the  copper  plate.  The  heat  transferred  to  the  plate  was  cal¬ 
culated  from  the  rise  in  the  temperature  of  the  water.  Copper  was  used 
because  it  is  easily  machined  and  has  a  high  conductivity.  The  same  pro¬ 
cedure  would  not  be  possible  with  rock  because  of  its  low  conductivity 
and  brittleness.  Other  investigators  have  studied  the  heat  transferred 
between  a  gas  jet  and  a  plane  surface  and  attempted  to  obtain  a  correla¬ 
tion  of  gas  properties  with  heat  transferred.  In  all  of  these  cases  the 
plane  surface  was  metallic,  however,  and  the  material  properties  remained 
relatively  constant  over  the  range  of  temperatures  studied. 

The  methods  described  above  are  not  applicable  for  rock  materials  be¬ 
cause  of  their  low  conductivity  and  the  difficulty  in  machining.  It  was 
therefore  necessary  to  develop  a  new  method  of  calculating  the  heat  transfer 

*Note:  A  recent  technical  note  titled  "Impingement  Heat-Transfer  Rates 

frum  Torch  Flames"  by  E.H.  Sriulte,  Jl.  of  Heat  Transfer,  pp.  231-233, 
May  1972  also  investigates  the  heat  transfer  from  a  flame  torch  to 
a  plane  surface. 
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from  a  flame.  Since  a  new  method  was  required  the  temperature  distribu-^ 
tion  obtained  from  the  heat  penetration  study  was  used  to  estimate  the 
heat  entering  the  rock  material.  The  problem  of  calculating  the  heat 
entering  a  body  from  the  temperature  history  inside  the  body  is  known  as 
the  inverse  problem  of  heat  conduction. 

Stolz  (Ref.  11)  and  others  (Refs.  12,  13,  14,  15)  have  solved  one¬ 
dimensional  inverse  problems.  Beck  (Ref.  16)  used  a  method  of  nonlinear 
estimation  to  solve  a  one-dimensional  nonlinear  inverse  heat  condition 
problem.  The  nonlinearity  was  due  to  temperature-dependent  material  pro¬ 
perties.  As  in  the  direct  method  the  complex  material  properties  of  rock 
and  introducing  more  than  one  dimension  makes  the  solution  of  the  inverse 
problem  very  difficult.  The  inverse  method  is  mathematically  ill-posed 
(Refs.  15,  17)  and  as  such  can  develop  oscillations  in  the  solutions. 

The  temperature  histories  obtained  experimentally  are  three-dimensional 
in  the  space  coordinates,  thus  giving  rise  to  a  three-dimensional  inverse 
problem.  Solving  the  three-dimensional  nonlinear  inverse  problem  is  more 
difficult  than  solving  the  direct  problem  of  heat  conduction. 

During  the  second  phase  of  this  investigation  further  consideration 
will  be  given  to  developing  a  method  for  converting  the  three-dimensional 
temperatures  into  surface  data.  A  possible  technique  for  eliminating  the 
difficulties  of  inverting  a  three-dimensional  nonlinear  problem  would  be 
to  devise  a  one-dimensional  experimental  test. 

For  calculation  of  specific  energy  of  breakage  for  the  fragmentation 
process  an  estimate  of  the  heat  energy  coupling  was  required.  An  approxi¬ 
mation  of  the  energy  entering  che  body  was  made  from  the  specific  heat 
and  the  measured  temperatures.  Calculation  of  the  heat  output  of  the  torch 
was  made  based  or.  the  heat  of  combustion  and  the  flow  rate  of  the  fuel. 


-* 

l 

t 

The  heat  transferred  from  the  flame  jet  into  the  rock  passes  through 
three  regions: 

1.  Region  I.  From  the  torch  to  the  gas  film  -  Forced  Convection, 

2.  Region  II,  Through  the  gas  film  to  the  surface  -  Conduction. 

3.  Region  III.  Through  the  rock  material  -  Conduction.  Figure  (1) 
illustrates  the  three  regions.  Anderson  and  Stresino  (Ref.  10)  state 
that  radiation  from  their  torches  to  a  metallic  surface  was  small  when 
compared  to  the  heat  transferred  by  convection.  The  same  assumption  was 
made  in  this  report  for  the  flame  jet  impinging  on  a  rock  surface.  Three 
factors  (Ref.  10)  which  complicate  the  analysis  of  a  flame  jet  impinging 
on  a  surface  are: 

"Reliable  values  of  the  effective  thermal  conductivity  in  the 
presence  of  dissociation-recombination  reactions  in  a  chemi¬ 
cally  complicated  gaseous  mixture  are  usually  not  available; 
the  heat-transfer  intensity  distribution  obtained  when  a  gas 
jet  impinges  on  a  flat  plate  is  imperfectly  understood;  and 
in  practical  cases  such  as  those  involving  flames  steep  tem¬ 
perature  gradients  often  exist  in  the  gas  stream." 

Metzger  (Ref.  18)  studied  the  heating  of  a  surface  by  a  high  velocity 
air  jet.  For  a  first  approximation  the  results  for  an  air  jet  can  be  used 
to  analyze  a  flame  jet. 

One  of  the  primary  factors  describing  the  efficiency  of  heat  transfer 
from  a  hot  gas  to  a  solid  is  the  value  of  the  heat  transfer  coefficient. 

An  equation  by  Metzger  (Ref.  18)  described  by  Carstens  (Ref.  19)  utilizes 
the  relationship  between  Stanton,  Reynolds,  and  Prandtl  numbers  and  the 
model  geometry  for  a  circular  nozzle: 

(Nst)|NRe,*-,'”(Np,)"-e3  -  o^(r)  '  0) 
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The  symbols  employed  are  defined  below.  The  above  equation  may  be  re¬ 
written  in  terms  of  basic  quantities  involved  in  the  dimensionless  fac¬ 


tors  as: 


-0.434 


This  may  be  rearranged  and  with  R  and  V  changed  to  units  of  in.  and 
ft/sec,  respectively,  and  the  transfer  coefficient  expressed  as: 


h  =  0.06225  Cng 


(vo)"'  (v) 


Carstens,  et  al . ,  (Ref.  19)  utilized  a  plane  front  equation  to  define 
the  temperature  distribution  within  the  rock,  and  the  quantity  of  heat 
flowing  into  the  rock.  The  latter  was  integrated  over  an  area,  approxi¬ 
mated  for  short  periods  of  time  and  summed  for  the  total  time  to  approxi¬ 
mate  a  transient  condition. 


Temperature  equation: 


.y  +  h/ict 
2vk£  kr  . 


The  rate  of  heat  transfer  is  given  by: 


q  =  2nR*h  (Tf  -  Ts) 
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The  total  heat  flow  is  calculated  from: 

t 


Nomenclature  for  the  above  equations  is  as  follows: 

Cp  Specific  heat  at  constant  pressure  of  exhaust  gas  -  Btu/lb-F 

D  Exit  diameter  of  je*  nozzle  -  ft 

g  Mass  velocity  of  exhaust  at  nozzle  exit  -  lb/hr-ft2 

h  Heat  transfer  coeff icier*  in  wall  film  -  Btu/hr-ft2-F 

k  Thermal  conductivity  of  exhaust  gas  -  Btu/hr-ft-F 

kR  Thermal  conductivity  of  rock  -  Btu/hr-ft-F 

Npr  Prandtl  number  -  dimensionless 

NRe  Reynolds  number  -  dimensionless 

Nst  Stanton  number  -  dimensionless 

q  Rate  of  heat  transfer  into  rock  surface  -  Btu/hr 
Q  Total  heat  transferred  into  rock  in  time,  t  -  Btu 
Rq  Radius  of  rock  wall  target  -  in. 
t  Time  -  sec 

Tp  Exhaust  gas  temperature  at  wall  -  °F 

Tr  Equilibrium  rock  temperature  prior  to  heat  addition  -  °F 

T  Temperature  in  rock  at  distance  y  from  surface  -  °F 

y  j 

Ts  Temperature  at  rock  surface  -  °F  / 

V  Exhaust  gas  velocity  at  nozzle  exit  -  ft/sec 

y  Linear  distance  into  rock  from  surface  -  ft  ' 

k  Thermal  diffusivity  of  rock  -  ft2/hr 

Pg  Exhaust  gas  density  -  lb/ft3 

u  Exhaust  gas  viscosity  -  lb/hr-ft 
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From  the  above  equations  it  can  be  stated  that  the  instantaneous 
rate  of  heat  transfer  into  the  rock  depends  upon  the  heat  transfer  co¬ 
efficient,  and  the  temperature  difference  between  the  rock  surface  and 
the  gas  temperature  at  the  wall.  When  a  flame  jet  is  initially  applied 
to  the  surface  of  a  rock  at  ambient  temperature,  the  temperature  of  the 
gas  film  at  the  surface  is  much  lower  than  the  exhaust  gas  temperature. 
Also,  as  the  rock  surface  is  heated,  the  temperature  difference  between 
the  rock  and  the  gas  film  decreases  with  time,  even  though  the  gas  film 
temperature  may  also  rise. 

The  rise  in  temperature  of  the  surface  of  the  rock  depends  in  turn 

upon  the  rate  at  which  heat  is  conducted  into  the  rock  from  the  surface, 

i.e.,  upon  kR  and  k,  both  of  which  may  also  vary  with  temperature.  The 

kR 

relationship  between  these  thermal  factors  is  k  =  -5—  where  CD  is 

p  k  K 

the  heat  capacity  of  the  rock. 

Inasmuch  as  the  spall  ability  and  tendency  to  generate  microfractures 
is  largely  dependent  upon  the  temperature  gradient,  a  low  value  of  conduc¬ 
tivity  would  appear  to  be  one  of  the  major  contributing  factors  to  rock 
spallability  and  microfracture. 

Rosenthal  (Ref.  20)  was  able  to  obtain  a  good  correlation  between  a 
moving  point  source  and  experimental  data  for  mild  steel  and  an  arc  welder. 
Surface  temperature  curves  with  high  gradients  and  a  high  peak  temperature 
were  obtained  (Fig.  2).  An  arc  welder  will  produce  a  heat  source  not  un¬ 
like  a  point  source  while  a  flame  jet  is  a  source  distributed  over  a  finite 
area.  The  size  of  the  area  is  dependent  on  the  standoff  distance  of  the 
torch  from  the  surface.  Rosenthal  made  a  quasi -stationary  assumption  and 
derived  an  equation  for  the  temperature  distribution  from  a  mv/ing  point 
source  into  a  half  space  as 

T  -  ’o  -A  ***  b  e-XVr‘  17) 
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where 

c  =  specific  heat 
k  =  thermal  conductivity 
q  =  rate  of  heat  transfer  into  rock  surface 
r*  =  ^i+ya+z2 

Tq  =  initial  or  reference  temperature  of  surface 
v  =  velocity  of  the  heat  source 
x  =  di.^tion  of  movement  of  the  heat  source 
e  =  x-vt 
\  =  cp/2k 
p  =  density 

The  physical  properties  of  Dresser  Basalt  and  the  speed  of  the  flame 
jet  was  used  in  equation  (7)  to  calculate  temperature  distributions.  A 
source  strength  of  60  cal /sec  was  chosen  to  produce  subsurface  temperatures 
similar  to  the  experimental  curves  obtained  for  the  flame  jet.  The  peak 
value  of  the  subsurface  temperature  decreases  rapidly  with  a  depth  of  1  mm 
in  Dresser  Basalt  (Fig.  3).  From  the  isothermal  diagrams  made  from  the 
temperature  distribution  (Figs.  4  and  5)  it  can  be  seen  that  for  a  point 
source  moving  on  Dresser  Basalt  the  high  thermal  gradients  all  occur  near 
the  source  and  the  temperature  falls  off  rapidly  within  the  first  milli¬ 
meter  of  the  surface.  While  the  point  source  does  not  represent  the  flame 
jet  accurately,  calculations  indicate  that  the  depth  of  penetration  into 
rock  material  is  very  shallow  and  the  thermal  gradients  are  highest  with¬ 
in  one  millimeter  of  the  surface.  This  is  borne  out  experimentally  when 
the  flame  torch  is  passed  over  a  spallahle  material  and  a  large  number  of 
particles  of  1  mm  and  smaller  are  produced. 
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salt 


Degrees  Centigrade 


Isothermal  Linesfor  Point  Source  -  Dresser  Basalt 


While  the  1  mm  subsurface  curve  is  much  sharper  than  the  experimental 
1  mm  curves  the  2  mm  analytical  and  experimental  curves  have  much  the  same 
shape.  The  similarity  of  the  deeper  curves  indicates  that  once  the  wave 
front  has  reached  a  sufficient  depth  it  has  dispersed  enough  that  it  makes 
little  difference  whether  it  was  produced  by  a  point  or  finite  area  source. 
Temperature  Measurements  -  Previous  Research 

Thirumalai  (Ref.  21)  employed  a  miniature  oxy-hydrogen  torch  as  a 
stationary  heat  source  to  study  the  temperature  distribution.  The  torch 
delivered  a  heat  flux  of  170  x  10^  watt/m^  (150  Btu/ft^sec).  Thermocouples 
were  embedded  from  the  back  side  of  the  slabs  of  rock  at  varying  depths 
from  the  test  surface  to  determine  thermal  gradients  (Table  2).  Sioux  quart 
zite  spalled  readily  and  showed  the  highest  thermal  gradient,  while  charcoal 
granite  spalled  only  slightly  with  partial  melting,  and  basalt  melted  at 
the  experimental  temperature.  The  latter  also  demonstrated  an  increase 
of  thickness  of  the  heat  layer. 

Experimental  Procedure 

A  surface  temperature  resistive  gage  was  bonded  with  ceramic  cement 
to  the  surface  of  a  1 1  x  1'  x  1"  rock  specimen.  Chromel -alumel  thermo¬ 
couples  (.005"  diam)  encased  in  Omegatite  were  placed  in  holes  (2  mm  diem) 
drilled  in  the  rock  slab  from  the  back  side  (Fig.  6).  One  hole  was  drilled 
through  the  surface  and  a  chromel -alumel  thermocouple  placed  in  the  hole 
to  expose  the  thermocouple  to  the  surface.  The  torch  was  passed  over  the 
rock  specimen  at  a  constant  speed  and  the  temperature  relative  to  the 
torch  position  recorded. 

Preliminary  surface  temperature  measurements  were  made  at  a  specified 
point  on  the  rock  surface  for  three  traverse  speeds,  three  offset  distances 
from  the  plane  of  the  sensors,  and  the  torch  distance  of  14  in.  from  the 


27 


CD 

in 


t o 
E 
S- 

<D 

-C 


x:  c 
h-  QJ  £ 
*r  £ 
CD  3  \ 
CD  C_J 
ftj  i- 
i~  CD  o 
CD 
> 

< 


OJVOr'CTiOOr-in 
^rocr^M  oroo^rs 
<\J  cm  i—  cm  ry»  cm  r—  <— 


fi  m  n 

r*—  CM  LT) 
r—  O  O 

vitvii  vii 


o  to 

vo  VO 

V  V 


CO 

<u 


<D 


<U  — ^ 
0>C\| 
ra  .— 
c.  — - 
<u 

>  CD 
<  C 

-a 


C\J 

a> 


-Q 

03 


QJ 

c. 


03 
C-  i 
O  I 
CL  < 


O) 

c 


C  f— 

E  Q- 

fT3  CO 

oJ  in  o 

CL 

»— 

CD  OO 

CTJO 

i. 

E  C 

3  CD 

3  *r- 

■M  C 

£  S- 

03  *r— 

•«—  3 

s- 

X  T3 

CD  3 

ra 

CLQ 

1, 

a; 

E  ** 

3  _J 

E 

s 

*r~  +-> 

x  <a 
ra  aj 

in 

2:  x 

m 

r—  <+- 

a> 

c 

ra  o 

-X 

4-) 

u 

C  -f-> 

*r* 

QJ  C 

XI 

E  QJ 

»- 

•r—  t— 

s-  *° 

TD 

a)  ra 

CD 

a.  s- 

XJ 

X  CD 

u. 

s~ 

aj 

>> 

03 


O 

O 

<d 

oc 


CL  fO 
CD  i- 
Q  4J 
C 
CD 
CJ 


OsJO^^tinrN^cNJ 
CM  LO  O')  07  LO  LO  C^i  CM 
CNJCNJi —  r-  CM  i —  r—  CM 


CD 


N 

-M 

fU 

3 

Cr 

x 

3 

O 

•r— 

00 


(U 


c 

ra 

s- 

o 


fT3 

O 

£ 

ra 

X. 

o 


LO  CO  CM  LO  f'-H  <V)  (\j 

COOOCOvNJNOM 


o  o  o 

LO  CM  CO 

co  co  co 


ro 

in 

03 

ca 

a> 

t/> 

in 

CL) 

S-. 

o 


00  CO  CO 


CO  07  .O 
VII  VII  VII 


J*. 

o 


u 

CD 

</) 


00 


o 

LO 


u 

rU 

in 

in 

<u 

u 

<11 

c 

C7 

ra 

* 


CJ 

+-> 

E 

QJ 

4-> 

E 

■M 

<a 

0) 

-C 

4-> 

-5- 

o 

Vf- 

r~ 

O 

oo  oo 

I—-  r— 

X 

•— 

CO  co 

o 

<a 

> 

A  A 

o 

r  — 

E  a> 

i  M- 

03  c 

|  O 

<4- 

C.  «r— 

O 

4“> 

X 

ra  CD 

3 

u  E 

r“m 

•r» 

M- 

C  in 
ra  3 

4-> 

JC  o 

<a  • 

U  3 

ai  Tiai  c 

_C  j= 

E  T- 

i- 

•r- 

+J 

o 

+-»  r— 

1  c 

in 

C  i — 

o 

c 

<D  ra 

03  o 

CD 

"D  Q.  C 

03CDLOLOLOCOOO 

inLocococoLor*vr'K 

r—  r—  LO  IQ  LO  07  CVJ  CM 


co  co  oo 

OO  CO  CO 


00  co 
co  CO 


o  -M 
c  air-  cn 
v-  >  QJ  C 

4J  in  ^  r^I 
C  <n  f—  r— 
03  <D  <a  ra 
4->  S-  *r-  O. 

<n  cn+->  in 
c  o  u 
o  t  it}  o 
o  a.  a.  z 

»-i  <sj  m  j#- 


23 


-  gradients  calculated  from  the  recorded  rate  of  temperature  ri 


Locations  of  Temperature  Sensors  for  Temperature  Measurements 


rock  (Figs.  7,  8,  and  9)  using  Missouri  Red  Granite.  Temperature  magni¬ 
tude  and  duration  increased  with  smaller  offset  and  decrease  in  table 
speed.  The  block  diagram  (Fig.  10)  illustrates  the  instrumentation  by 
which  surface  temperatures  were  obtained. 

Subsurface  temperature  measurements  were  made  on  the  following  rocks 
at  the  depths  indicated  for  one  traverse  speed  and  three  offset  conditions: 

A.  Dresser  Basalt  1  mm,  2  mm,  3  mm,  4  mm,  5  mm,  6  mm 

B.  Missouri  Red  Granite  1  mm,  ?  mm,  4  mm 

C.  Charcoal  Grey  Granite  1  mm,  2  mm,  4  mm 

D.  Buena  Black  Granite  1  mm,  2  mm,  4  mm 

Typical  curves  for  Dresser  Basalt  (Figs.  11  and  12)  show  that  the  signi¬ 
ficant  temperature  rise  occurs  only  near  the  surface.  Other  recorded 
temperature  curves  appear  in  Appendix  B.  The  block  diagram  of  Figure  13 
illustrates  the  test  setup  used  to  obtain  the  subsurface  temperature  curves. 
The  electronic  equipment  (Fig.  14)  consisted  of  two  Tektronix  X549  dual  trace 
storage  oscilloscopes,  four  ice  reference  junctions,  and  a  specimen  of  rock 
with  four  thermocouples. 

Heat  Source.  The  heat  source  used  in  the  tests  described  in  this  chap¬ 
ter  was  an  oxygen  kerosene  FSJ-6  stone-shaping  flame  jet  torch.  Stoichio¬ 
metric  calculations  (manufacturer's)  give  the  maximum  combustion  energy 
as  450,000  Btu/hr  for  complete  combustion.  The  flame  jet  had  an  orange- 
yellow  flame  signifying  incomplete  combustion. 

Although  the  flow  on  the  kerosene  and  oxygen  was  closely  controlled 
the  temperature  data  indicated  that  there  were  differences  in  the  torch 
temperature  for  different  tests  (Figs.  15  and  16).  The  gas  film  tempera¬ 
tures  obtained  when  the  torch  was  operated  at  standoff  distances  of  8,  10, 
and  14  inches  (Figs.  15  and  16)  show  considerable  scatter-.  The  leveling 
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Surface  Temperature  at  a  Fixed  Point  due  to  a  Moving  Heat  Source  -  (Mo.  Red  Granite) 


Fig.  8.  Surface  Temperature  at  a  Fixed  Point  due  to  a  Moving  Heat  Source  (Mo.  Red  Granite) 


Distance  of  Torch  f rom  Tiermocoupl 


Subsurface  Temperature  Measurements  -  Dresser  Basalt 


Distance  of  Torch  from  Thermocouple  (cm) 
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Distance  of  Torch  from  Thermocouple  (in.) 

Fig.  12.  Subsurface  Temperature  Measurements  -  Dresser  Basalt 


.  Subsurface  Temperature  Measurement  Instrumentation 


14. 


Subsurface  Temperature  Measurement  Test  Setup 
shown  with  Acetylene  Torch. 
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off  of  the  curves  near  zero  offset  was  expected  (Ref.  10).  The  highest 
temperatures  did  not  always  occur  at  zero  offset,  probably  resulting 
from  one  or  more  of  the  following  factors: 

1.  Carbon  builds  up  in  the  combustion  chamber  and  gas  accelerating 
nozzle  due  to  incomplete  combustion  of  the  fuel. 

2.  The  torch  could  not  be  adjusted  to  exactly  the  same  operating 
conditions  each  time  it  was  lighted. 

3.  Slight  misalignment  of  the  torch  at  standoff  distances  of 

10  and  14  inches  would  be  magnified  on  the  surface  of  the  rock. 

4.  The  large  standoff  distances  (10-14  inches)  could  allow  a 
surface  gas  layer  to  build  up  around  the  stagnation  point. 

The  problems  of  alignment  and  large  surface  gas  layers  should  be  mini¬ 
mized  when  the  torch  is  moved  closer  to  the  rock  in  future  tests.  The 
carbon  buildup  in  the  torch  has  been  slowed  by  preventing  unburned  carbon 
from  remaining  in  the  torch  during  shutdown.  Experience  in  operating  the 
torch  plus  repetition  of  some  of  the  tests  reduce  the  error  due  to  im¬ 
proper  adjustment  of  the  torch. 

The  approximate  upper  and  lower  bounds  of  the  surface  temperature  and 
approximate  gradients  within  1  mm  of  the  surface  for  Missouri  Red  Granite 
were  projected  from  measurements  (Fig.  17). 

Numerical  Procedures 

The  maximum  temperatures  for  the  time  histories  were  calculated  for 
varying  offset  and  depth  (Fig.  18).  For  calculations  the  slab  was  divided 
(Fig.  19)  to  represent  strips.  An  average  temperature  of  a  strip  was  cal¬ 
culated  and  the  heat  required  to  raise  the  temperature  of  the  strips  to 
the  average  temperatures  was  determined. 
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A —  Surface  gage  temperature 
□ —  Gas  film  temperature 
O  —  Subsurface  temperatures 


Heat  source  -  Flame 

Speed  -  2.1  in/sec;  5.33  cm/sec 

Distance  -  14  in. 

Offset  -  0  in.;  0  cm 


w 

w 


w 


Depth  (cm) 


Approximate  Temperature  Distribution  -  Missouri  Red  Granite 


*Note:  Gas  Film  Temperatures 

Fig.  13.  Subsurface  and  Gas  Film  Temperatures  (Degrees  -  C) 


The  total  heat  was  then  the  sum  of  the  heat  that  entered  the  strips. 

<Wi ■  j,  "i  ■  j,  eiCjMW  <8> 

where 

C,.  *  specific  heat  at  temperature 
A.j  =  cross-sectional  area  of  the  strip 
T.j  =  maximum  temperature  of  the  strip 
Tq  -  initial  temperature  of  the  slab 
L  =  length  of  the  slab 
p  =  density  of  the  slab 
n  =  number  of  strips 

The  coupling  effect  can  then  be  expressed  as  a  percent  of  the  total 
available  heat  energy  from  the  torch. 

(l  — ^  x  100  =  HEC  (heat  energy  coupling)  (9) 

\  ^torch  / 
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The  curves  for  Dresser  Basalt  were  chosen  to  illustrate  the  calcula¬ 
tion  of  heat  energy  coupling,  for  a  table  speed  of  4.5  in. /sec  and  a  dis¬ 
tance  from  the  torch  to  the  rock  of  14  in.  The  choice  of  Dresser  Basalt 
was  made  because  temperatures  had  been  recorded  at  six  subsurface  depths 
and  some  pertinent  material  properties  could  be  found  in  References  22  and 
23.  The  surface  temperature  measurements  for  basalt  were  not  completed 
but  were  estimated.  Three  cases  were  chosen: 

Case  I.  Temperature  at  surface  assumed  to  be  between  the  gas 
film  temperatures  and  the  1  mm  subsurface  temperature. 

T  =  250°C 

Case  II.  Temperature  of  surface  assumed  to  be  at  the  gas  film 
temperature. 

Case  III.  Temperature  of  surface  assumed  to  be  the  same  as  1  mm 
subsurface  temperature. 

The  density  was  assumed  to  be  constant  for  all  cases  and  equal  to  2.70  g/cc. 
The  heat  input  was  considered  to  be  centered  along  a  line  six  inches  long 
and  moving  with  the  speed  of  4.5  in. /sec  for  a  distance  of  1  cm.  The 
rock  heat  capacity  was  taken  as  960  joules/kg-°K  for  temperatures  below 
100°C  and  as  980  joules/kg-°K  for  temperatures  less  than  160°C.  It  was 
further  assumed  that  all  the  materials  within  the  6  in.  by  1  cm  area 
directly  under  the  source  reached  their  peak  temperatures  simultaneously. 

<Wh  °  11  ■SM  J°u,es 
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For  Case  I. 


Qtotal  =  1068  jou1es 
HEC  =  9.3% 

For  Case  II. 

Qtotal  "  2056  j0u,es 
HEC  =  17.9% 

For  Case  III. 

Qtotal  *  650  j°u'es 
HEC  =  5.7% 

The  low  value  of  5.7%  coupling  efficiency  appears  to  be  the  most 
acceptable,  although  estimates  by  various  investigators  give  values  as 
high  as  7  percent  for  an  optimum  standoff  of  about  6  inches.  For  an  in¬ 
verse  square  relationship  the  coupling  efficiency  for  14  in.  standoff 
would  be  about  2  percent.  Hence,  the  above  approximate  calculations  are 
reasonably  accurate  in  view  of  the  acute  difficulties  involved  in  making 
near  surface  temperature  measurements  in  rock. 


47 


CHAPTER  III 
THERMAL  STRESSES 

Introduction 

The  effects  of  thermal  stresses  induced  by  heat  sources  applied  to 
the  surface  of  rocks  and  other  brittle  materials  are  determined  by  the 
heat  flux,  diffusivity,  conductivity,  thermal  coefficient  of  expansion, 
fractures  and  flaws,  rate  of  travel  of  moving  sources,  strength  of  rock, 
and  related  factors.  One  of  the  important  phenomena  which  has  not  been 
considered  in  the  literature  is  the  effect  of  stress  relief  upon  stress 
distribution,  failure  mechanisms,  heat  transfer,  and  temperature  distri¬ 
bution  after  fractures  are  induced  in  the  rock  surface,  or  after  spalla¬ 
tion  occurs. 

One  of  the  planned  objectives  of  the  theoretical  analysis  was  to  find 
methods  of  numerical  evaluation  for  solutions  of  equations  for  thermal 
stress  distributions  in  semi -infinite  solids  due  to  moving  heat  sources. 
However,  existing  solutions  are  for  point  sources  that  involve  mathe¬ 
matical  discontinuities  which  cannot  be  programmed  in  the  computer  to 
obtain  accurate  calculated  values.  Also,  the  calculation  of  three- 
dimensional  stress  distributions  from  known  or  measured  temperature  distri¬ 
butions  was  found  to  require  so  much  detailed  programming  as  to  be  pro¬ 
be  :  tivo. 

Tht  results  achieved  in  the  analysis  which  was  completed  will  serve 
a  a  is-  .e  *  rr  the  future  development  of  possible  simplified  methods  of 
calcu'«tion  of  tnermal  stresses.  Existing  solutions  for  temperature 
distribution  for  a  moving  heat  source  are  useful  in  determining  the  posi¬ 
tion  of  critical  thermal  gradients.  Although  the  state-of-the-art  of 
quantitative  analysis  of  thermal  shock  is  limited,  concepts  involved  are 
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likewise  useful  in  explaining  observed  phenomena  in  rocks.  Steady  state 
models  employing  the  concept  of  bending  stresses  in  a  plate  to  simulate 
a  heated  surface,  and  a  hemispherical  heat  source  model  are  employed  to 
define  some  of  the  effects  of  both  surface  and  subsurface  thermal  gra¬ 
dients.  In  general,  rock  temperatures  which  cause  spallation  are  rela¬ 
tively  low,  and  the  behavior  of  the  rock  is  considered  to  be  elastic. 

Moving  Heat  Sources 

The  following  describes  the  theoretical  study  of  thermal  stresses  in 
rocks  due  to  a  moving  heat  source.  The  physical  problem  consists  of  ob¬ 
taining  the  stress  distribution  in  a  rock  semi -space  due  to  a  heat  source 
of  known  strength,  w,  moving  with  a  constant  velocity,  v,  along  the  x.i-axis 
on  the  bounding  plane  x3  =  0,  as  shown  in  the  figure  below. 


Fig.  20. 


Point  Source  Moving  on  the  Bounding  Plane  of  a  Semi-infinite  Rock 


In  this  analysis,  the  heat  source  is  assumed  to  be  moving  steadily 
and  continuously  on  the  rock  surface,  and  the  resulting  temperatures 
and  displacements  are  small.  The  dynamic  effects  of  thermal  inertia  can 
then  be  neglected  and  the  analysis  set  within  the  framework  of  the  linear 
theory  of  uncoupled  thermoelasticity.  For  this  case  the  stress  problem 
becomes  explicit,  that  is,  the  equations  for  the  stresses  are  obtained 
from  the  known  temperature  distribution  derived  by  solving  a  well-defined 
heat  conduction  problem. 

Temperature  Distribution.  Rocks  are  complex  materials  with  inherent 
inhomogeneity  and  anisotropy  as  well  as  temperature-dependent  properties. 
The  problem  of  thermal  stress  analysis  in  rocks,  even  under  the  above 
mentioned  assumptions,  is  an  extremely  difficult  one.  An  exact  solution 
of  the  temperature  problem  in  itself  is  a  formidable  task  because  it  in¬ 
volves: 

1.  A  three-dimensional  formulation 

2.  Non homogeneous  and  anisotropic  material 

3.  Temperature-dependent  properties,  and 

4.  Transient  analysis  as  a  result  of  the  motion  of  the  heat  source 

It  should  be  noted  that  the  inclusion  of  each  of  these  factors  in¬ 
volves  a  different  formulation  of  the  problem.  Moreover,  the  problem  is 
not  linear  in  that  the  solution  to  an  equation  involving  two  or  more 

of  the  above  factors  cannot  be  obtained  by  combining  the  solutions  to 
the  problems  involving  each  of  these  factors  separately.  That  is,  there 
does  not  exist  a  general  method  for  the  solution  of  nonlinear  partial 
differential  equations. 

Numerous  analytical  methods  have  been  developed  for  solving  the  non¬ 
linear  equations  of  unsteady  state  heat  conduction,  each  with  its  own 


50 


limitations  and  rather  narrow  field  of  application  (Refs.  24-28).  Many 
of  these  methods  are  attractive  only  theoretically  as  their  practical 
applications  lead  to  cumbersome  mathematical  operations  that  cannot  al¬ 
ways  be  performed  without  very  special  techniques  or  simplifying  assump¬ 
tions  (Refs.  20,  29-33).  In  order  for  one  to  be  able  to  develop  a  work¬ 
able  solution  to  thermal  stress  problems  in  rocks,  simplifying  assumptions 
are  introduced  based  on  such  practical  considerations  as  the  behavior  of 
rock  properties  and  the  main  objective  of  the  research. 

On  a  macroscopic  scale,  the  rock  body  may  be  treated  as  homogeneous 
material  because  of  its  large  mass  and  random  crystal  orientation.  Also, 
due  to  the  manner  of  their  formation,  in  general,  rocks  exhibit  different 
properties  in  directions  parallel  and  perpendicular  to  structure  such,  as 
bedding,  but  for  a  first  approximation  may  be  treated  as  isotropic  solids 
(Ref.  34).  The  main  objective  of  this  phase  of  research  was  to  study  the 
influence  of  the  temperature  dependence  of  the  thermal  and  elastic  proper¬ 
ties  of  rocks  on  the  stresses  produced  by  a  moving  heat  source.  This,  com¬ 
bined  with  the  fact  that  in  the  literature  the  temperature-dependent  rock 
properties  are  almost  invariably  obtained  by  assuming  the  material  to  be 
homogeneous  and  isotropic  justifies  the  using  of  simplified  analysis 
(Refs.  23  and  35). 

Under  these  assumptions,  the  mathematical  equations  for  the  heat  con¬ 
duction  problem  are  as  follows.  The  governing  equation  is: 


3T  -  9T  9T  .  9  9T 

pC^  3T  ~  3x j  3x  +  9x2  3X 

1  2  2 

+  -~  k(T)  |L  +w«(x,-vt,0,  0;  t) 


(10) 
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Boundary  Conditions: 


9T 

1*  (xi,  x2,  0;  t)  =  0,  if  it  is  assumed  that  the 

bounding  plane  x3  =  0  is  thermally  insulated. 
Also,  note  that  T(vt,  0,  0;  t)  -* 

2.  T(R;  t)  =  Tq,  R  -  ®. 

Initial  Condition: 

T(R;  0)  =  T0. 

In  these  equations,  the  following  notations  are  used: 


01) 


X  1  *  )  X3 

R  =  (xf  +  x|  +  xi)* 
t 

T  —  T(Xi ,  x2,  x3;  t) 
To 

PC  (T) 

k(T) 

w 

<$(Xl-Vt,  0,  0;  t) 


=  space  variables,  cm 
=  space  norm  in  the  x..  space 
=  time  variable,  sec 
=  temperature,  °C 

=  initial  or  reference  temperature 

=  volume  heat  capacity  ,  — 

cnr-°C 

=  thermal  conductivity, 

=  rate  of  heat  source. 


02) 


03) 


cal 

sec 


cal 

sec-cm-°C 


=  Dirac  Delta  function  defined  as: 

=  0,  Xi  f  vt 

=  1,  Xi  =  vt,  x2  =  x3  =  0  at  time  =  t  sec 
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Quasi -Stationary  Approximation 

The  mathematical  analysis  of  the  system  of  equations  (10-13)  is  greatly 
simplified  if  it  is  assumed  that  the  state  of  temperature  is  quasi- 
stationary.  Bornefeld  (Ref.  36)  experimentally  established  that  the 
temperature  distribution  in  a  solid  heated  by  a  moving  heat  source,  except 
during  the  initial  transient  period,  remains  stationary  when  referred  to 
a  coordinate  system  moving  with  the  source.  Based  on  this  assumption, 
Rosenthal  (Ref.  20)  •  "ived  his  classic  theory  of  moving  heat  sources  for 
homogeneous,  isotropic  materials  having  constant  thermal  properties. 

The  governing  heat  conduction  equation  for  this  case  is 

V’T  ■  7  It  <,4> 


where 

V2  --  +  :-y-r  +  ,  the  Laplacian ,  and  (15) 

dXi  oaz  oX  3 
l 

k  =  ,  the  diffusivity. 

If  a  new  coordinate  system,  moving  with  the  source  is  introduced, 

then 


Ci  "  *i  -  vt,  C2  =  x2,  C 3  =  x 3  . 

Hence,  the  governing  equation  in  the  moving  system  becomes 
32T  J2T  3ZT  J  31  v  9T 

acT  Wz  acT  <  at  '  <  aci  ' 

For  the  quasi-stationary  temperature  distribution 
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and  the  governing  equation  reduces  to 


92T  52T  ,  3^T  .  y  dT 

3^i  3^2  as!  3Ci  ‘ 

This  equation  is  further  simplified  by  putting 
T  =  Tq  +  exp(^r)  '/•'(Si.  S2,  €3) 

After  substitution,  the  following  equation  is  obtained: 


or  in  the  usual  symbolic  form 


Tho  function  i|;  =  iKS^)  is  to  be  determined  from  the  above  equation  and 
the  boundary  conditions. 

Numerical  Techniques 

The  above  discussion  suggests  that  one  approach  for  the  solution  of 
the  temperature  problem  in  rocks  is  to  oversimplify  the  problem  and  then 
use  one  of  the  applicable  analytical  methods.  An  alternative  approach 
is  to  obtain  the  solution  numerically.  This  approach  has  the  advantage 
that  it  does  not  require  oversimplification.  The  obvious  disadvantage 
is  that  solution  can  he  obtained  for  a  finite  region  and  finite  time 
only  because  the  numerical  treatment  of  the  problem  imposes  restrictions 
on  computing  time  and  storage. 

The  most  common  methods  of  obtaining  a  numerical  solution  are: 

1.  The  finite-element  techniques,  and 

2.  The  finite-difference  techniques 


(17) 

(18) 


The  steep  temperature  gradients  in  rocks  (Ref.  37)  require  very  small- 
sized  elements  (small  grid  size  in  the  case  of  difference  methods).  The 
use  of  the  finite-element  techniques  was  deemed  impracticable  initially, 
due  to  the  large  storage  and  time  requirements  and  nonavailability  of  a 
three-dimensional  dynamic  code.  Effort  has  therefore  been  concentrated 
on  finite-difference  techniques. 

Three  well-known  finite-difference  approximations  are  (Refs.  38  and  39): 

1.  Forward  difference  (explicit  scheme), 

2.  Backward  difference  (implicit  scheme),  and 

3.  Crank-Nicol son  (Ref.  40)  and  other  combination  schemes  (combina¬ 
tion  of  the  first  two  schemes). 

The  explicit  schemes  suffer  because  of  stability  criteria,  hence,  pro¬ 
hibitively  small  time  increments  must  be  used.  The  implicit  methods  lead 
to  a  set  of  simultaneous  equations  which  must  be  solved  at  each  time  step. 
Thus,  these  techniques,  while  otherwise  efficient,  soon  become  impractical 
as  the  number  of  nodes  increases.  The  combination  schemes  have  a  higher 
order  accuracy  compared  to  the  other  two  types  of  approximations,  but 
they  too  suffer  from  the  last  mentioned  limitation  and  cannot  be  effec¬ 
tively  employed  for  most  problems  involving  large  numbers  of  nodes. 

After  an  extensive  survey  of  difference  methods  was  made,  the  alter¬ 
nating  direction  implicit  (ADI)  technique  of  Douglas  (Ref.  41)  (discovered 
independently  by  Brian  [Ref.  42])  was  chosen  for  the  solution  of  the  tem¬ 
perature  problem.  This  technique  was  first  introduced  by  Peaceman  and 
Rachford  (Ref.  43)  and  Douglas  and  Rachford  (Ref.  44).  Being  uncondition¬ 
ally  stable,  as  the  backward-difference  and  the  combination  schemes,  the 
ADI  techniques  require  the  solution  of  many  small  independent  systems  of 
simultaneous  equations  instead  of  one  large  system.  Furthermore,  each 
system  is  tridiagonal  and  a  very  efficient  method  for  solving  the  tri¬ 
diagonal  system  (Thomas  algorithm)  without  iteration  is  well  known  (Ref.  43). 
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Jump  Condition.  Since  the  heat  source  is  acting  on  a  finite  area 
of  the  bounding  plane  x3  =  0,  the  boundary  condition  on  this  plane  in¬ 
volves  a  discontinuity.  In  equation  (10)  the  assumption  that  the  source 
is  acting  at  a  single  point  was  introduced  primarily  to  simplify  the  pre¬ 
liminary  numerical  analysis.  For  this  case,  the  boundary  discontinuity 
gives  rise  to  a  step  function  as  shown  below  by  the  jump  condition  (Eq.22). 

For  programming  convenience,  the  governing  equation  (10)  is  replaced 
by  the  equivalent  system  of  the  homogeneous  equation: 


with  the  continuity  requirement 


I 


T ( x 7  -  vt,  0,  0;  t)  =  T(xt.-  vt,  0,  0;  t) 


(23) 


The  (-)  and  (+)  superscripts  denote,  respectively,  the  Xi-coordinates 
of  the  points  just  to  the  left  and  just  to  the  right  of  the  point  x  =  vt. 

The  presence  of  the  jump  condition  greatly  complicates  the  analysis 
iri  some  cases.  Special  finite-difference  analogs  muct  be  obtained  to 
approximate  this  type  of  jump  condition. 

For  the  investigation  of  the  accuracy  and  consistency  of  the  differ¬ 
ence  analogs  that  were  derived,  it  was  necessary  to  solve  a  simple  prob¬ 
lem  involving  a  similar  jump  condition.  A  two-dimensional  problem  of 
heat  flow  in  a  thin  plate  of  constant  thickness  g,  heated  by  a  source 
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moving  along  one  of  its  edges  was  selected  (Fig.  21).  The  analytical 
solution  for  this  problem,  based  on  the  assumption  of  constant  proper¬ 
ties  and  a  quasi-stationary  state  is  given  by  Rosenthal  (Ref.  20). 


Fig.  21.  Point  Source  Moving  along  an  Edge  of  a  Thin  Plate 


The  system  of  governing  equations  for  this  problem  involves  the 
governing  homogeneous  equation 


/  32T  +  a2T\ 


the  jump  condition 


-k  |y  (x-  -  vt,  0;  t)  +  k  |y  (x+  -  vt,  0;  t)  =  w 
and  the  continuity  requirement 


T(x~  -  vt,  0;  t)  =  T(x+  -  vt,  0;  t) 
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Finite-Difference  Approximation.  In  the  figure  below  let  the  source 

A  v 

be  at  node  5  at  time  t.  Let  the  grid  ratio  =  v  so  that  in  one  time 
step  the  source  advances  by  one  space  step  in  the  x-direction. 

w 


_ J 

_ 

6 

5 

7 

3 

4 

1 

2 

y 

Fig.  22.  Finite-Difference  Grid  for  the  Jump  Condition  Approximati 


Let  Ax  =  Ay  =  h. 

The  use  of  the  three  point  forward  difference  formulas  yields 

!R-k<-3  '"'T'-T.l 

and 


3T 

3y 


1 

2h 


(-3  T6  +  4  T3 


Ti). 


Since  the  source  has  moved  from  node  6  to  node  5,  the  continuity 
condition  at  node  5  is  approximated  by 

T6  =  T5 

and  substitution  of  the  last  three  equations  in  the  jump  condition 
(25)  yields: 
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f-  +  T7  +  y  (Ta  -  T„)  +  y  (T2  -  T,) 


In  the  case  of  a  stationary  jump,  the  effect  of  the  discontinuity 
gradually  dies  out  with  each  iteration  and  the  number  of  minimum  itera¬ 
tions  required  to  smooth  out  the  discontinuity  is  proportional  to  the 
grid  size  (Ref.  38).  The  problem  of  a  moving  heat  source  involves  a 
"moving  discontinuity".  To  smooth  out  this  type  of  discontinuity,  a 
very  large  number  of  iterations  are  required  at  each  time  increment.  The 
approximation  (27)  was  found  to  be  impractical  for  the  problem  under  con¬ 
sideration,  as  the  results  obtained  by  using  a  practically  allowable 
smallest  value  of  h  and  a  large  number  of  iterations  were  unacceptable. 

Other  approximations  were  obtained  by  replacing  the  jump  condition 
by  a  flux  boundary  condition.  One  approximation  of  this  type  was 


'k  ay ! 5  “  w  ’ 


and  the  other,  based  on  an  average  value  was 


1/2  (  -k  iyl  5  “k  |y|  6  }  =  w  . 


Again,  in  both  of  these  cases,  the  same  computational  conditions  as 
for  equation  (27)  were  found  to  be  true. 

Stress  Distribution.  As  mentioned  in  the  previous  section,  the 
stress  problem  in  the  case  of  uncoupled  theory  of  linear  thermoelasticity 
is  an  explicit  one.  The  method  proposed  for  the  analysis  of  the  thermo- 
elastic  problem  under  consideration  was  to  obtain  the  temperatures  by 
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solving  the  heat  conduction  problem  and  then  use  these  values  to  solve 
the  stress  problem  numerically  by  using  one  of  the  available  finite  ele¬ 
ment  codes.  This  method  of  analysis  was  chosen  after  a  detailed  survey 
of  existing  techniques  for  the  solution  of  nonlinear  problems  of  transient 
thermal  stresses  (Refs.  45-50).  It  is  important  to  note  that  the  term 
"nonlinear"  here  refers  to  the  elastic  solution  of  problems  having  vari¬ 
able  material  properties  and  not  to  the  conventional  usage  for  plastic 
analysis. 

The  problem  of  a  semi-infinite  elastic  solid  with  a  plane  boundary, 
commonly  known  as  the  problem  of  Bossinesq  and  Cerruti,  has  been  the  ob¬ 
ject  of  study  by  numerous  researchers.  On  the  other  hand,  few  problems 
are  solved  with  nonlinear  material  properties.  Some  of  these  are  given 
in  References  51-54.  Analyses  given  in  these  references  are  limited  to 
one-or  two-dimensional  problems  only. 

Nowacki's  Solution.  For  the  purpose  of  preliminary  analysis,  solu¬ 
tion  to  a  semi-infinite  rock  semi-space  problem  was  attempted  assuming 
the  material  to  be  homogeneous,  and  isotropic  with  temperature  independent 
properties.  A  closed  form  solution  for  this  problem  is  given  by  Nowacki 
(Ref.  49)  and  is  described  in  Appendix  C. 

From  the  solution,  given  by  equation  (C-2)  it  is  seen  that  the  tempera¬ 
ture  at  the  origin,  R  =  0  tends  to  infinity.  Similarly,  the  stresses 
in  the  infinite  space  become  infinite  at  R  =  0.  The  expression  (C-4)  for 
the  stresses  o  involves  an  improper  double  integral  which  is  too  com¬ 
plicated  to  be  performed  analytically.  An  attempt  to  evaluate  this  inte¬ 
gral  numerically  using  Simpson's  rule  together  with  Newton's  3/8  formula 
(Ref.  55)  was  unsuccessful  as  the  function  p ( n i ,  nz,  0)  =  -o33(ni,  nz >  0) 
becomes  singular  at  the  origin,  R  =  0. 
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As  an  alternative,  the  following  approach  was  taken.  From  equation 
(C-2)  temperatures  were  obtained  for  different  values  of  R  =  R(£ls  £2,  £3)* 
omitting  the  pole  R  =  0  where  T  -*■  An  attempt  was  made  to  use  these 
temperature  values  explicitly  in  the  three-dimensional  finite  element  code 
ELAS  (Ref.  56).  The  steep  temperature  gradients  in  the  vicinity  of  the 
source  required  extremely  small-sized  elements,  0.01  cm-side  cubes,  for  the 
problem  attempted.  A  sufficiently  large  volume  must  be  considered  in  order 
to  accurately  approximate  the  boundary  conditions  of  zero-stress  at  points 
far  removed  from  the  heat  source.  Thus,  a  grid  consisting  of  extremely 
small -sized  elements  surrounding  the  source  and  progressively  increasing¬ 
sized  elements  away  from  the  source  is  required.  Setting  up  of  such  a 
grid  in  three  dimensions  poses  serious  problems.  The  total  number  of  nodes 
and  elements  becomes  prohibiti vely  large  even  when  only  a  small  block  is 
to  be  analyzed.  This  approach  thus  involved  serious  computer  storage  and 
time  considerations  and  was  deemed  impracticable  in  view  of  the  available 
computer  facilities. 

Approximate  Models 

Other  thermal  stress  models  offer  a  basis  for  indicating  some  of  the 
important  phenomena  involved  in  spallation  such  as  temperature  gradient 
for  a  moving  source,  the  plate  bending  model  and  a  hemispherical  source 
model,  the  latter  two  being  steady  state  in  character.  Some  of  the  con¬ 
cepts  of  thermal  shock  appear  to  have  application. 

Temperature  Distribution  -  Moving  Source.  Rosenthal  (Ref.  20)  has, 
in  addition  to  equations  for  two-dimensional  heat  flow  for  a  quasi¬ 
stationary  state,  developed  a  three-dimensional  equation  which  is  useful 
in  an  analysis  of  temperature  effects  of  moving  heat  sources  on  rock  sur¬ 
faces  . 
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The  expression  for  three-dimensional  heat  flow  may  be  obtained  from 


equation  (10)  to  yield  in  polar  coordinates 


d fi  + 

dR2 


where 


2_  d£  /  v  \ 
R  dR  “  \2k  / 


*s0, 


R  =  (C2+y2+z2)'5 

S  =  x  -  vt  . 


Also, 


hence, 


d2ijj  +  2_  dij)_  _  1_  d2Ri|> 
dR2  R  dR  R  dR2 


d  2RiT>  /  y\ 
dR2  \2k  ) 

The  boundary  conditions  are 


Ri|/  =  0 


-►  0  for  £  -*+'*> 


gy  -  0  f  or  y  -  ± 


gy  0  for  z 


The  temperature  distribution  for  an  infinitely  thick  plate  is  then 


given  by 


T  -  \  -  w  H-fc (e+R)) 


For  a  point  source  the  temperature  becomes  infinite  at  R  =  0. 
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Rosenthal  (Ref.  20)  employed  equation  (28)  to  describe  the  tempera¬ 
ture  distribution  in  thick  plates  subjected  to  the  heat  generated  in 
mild  steel  by  the  arc  of  a  welding  electrode  (Fig.  2). 

The  characteristics  of  the  temperature  distribution  are  given  (Ref.  20) 
as  follows: 

1.  The  temperature  rise  is  much  steeper  than  the  fall,  creating  a 
very  high  thermal  gradient  in  the  surface  of  the  material  in 
front  of  the  source. 

2.  Because  heat  does  not  propagate  instantly  the  temperature  maxi¬ 
mum  envelope  lags  behind  the  heat  point,  the  lag  being  a  func¬ 
tion  of  source  speed  and  diffusivity. 

3.  The  curve  n-n  (Fig.  2)  separates  the  areas  of  rising  and  fall¬ 
ing  temperatures. 

Thermal  Stress-Plate  Bending  Model.  Thirumalai  (Ref.  37)  treats  a 
thin  heated  layer  of  rock  on  the  surface  of  unheated  rock  as  a  thermal 
inclusion  in  the  form  of  a  plate.  The  boundary  stresses  of  the  plate 
are: 


where 


°>oc  *  °yy  ‘  -EaaT/O-v) 


zz 


=  0 


(33) 


E  =  Young 's  modulus 
a  =  coefficient  of  thermal  expansion 
AT  =  temperature  difference  per  unit  depth 
v  =  Poisson's  ratio 
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The  temperature  difference  and  restraint  of  the  plate  cause  a  bending 
stress,  Mi,  per  unit  length: 


M,  .  f  oxx  ZdZ  = 


/  otATzdz 


It  may  be  shown  (Ref.  37)  for  the  assumptions  made  that  for  a  given  rock, 
failure  by  thermal  spalling  depends  upon  the  associated  shear  strain. 

Thermal  Stress-Hemispherical  Source.  Both  the  thermal  shock  and 
thermal  plate  bending  models  attribute  failure  in  spallation  to  the  stress 
accompanying  shear  strain.  However,  each  assumes  uniform  heating  over  a 
surface  and  a  constant  temperature  gradient  to  a  specified  depth.  While 
a  moving  source  and  three-dimensional  models  are  desirable,  closed  solu¬ 
tions  have  been  obtained  only  for  point  sources  and  these  have  not  proven 
amenable  to  computer  analysis  because  of  mathematical  discontinuities. 

A  steady  state  hemispherical  source  model  solved  by  Sternberg  and 
McDowell  (Ref.  57)  sheds  some  further  light  upon  the  effects  of  source 
boundaries,  stress  gradients  on  the  surface  and  at  depth. 

Thus,  for  a  hemispherical  temperature  distribution  over  a  circular 
region  of  exposure  0  <  r  <  a,  z  =  0,  the  heat  flux  is 


f(Q)  =  c 


[’  -(5J J 


The  solution  to  the  steady-state  thermoelastic  problem  described 
is  found  as  follows.  Let  Ai  be  the  positive  root  of  the  equation 


_ +  ZLS  1 

a2  +  A  A 
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(37) 


ifr  :  ^i«fcvssaa^ 


and  define  n  by 


Here,  3  *  ~  G(l+v)  and  G  Is  the  shear  modulus  of  the  material.  In 
equation  (38),  ur  and  uz  denote  the  displacements  in  the  r  and  z  directions 
and  arr  and  denote  the  radial  and  tangential  stresses,  respectively.* 


*  Note:  A  doctoral  dissertation  titled  "Thermal  Fracturing  of  Hard  Crystal¬ 
line  Rock"  by  P.J.  Lauriello,  Rutgers  University,  October  1971, 
and  presented  as  a  paper  at  the  14th  Rock  Mechanics  Symposium, 
discusses  the  thermal  stress  distribution  in  a  semi-infinite  region 
due  to  radiation  of  a  finite  spot  on  the  surface.  The  source  is 
stationary;  however,  it  appears  that  the  solution  could  be  modified 
to  study  the  effect  of  moving  sources.  Lauriello  has  also  modified 
his  solution  for  localized  dielectric  heating  at  the  semi-infinite 
region  and  for  convective  heating  of  a  finite  spot  on  the  surface. 
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The  thermal  gradients  (Fig.  23)*  are  greatest  at  the  surface  of 
the  material  near  the  edge  of  the  heat  source.  The  radial  and  tangential 
stresses  are  likewise  of  highest  magnitude  in  the  same  locality.  It  is 
of  interest  to  note  that  the  tangential  stresses  become  compressive  at 
positions  only  a  short  distance  under  the  heat,  source. 

The  radial  and  tangential  stresses  at  the  surface  are  principal 
stresses  and  the  shear  stress,  equal  to  one-half  their  algebraic  differ¬ 
ence,  is  greatest  at  a  point  just  under  the  edge  of  the  heat  source.  As¬ 
suming  that  the  stress  pattern  is  somewhat  similar  for  a  moving  heat  source, 
both  high  tensile  and  high  shear  stresses  exist  behind  the  heat  front,  parti 
cularly  if  the  rise  in  temperature  is  rapid.  Thus,  it  would  be  reasonable 
to  hypothesize  that  spallation  is  due  to  both  tensile  and  shear  failure, 
with  most  of  the  cracks  originating  at  or  near  the  surface. 

Thermal  Shock.  Thermal  shock  describes  the  phenomena  associated 
with  sudden  transient  changes  in  temperature  of  a  body,  as  differentiated 
from  conditions  of  static  or  slowly  changing  thermal  stress.  No  defini¬ 
tion  appears  to  be  given  in  the  literature  of  the  rate  of  stress  change 
which  describes  the  state  of  thermal  shock.  However,  an  index  of  merit 
is  used  to  indicate  the  relative  resistance  to  thermal  shock  of  a  material. 
It  is  of  particular  interest  in  conjunction  with  sudden  heating  or  cool¬ 
ing  of  brittle  materials  because  of  cracking,  spallation  and  failure. 

*  These  results  were  obtained  using  the  following  values  for  Dresser  Basalt; 
a  =  3  cm 
c  =  600° c 

,«  -  ■]  n-4  -  - 

“  10  cm-°C 

v  =  0.24 
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Isothermal  Lines,  °C 
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Normalized  Radial  Stress,  _!j_ 

G 
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Fig.  23.  Temperature  St res  .  Distribution  -  Hemisf 
Steady  State 
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The  rates  of  traverse  employed  in  the  experimental  work  in  this 
project  were  as  high  as  4.5  in. /sec  with  high  flame  temperatures  and  rela¬ 
tively  high  surface  film  temperatures.  Also,  if  water  is  used  to  cool 
the  rock  after  the  passage  of  the  flame  jet,  the  rock  surface  experiences 
a  transient  thermal  shock. 

If  one  surface  of  a  flat  plate  of  uniform  temperature  is  suddenly 
immersed  in  a  medium  of  lower  temperature,  the  following  equation  can  be 
written  for  the  stress  in  a  plate  with  elastic  temperature-independent 
properties  (Ref.  58): 


(39) 

'  0 

where  a*  is  the  ratio  of  the  stress  developed  to  that  which  would  be  de¬ 
veloped  if  the  thermal  expansion  were  unrestrained.  Also, 


where 

a  =  stress 
v  =  Poisson's  ratio 
E  =  Young  '  s  modulus 
a  =  coefficient  of  thermal  expansion 
T  =  ave.  temp,  across  date  thickness 
T  =  temperature  at  a  point 
Tq  -  ambient  temperature  of  plate 


(40) 


For  thin  plates  it  has  been  found  that  the  surface  stress  is  strongly 
dependent  upon  the  heat  transfer  coefficient  h  and  the  heat  conductivity  k, 
the  stress  being  proportional  to  h/k.  For  a  plate  of  a  given  thickness 
it  can  be  shown  that  the  maximum  shock  temperature  which  it  can  withstand 
is  proportional  to  kn^/Eoh,  where  is  the  breaking  stress.  The  tests 
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to  determine  the  relative  shock  merits  of  materials  consisted  of  subject¬ 
ing  discs  2  in.  diam  by  1/4  in.  thick  to  repeated  thermal  shock  cycles 
until  they  failed,  i.e.,  they  were  heated  in  a  furnace  and  quenched  in 
water. 


For  thick  plates  the  surface  temperature  can  be  changed  long  before 
the  temperature  of  the  inner  part  of  the  body.  For  heat  immersion  the 
surface  layers  cannot  expand  and  under  these  conditions  of  constraint 
the  surface  stress  developed  is  EaTQ/(l-v)  which  is  relatively  independent 
of  the  heat  transfer  coefficient  or  conductivity  of  the  material,  particu¬ 
larly  if  the  latter  is  small.  For  poor  heat  conductors  such  as  rock  only 
the  material  near  the  surface  is  heated  and  the  high  compressive  stresses 
are  limited  to  only  the  surface  layers.  The  better  the  conductivity  the 
greater  the  depth  of  material  that  is  under  compressive  stresses. 

The  above  failure  mechanisms  are  valid  only  for  the  rupture  of  plates, 
which  is  caused  by  the  differential  expansion  of  the  inside  and  outside 
layers  of  the  plate.  The  mechanism  of  spallation  of  the  surface  of  a 
thick  solid  involves  a  different  combination  of  parameters.  That  is, 
in  rock  spallation  a  maximum  size  of  chip  or  fragment  is  desireo.  It 
should  be  noted  that  in  a  heating  process  the  whole  surface  of  a  solid 
is  in  compression,  and  spallation  occurs  by  means  of  shear.  On  the  other 
hand,  when  the  surface  of  a  body  is  cooled,  it  fails  either  in  tension  or 
shear.  For  rocks  the  latter  process  appears  to  be  advantageous  because 
the  tensile  strength  of  rock  is  very  low. 

In  a  study  of  the  spallation  of  ceramic  materials  used  as  liners  for 
rocket  nozzles  and  surfaces  in  similar  environments,  Newman  (Ref.  59) 
points  out  that  when  severe  gradients  are  present  thermal  stresses  arise 
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with  boundary  conditions 


eU.O)  =  gq 

(H) ^  -  o  - » 

The  solution  for  (42)  is 


(42) 


where 


sc  +  e 
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The  shear  strain  may  be  approximated  as  follows  (Ref.  60): 


Fig.  24.  Diagram  Showing  Shear  Strain  Developed  by  a  Temperature 
Gradient  (after  Norton,  Ref.  60) 


The  surfaces  ab  and  a 1 b '  (Fig.  24)  are  isothermal  planes  parallel 
to  the  surface  at  temperatures  T  and  T  +  6T. 

s  ■  S“  § 

where  S  is  a  characteristic  size  of  an  element  of  brittle  material. 

It  is  recognized  that  the  above  equation  serves  largely  as  an  il¬ 
lustrative  method,  and  is  used  here  to  give  only  a  first  approximation 
for  the  stresses  occurring  at  the  heat  front  in  the  rock  surface. 

Using  this  model,  Newman  points  out  that  a  plot  (Fig.  25)  of 
f.  vs.  for  p.  -  100  as  a  function  of  t,  shows  that  the  gradient 

which  will  influence  failure  is  given  by  the  envelope  to  the  individual 
stress  curves.  This  model  shows  that  for  small  t  the  heat  penetration 
is  shallow,  but  the  stress  gradient  is  high.  If  one  assumes  a  minimum 


72 


Reduced  Temperature  Gradient 


Fig.  25. 


Reduced  Temperature  Gradient  as  a  Function  of  Both 
Reduced  Distance  and  Reduced  Time  (after  Newman,  Ref.  60) 


stress  gradient  below  which  fracture  will  not  occur,  then  few  fractures 
originate  at  any  appreciable  distance  below  the  surface.  Likewise,  the 
thermal  strain  energy  required  to  maintain  crack  propagation  diminishes 
rapidly  with  depth. 

The  value  of  strain  along  the  envelope,  which  influences  failure, 
is  g i ven  by 


The  probability  of  crack  initiation  is  assumed  to  follow  the  statisti¬ 
cal  model  developed  by  Weibull  (61),  also  called  the  "risk  of  rupture 
criterion"  for  thermal  shock: 


R  =  /  f(oi)  dVi  (44) 

where  o •  is  the  stress  in  a  volume  element  dVi  and  the  integration  is  per¬ 
formed  for  the  whole  volume  under  stress. 

This  may  also  be  expressed  in  terms  of  strain  (Ref.  59)  for  the  prob¬ 
ability  of  a  single  crack  in  the  body: 


F(e)  =  1  -exp 


-A  /  n(e)  dx 


(45) 


L  0  J 

where  is  the  maximum  crack  depth  and  n(e)  is  the  Weibull  material  func¬ 
tion.  The  probability  that  second  and  third  cracks  will  form  to  release 
a  fragment  is  the  probability  that  they  are  near  enough  to  intersect  the 
other  cracks.  Employing  the  Weibull  distribution  n(c)  =  Km  Newman  shows 
that  the  depth  of  crack  to  be  expected  is  dependent  upon  the  rate  of  heat¬ 
ing,  being  lower  for  higher  heating  rates  (Fig.  26)  where 
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For  useful  application  the  Wei  bull  function  must  be  known,  and  it 
has  been  found  that  only  one-dimensional  models  can  be  scaled  and  this 
scaling  is  restricted  to  small  temperature  gradients.  None  of  the  analy¬ 
ses  given  by  Newman  and  others  consider  the  effect  of  stress  relieving. 

For  example,  some  ceramics  are  more  stress  resistant  where  flaws  have 
been  del iverutely  introduced.  However,  small  local  flaws  may  contri¬ 
bute  to  spallation. 

Crack  Propagation 

Sih  and  Liebowitz  (Ref.  62)  show  for  a  condition  of  steady  heat  flow 
that  when  there  is  a  crack  in  a  solid  whose  thermal  conductivity  is  dif¬ 
ferent  from  the  surrounding  material,  the  local  temperature  gradient  is 
intensified.  The  increase  in  temperature  gradient  is  accompanied  by  in¬ 
creased  thermal  stresses  which  may  cause  propagation  of  imperfections. 

It  has  been  shown  further  that  for  steady  state  temperature  distributions 
the  stresses  at  the  tips  of  Griffith  cracks  are  the  same  as  those  induced 
by  mechanical  forces  as  described  by  the  l//r  (r=crack  tip  radius)  singu¬ 
larity.  Hence,  stress  intensity  factors  and  crack  prediction  caused  by 
heat  flow  may  be  analyzed  in  conformance  with  the  Griff ith-Irwin  theory 
of  isothermal  crack  growth  and  propagation. 

An  energy  balance  equation  can  be  set  up  for  fracture  problems  where 
inertia  is  not  negligible  (Ref.  63).  If  it  can  be  assumed  that  the  energy 
dissipation  around  the  crack  tip  is  small,  the  energy  required  for  such 
dissipation  is  given  by: 

u  -  v  -  T  =  tc  (46) 

where 

u  =  time  rate  of  external  work 

v  =  time  rate  of  elastic  energy 

T  =  time  rate  of  kinetic  en<  rgy  for  a  small  region  at  the  crack  tip 

I.  =  rate  of  crack  closure  energy 

c  J 


For  application  to  fracture  propagation  caused  by  thermal  stress  the 
external  work  is  furnished  by  the  heat  transferred  to  the  rock,  part  of 

which  is  converted  to  elastic  energy.  In  crack  formation  there  must  be 

sufficient  thermal  and  elastic  energy  present  to  initiate  a  crack  and  to 
cause  it  to  propagate. 

If  D  is  the  rate  of  dissipative  energy,  the  Griffith  criterion  states 
that  the  crack  will  propagate  in  the  direction  where  the  E  /D  ratio  is  a 
maximum.  While  the  initial  direction  of  a  fracture  from  a  notch  is  deter¬ 
mined  by  several  factors,  the  plane  of  maximum  tensile  or  shear  stress  may 

be  the  plane  of  fracture,  which  may  be  applicable  in  deep  kerfing. 

Strain  Failure.  Thirumalai  (Ref.  37)  notes  that  the  shear  strain  of 
rock  at  failure  may  be  taken  as  the  sum  of  the  elastic  eg,  and  inelastic, 
ei ,  strains: 

ef  "  ee  +  ei  '47^ 

and  the  degree  of  proportionality  of  elastic  strain  is  given  by 

ep  =  100  ee^Ee  +  £i ^  (4S) 

Measurements  (Ref.  37)  of  shear  strain  parameters  were  performed  on  a 
punch  type  shear  apparatus  and  the  inelasticity  was  determined  by  stress 
recycling.  The  shear  strength  of  quartzite  and  basalt  (Table  3)  were 
found  to  be  the  highest,  and  the  shear  strain  to  cause  failure  in  the 
basalt  was  higher  than  that  for  the  granites.  This  factor  plus  the  low 
thermal  gradients  occurring  in  basalt  partially  explain  its  low  spall- 
ability. 
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Table  3. 


Total  Strain  at  Shear  Failure  and  Degree  of 
Elasticity  for  Three  Test  Rocks* 


Sioux  Quartzite 


Specimen 

number 

Shear  Strength 
psi  N/mm2 

ef 

ee 

percint 

1 

11,880 

83.4 

0.0287 

2 

10,897 

76.5 

.0287 

3 

10,897 

76.5 

.0325 

0.0239 

74 

4 

7,948 

55.8 

.0287 

5 

11,866 

83.3 

.0315 

.0239 

76 

6 

9,914 

69.6 

.0287 

.0210 

74 

7 

12,863 

90.3 

.0315 

.0249 

79 

Mean  value. , 

.0300 

76 

Standard 

deviation . . 

.0017 

2.36 

Charcoal  Granite 

1 

7,550 

53.0 

0.0315 

0.0216 

69 

2 

0,017 

63.3 

.0315 

.0226 

69 

3 

7,749 

54.4 

.0296 

.0195 

66 

4 

7,265 

51.0 

.0306 

.0210 

69 

5 

7,749 

54.4 

.0315 

.0210 

67 

Mean  value. . 

.0309 

68 

Standard 

deviation. . 

.0008 

1.41 

Dresser  Basalt 

1 

9,117 

64.0 

0.0363 

2 

12,849 

90.2 

.0363 

0.0277 

76 

3 

9,900 

69.5 

.0363 

.0268 

74 

4 

13,381 

97.1 

.0344 

.0258 

75 

5 

12,849 

90.2 

.0344 

.0254 

74 

6 

9,900 

69.5 

.0344 

.0254 

74 

Mean  value. . . . 
Standard 

deviation.  - 


.0354  75 

.0010  0.89 


*  The  tests  were  run  without  stress  cycling. 
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CHAPTER  IV 


MECHANICAL  FRAGMENTATION 


Introduction 

Analytical  studies  of  mechanical  fragmentation  of  rocks  for  different 
indentors  have  been  made  by  using  the  slip-line  method.  Four  indentors, 
two  yield  criteria  and  three  interface  conditions  between  rock  and 
indentor  were  studied.  They  are 

A.  Indentors 

1.  Indentor  No.  1  -  Wedge  indentor 

2.  Indentor  No.  2  -  Indentor  with  blunt  point 

3.  Indentor  No.  3  -  Cylindrical  indentor 

4.  Indentor  No.  4  -  Indentor  with  round  point 

B.  Yield  Criteria 

1.  Maximum  shear  theory 

2.  Internal  friction  theory  (Coulomb's  yield  criterion) 

C.  Interface  Conditions 

1 .  Smooth  surface 

2.  Rough  surface 

3.  Frictional  surface 

General  Theory 

The  slip-line  theory  was  developed  to  analyze  plane  strain  deformation 
of  a  rigid  perfectly  plastic  isotropic  solid  (Ref. 64).  The  equilibrium 
equations  for  plane  plastic  flow  are 


=  0 


ax  ay  u 
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(49) 


For  the  maximum  shear  theory,  let 

oxx  =  P  +  K  cos20 

°yy  =  p  -  K  COS20  (50) 

t  =  K  sin2G 

where  p  ■  1/2^°xx+ayy)  is  the  hydrostatic  part  of  the  stress  tensor  and 
e  is  the  angle  measured  from  the  x-axis  to  the  direction  of  the  algebraically 
maximum  principal  stress  shown  in  (Fig.  27). 

Substitution  of  equation  (50)  in  equation  (49)  gives  partial  dif¬ 
ferential  equations  of  the  hyperbolic  type  which  can  be  solved  by  the 
method  of  characteristics  (Ref.  64,65).  The  characteristic  (slip-line) 
directions  are 

3x  =  tan<6±?)  (51) 

which  show  that  the  two  families  of  slip-lines,  a  arid  6,  form  an  orthogonal 
network.  Along  the  slip-lines 

P  ±  2K0  =  constants  (52) 

Similarly,  for  the  internal  friction  theory  a  stress  circle  is  tangent 
to  the  envelope  and  the  governing  stress  equations  are  hyperbolic  (Fig.  ?8). 
The  direction  of  the  slip-lines  is  given  by  the  orientation  of  the  planes 
specified  by  the  points  of  contact  of  the  stress  circle  with  the  yield 
envelope.  Two  families  of  slip-lines  exist  which  are  inclined  at  angles 
of  i  (J-  -  £■)  to  the  direction  of  the  algebraically  maximum  principal  stress. 
Along  the  slip-lines 


3^-  =  tan(e  ±  p),  where  u  =  jj-  -  £ 


(53) 


and 


j  cot  <p  ln(o)  ±  6  =  constants  (54) 

where 

°  =  ^cxx+0yy)  +  K  cot  ♦ 

In  equations  (52)  and  (54),  o  is  a  function  of  6.  It  is  necessary 
to  know  how  0  varies  along  the  slip-lines  in  order  to  obtain  o.  The 
solutions  which  assume  that  0  is  constant  along  one  or  both  families 
of  slip-lines  are  frequently  applied  to  this  problem  (Ref.  66,67)  and  are 
used  in  the  following  analyses. 

Slip-line  Fields.  Two  possible  slip-line  fields,  as  proposed  by 
Prandtl  and  Hill  (64)  are  shown  in  Fig.  29.  Hill  takes  into  account 
the  surface  displacement  next  to  the  wedge.  The  position  of  the  dis¬ 
placed  surface  is  determined  based  on  the  assumption  of  constant  volume. 
Hill's  proposed  slip-line  field  was  applied  in  this  study. 

The  proposed  slip-line  field  (Fig.  29)  applied  to  the  smooth  surface 
of  wedge  i ndentors , can  be  divided  into  three  regions.  In  regions  I  and  III, 
the  two  families  of  slip-lines  are  straight  lines  and  intersect  at  the 
constant  angle  (tj-  -  |-)*,  i.e.,  G  is  constant  and  the  stresses  are  constants 
throughout  the  region.  They  are  the  "constant  regions".  In  region  II, 

0  is  constant  along  one  family  of  slip-lines,  the  other  family  is  composed 
of  circles  or  exponential  spirals,  which  is  the  "centered  fan"  or  "radial 
shear  region".  Every  proposed  slip-line  field  used  in  this  study  will  be 
the  combination  of  these  three  basic  regions. 

A-  TT 

For  maximum  shear  theory  the  angle  is  j  ,  i.e.,  <j>  ~  0. 
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Fig.  27.  Mohr  Diagram  for  Maximum  Shear  Theory 


Fig.  28.  Mohr  Diagram  for  Internal  Friction  Theory 


Fig.  29.  Two  Possible  Slip-line  Fields  Proposed  by  Prandtl  and  Hill 
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Interface  Conditions.  These  are  considered  for  three  conditions. 


(A)  Smooth  Surface: 

The  proposed  slip-line  fields  for  the  smooth  surface  for 
indentors  are  shown  in  Figures  30(a  &  c)  and  31  (a  &  c),  (at  n  =  45°).  In 
this  case,  the  algebraically  maximum  stress  acts  normal  to  the  bit-rock 
interface  and  is  assumed  to  be  constant  along  the  length  of  the  edge. 

The  stress  at  the  interface  can  be  expressed  as: 


ap  =  2  K(1  +  4>)  =  normal  stress 

(55a) 

rn  =  0  =  shear  stress 

for  the  maximum  shear  theory  and  as 

°n  =  2  K 's'ln  'l^^^  exp(2iJ'tan<j>)  -  K  cot  $ 

(55b) 

Tn  =  0 


for  the  internal  friction  theory 

where  ij >  is  the  direction  change  of  one  slip-line  from  the  interface  to  the 
rock  surface. 

(B)  Rough  Surface: 

The  proposed  slip-line  fields  for  a  rough  surface  are  shown  in 
Figures  30(b  &  d),  31 (b  &  d),  and  39.  Since  two  families  of  slip-lines 
always  intersect  at  the  same  angle  (tt/4  -  <p/2)  or  (rr/4  +  <j>/2),  if  8  is 
greater  than  (tt/4  -  <f)/2),  then  a  "false  nose"  region  will  be  created.  The 
"false  nose"  region  is  a  constant  region  and  moves  downward  with  the  bit. 
If  8  is  less  than  tt/4  -  <j>/2,  then  the  interface  coincides  with  a  slip- 
line  of  the  shear  region.  The  stress  at  the  interface  (or  the  edge  of  the 
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false  nose,  if  6  is  greater  than  ir/4  -  <t>/2)  can  be  expressed  as; 


on  =  K  (1  +  2*) 
Tn  -  K 


(56a) 


for  the  maximum  shear  theory  and 


an  =  K  cot  4>  [(1  +  sin4>)  exp(2ijjtan<f>)  -  1] 
t  *  (1  +  sin<j>)  exp(2i^tan4)) 


(56b) 


for  the  internal  friction  theory. 


(C)  Frictional  Surface: 

In  this  case,  the  slip-line  field  is  between  the  smooth  and  rough 
conditions.  Assume  n  is  constant  along  the  indentor  surface.  From  the 
smooth  to  the  rough  condition,  n  will  change  from  (tt/4  +  <j>/2)  to  zero.  In 
proposed  fields  for  the  different  indentors  (Fig.  30(a  &  b)  and  Fig.  31 (a  &  b)) 
if  the  3  angle  is  greater  than  (tt/4  -  <j>/2) ,  then  the  frictional  surface  has 
the  same  slip-line  field  as  the  rough  surface. 

The  stresses  at  the  interface  at  3  <  (rr/4  -  <p/ 2)  will  be: 


o  =  K  (1  +  2ip  +  sin2n) 
i  =  K  cos  2n 


(57a) 


for  the  maximum  shear  theory  and 


=  --"coif"^  lsik  '  Sin  (*  ~  2n)J  exP(2^tan4,)  -  K  cot  * 


=  K  cos(4>  -  2n)  exp( 2tptarnf>) 


(57b) 


for  the  internal  friction  theory. 
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Indentor  No.  1 
0°  <  n°  <  45° 

<M  0° 


Indentor  No.  1 
B  >  45° 

0°  <  n  <  45° 

<P  =  0° 


Indentor  No.  2 
0°  s  n  s  45° 

*  =  0° 


Indentor  No.  2 


Indentor  No.  3 

0°  <  n  <  45° 
<P  *  0 


Indentor  No.  3 
n  ■  0° 

4>  =  0° 


Indentor  No.  4 


The  force-penetration  characteristic  can  be  obtained  by  integrating 
the  interface  stresses  along  the  tooth  profile  (the  formulas  for  different 
indentors  and  surface  conditions  are  listed  in  Appendix  D-Ref.  68).  The  work 
accompanying  the  penetrating  force  has  been  integrated  for  unit  depth 
using  the  equation 


FdH 


Numerical  Results 

The  numerical  results  for  four  indentors  were  obtained  for  the 
maximum  shear  theory  (Fig.  32  to  38),  and  for  the  internal  friction 
theory  (Fig.  40  to  Fig.  45). 

The  force-depth  curves  for  indentor  No.  2  (Fig.  32)  are  nearly  linear, 
and  the  magnitude  of  the  penetration  force  is  proportional  to  the  width 
d  of  the  bit  tooth.  As  d  approaches  zero,  the  curves  become  linear, 
which  is  the  condition  for  indentor  No.  1.  The  force  necessary  to 
penetrate  the  rock  increases  rapidly  for  both  increase  in  friction  and 
the  width  cf  the  nose  of  the  indentor.  The  penetration  forces  for  dif¬ 
ferent  interface  conditions  (varying  n)(Fig.  33)  increase  nonli nearly 
with  respect  to  n-  This  nonlinear  relationship  is  also  true  for 
different  indentors  and  wedge  angles (Fig.  35  and  Fig.  41). 

The  force-angle  and  work-angle  values  for  a  unit  depth  penetration 
(Fig.  34)  increase  rapidly  with  an  increase  in  4>  >  45°.  Static  energy 
was  calculated  by  integrating  the  current  force  at  depth  H  from  zero  to 
a  unit  depth.  For  the  force-penetration  characteristic  for  the  cylindrical 
indentor  (Fig.  35)  as  H  increases  the  force  approaches  a  constant  value 
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Force  (F/Kb) 
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Force  (F/Kb) 


and  the  assumed  slip-line  field  is  doubtful  above  a  certain  depth  of 
penetration.  The  dashed  line  sections  show  the  uncertain  region  of  the 
slip-1 ine  field. 

The  force-depth  relationships  for  round  nosed  indentors  (Fig.  36  &  38) 
No.  3  and  No.  4  for  different  values  of  (3  and  n  increase  rapidly  with 
the  initial  penetration,  but  level  off  rapidly  for  small  values  of  6. 

The  force-angle  and  work-angle  for  a  unit  depth  of  penetration  (Fig.  38)  are 
less  than  for  wedge  shaped  bits  for  large  values  of  8.  The  variation  of 
the  penetration  force  from  a  wedge  indentor  (r=0)  to  a  round-point  dull 
indentor  shows  that  an  increasing  amount  of  energy  is  consumed  as  a 
wedge  indentor  wears  to  a  dull  rounded  point. 

Figures  40  to  45  were  plotted  for  the  internal  friction  theory  using 
the  same  coordinates  as  for  the  maximum  shear  theory,  which  shows  that  the 
penetration  force  increases  exponentially  with  respect  to  ip  as  well  as  to 
41- 

Since  ip  is  a  function  of  n»  the  penetration  force  also  increases  with 
increasing  interface  frictiun  (Figs.  41  and  43).  The  rapid  increase  of 
force  with  the  internal  friction  ($)(Fig.  46)  is  critical  above  values 
of  30°  to  40°  depending  upon  the  value  of  6. 

Energy  Calculations 

From  the  stress  field  study,  for  a  given  eiergy  per  blow  of  a  pneumatic 
hammer  the  depth  of  penetration  can  be  predicted  from  the  energy-depth 
curve  (Fig.  47).  To  illustrate  effects  of  parameters  the  following  values 
are  selected: 
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0  =  40 
b  =  2" 


n  =  0° 

3  =  0° 

r  =  1/32  in.  and  r  =  0.  in. 

The  above  values  are  applicable  to  the  bits  used  in  the  experimental  part 

of  this  work.  If  the  compressive  strength  a  is  known  for  the  relationship 
a  costp  0 

K  =  Tj-p— -  -y  ,  the  K  values  for  the  following  rocks  are  obtained: 


Rocks 

K(psi) 

Charcoal  Grey  Granite 

7,759 

Kitledge  Pink  Granite 

5,662 

Dresser  Basalt 

10,112 

Black  Granite  (Buena  Gabbro) 

7,549 

Missouri  Red  Granite 

5,380 

Sioux  Quartzite 

11 ,347 
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Force  (F/Kb) 


70°  60° 


50' 

S1 


Fig.  38.  Force-Angle  Diagram 
and  <j>  -  0°. 


Indentor  No.  1 

<J> 


0 


T  '  7 


Indentor  Nd.  2 


6  > 


TT  (p 

V  '  7 


a  > 


I 


Indentor  No.  4 


Force  (F/Kb) 


3(-Force) 

i _ I _ i _ i _ i _ i _ i _ i _ l _ l 

90°  80°  70°  60°  50°  40°  30°  20°  10°  0° 

3(-Work) 

Fig.  42.  Force-Angle  Diagram  for  Indentors  No.  1  and  No.  2  for  H  =  1. 
n  =  0®,  and  $  =  40° 
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The  jackhammer  output  energies  per  blow  (65%  efficiency)  are: 


Throttle  Setting 
1 
2 

3 

4 


Energy  (ft-lb/blow) 
8.16 
16.35 
18.10 
19.15 


Now  let  us  assume  that  the  dimensionless  length  equals  1/32  inch. 
For  a  given  K,  b  and  W  (energy/blow) ,  jjjg-  can  be  calculated  and  from 
Fig.  47  the  penetration  depth  is  obtained.  For  example,  for  Grey 
Granite,  Throttle  # 3 

K  =  7,759  psi,  W  =  18.10  ft-lb/blow 
b  -  2  in.,  r  =  1/32  in.  and  0  in. 


Then, 


W  .  18.1  x  12  x  32  _  „ 

Kb  ’  7759  ..  „  ..  "  ~ 

32~x~  3?  x  2  x  32 


From  Fig.  47,  for  =  14.33, 


H  -  .21  x  1/32  =  .00656  in.,  for  r  =  1/32  in. 
H  =  .36  x  1/32  =  .01124  in.,  for  r  =  0  in. 


The  other  penetration  depths  for  different  rocks  and  conditions  are  listed 
in  Table  4. 
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Table  4.  Single  blow  penetration  depth 


Depth  in. 

(r  =  1/32 

in.) 

Depth  in. 

O 

II 

in.) 

ROCKS 

Thro. 
No.  1 

Thro. 
No.  2 

Thro. 

No.  3 

Thro. 

No.  4 

Thro. 
No.  1 

Thro. 

No.  2 

Thro. 

No.  3 

Thro . 
No.  4 

Grey  Granite 

.0038 

.0061 

.0066 

.0068 

.0075 

.0106 

.0112 

.0116 

Kit! edge  Pink 
Granite 

.0045 

.0077 

.0083 

.0086 

.0086 

.0125 

.0131 

.0136 

Basalt 

.0030 

.0051 

.0054 

.0056 

.0064 

.0094 

.0098 

.0101 

Buena  Black 
Granite 

.0039 

,0063 

.0067 

.0070 

.0076 

.0108 

.0114 

.0117 

Missouri  Red 
Granite 

.0047 

.0079 

.0085 

.0088 

.0090 

.0128 

.0135 

.0139 

Sioux  Quart¬ 
zite 

.0028 

.0045 

.0051 

.0052 

.0063 

.0086 

.0094 

.0096 

If  depths  of  penetration  are  known  then  the  required  energy  can  be 
calculated.  The  region  of  triaxial  compression  (false  nose)  extends 
beyond  the  depth  calculated  in  Table  4.  The  actual  depth  can  be  determined 
experimental ly. 

The  above  studies  of  a  single  bit  penetration  on  a  semi-infinite  rock 
give  preliminary  information  on  surface  chipping  of  rocks  for  different 
indentors.  For  evaluation  of  disintegration  phenomena  and  the  process  of 
kerf  ridge  removal,  further  studies  on  rock  properties,  interface  friction, 
dynamic  fracture,  indexing,  and  displacement  field  coupling  will  be 
necessary. 
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CHAPTER  V 


SPECIFIC  ENERGY 


Introduction 

Although  a  large  number  of  attempts  have  been  made  to  ascertain  the 
basic  relationship  between  the  particle  size  distribution  of  crushed 
and  ground  material  and  the  energy  required  to  accomplish  the  degradation, 
empirical  equations  have  been  established  which  apply  for  the  most  part 
only  to  the  particular  experiments  performed. 

Fundamentally,  the  amount  of  energy  consumed  in  fracture  and  crush¬ 
ing  is  proportional  to  the  new  surface  area  formed,  provided  the  original 
material  is  homogeneous  and  isotropic.  Likewise,  if  this  law  is  to  be  ap¬ 
plied  in  the  form  of  the  well-known  Rittinger  equation,  then  the  particles 
must  all  be  assumed  to  be  the  same  shape  and  free  from  cracks  resulting 
from  the  crushing  process. 

Crushing  Laws 

Charles  (Ref. 69  )  has  shown  that  both  the  Kick  and  Rittinger  laws  of 
crushing  are  special  cases  of  the  expression 

/•Xz  n 

E  =  f  -  C  dx/xn  (58) 

where  n  =  1  and  2  for  the  above  two  laws,  respectively,  arid  x  is  the 
particle  dimension. 

He  also  pointed  out  that  for  a  given  cumulative  distribution 

y  =  f(x)  (59) 

where  y  =  weight  percent  less  than  a  size  x,  the  energy  required  to 
reduce  a  weight,  dy,  from  size  xm  to  x  is 


dE,  = 


C  dxj/xi(dy) 


(60) 
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The  total  energy  is  given  by 


'■/  7 

•'ft  ~  ft 


100 


dEidy 


(61) 


or 


k  x 

-//  (-C  dxi/xi )f 1 (x)dx  (62) 


x„  xm 
o  m 

where  k  is  the  size  modulus,  taken  as  the  value  of  x  (largest  size)  at 
the  100  percent  intercept,  and  C  is  a  constant.  For  the  Schuhman  equation 

y  =  100(x/k)“  (63) 

the  energy  relationship  may  be  derived  from  the  above  for  large  xm  and 


small  xQ  as 


c  _  Ca  - 1 -n 

E  '  (n-1 )  (a~n+T)'  k 


(64) 


That  is,  for  constant  a,  C  and  n,  a  plot  of  E  versus  k  permits  the  evalua¬ 
tion  of  n  for  a  particular  process  and  for  crushing  equipment  design. 

On  the  basis  of  single  fracture  processes  and  the  occurrence  of  flaws 
according  to  the  Poisson  distribution  law  for  edge  flaws,  area  flaws  and 
volume  flaws  (Ref.  70)  Gilvarry  derived  an  expression  for  cumulative 
fractions  for  geometrically  similar  fragments  as: 


y  =  1 


■4M  $Y*(jY] 


(65) 


where  k,  j,  &  are  constants  with  units  of  length  denoting  the  spacing  of 
the  three  types  of  flaws,  respectively,  noted  above.  For  small  x  the 
(x/k)  term  becomes  dominant  and  corresponds  to  the  hypothesis  that  "the 
dominant  mechanism  producing  fragmentation  is  the  activation  of  edge 
flaws"  (Ref. 70).  The  activation  of  surface  and  volume  flaws,  then,  is 
important  only  when  a  small  number  of  large  fragments  is  produced. 
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It  should  be  noted  that  in  addition  to  the  assumption  that  all  particles 
are  geometrically  similar,  the  Poisson  distribution  function  does  not  al¬ 
low  for  flaw  depletion. 

In  a  recent  analysis  of  thermal  crushing  of  rock,  Thirumalai  and 
Cheung  (Ref.  71)  assumes  that  the  volume  flaws  dominate  the  fragmentation 
of  spherical  specimens  of  rock  subjected  to  sudden  intense  surface  heat¬ 
ing,  giving 

y-(f)‘  (66) 

for  x  <  <  l.  Thus,  equation  (66)  reduces  to 


A  study  of  the  empirical  laws  of  crushing  and  grinding,  and  experi¬ 
mental  results  obtained  for  various  types  of  crushing  show  that  the  amount 
of  energy  required  is  a  function  of  the  properties  of  the  material  being 
crushed,  the  method  of  crushing,  the  size  of  the  particles  produced,  hetero 
geneity,  distribution  of  flaws,  and  several  other  factors. 

Specific  Energy  and  Resistance  to  Fragmentation 

Some  basic  factors  in  the  investigation  of  thermal -mechanical  rock 
fragmentation  i)  the  energy  required  for  fracture,  (2)  energy  distri 

bution  in  fr  oups  (screen  fractions)  according  to  size,  and  (3)  a 

simple  test  wh  n  be  employed  to  predict  resistance  to  fragmentation. 

A  drop  crushing  test  utilized  by  Gross  and  Zimmerly  (Ref.  8)  has  been 
employed  by  several  investigators  to  approximate  the  energy  required  to 
crush  brittle  materials.  The  apparatus  consists  essentially  of  a  steel 
mortar  (anvil)  with  a  plunger  or  piston.  Cylindrical  specimens  of  rock 
are  placed  in  the  mortar  under  the  piston  and  a  weight  is  dropped  on  it 
(fig.  48).  The  energy  absorbed  by  the  system  is  subtracted  from  the  total 
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Fig.  48.  DROP  TEST  MACHINE 


drop  weight  energy  and  the  remainder  is  assumed  to  be  used  in  crushing. 
The  crushed  material  is  screened  (Appendix  E)  and  the  Rittinger  constant 
is  evaluated  by  the  equation 


E 

r 


%  i  _  1 00 " 

$  .d .  3 

II  o  J 


(68) 


where 


Er  =  crushing  energy  per  unit  volume 

Kr  =  Rittinger  constant 

7o.  -  percent  of  given  particle  size 

d.  =  average  particle  diameter  of  i  th  screen  fraction 

dQ  =  diameter  of  unbroken  specimen 

S.  =  shape  factor  for  i  th  screen  fraction  (S.  =  1.0  for  drop  tests) 

Equation  (68)  is  derived  from  Rittinger' s  statement  that  the  total  energy 
required  for  crushing  is  proportional  to  new  surface  area  formed.  It  as¬ 
sumes  that  all  of  the  particles  formed  are  of  the  same  shape.  To  adjust 
for  irregular,  flat  or  elongated  particles  an  average  or  effective  particle 
shape  factor  (Appendix  E)  is  introduced,  its  value  being  less  than  one. 

In  practice  its  value  is  determined  from  observation  of  particle  shapes 
in  each  screen  size  fraction  (Appendix  E). 

The  theories  that  relate  energy  to  surface  area  assume  (Ref.  72)  that 
the  particles  are  spherical  or  cubic  in  shape.  Since  some  of  the  particles 
collected  from  the  thermal -mechanical  fragmentation  tests  were  flat  chips, 
it  was  necessary  to  apply  a  particle  shape  correction  factor. 

Assuming  that  all  particles  were  cubes  having  side  d  gives  a  surface 
area  of  6d2.  If  the  cube  is  split  in  half,  the  surface  becomes  6d2  +  2d?. 
For  each  additional  split  2d '  surface  area  is  created.  Thus,  cutting  the 
cube  by  a  number  of  parallel  planes  results  in  an  increased  surface  area 
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equal  to  twice  the  number  of  resulting  parallelopipeds  less  one  times 
the  square  of  the  diameter.  Expressing  the  average  thickness  of  the 
particles  as  a  percentage  of  the  observed  diameter  will  give  the  new  sur 
face  area  as: 

-  l^d2 

where 

d  =  observed  diameter 

6  =  ■— r  —  ,  thickness  of  particles  in  percent  of  observed  diameter 
Letting  c  =  ^  the  equation  becomes: 

surface  =  d2 ( 4  +  2c) 


surface  =  6d2  +  2 


For  the  divided  cube,  the  specific  surface  (surface  per  unit  weight) 
becomes: 

S  _  Sn  _  (4  +  2c)d2  _  (4  +  2c) 

w  weight  ~  d3p  dp 


where 


Sw  =  specific  surface 

S  =  surface  of  divided  cube 
n 

p  =  density 


From  statistics,  the  corrected  diameter  (harmonic  mean  diameter)  is 
expressed  by 

h  -  6 

a  =  — 

cor  p5 

w 

where 

d  =  corrected  diameter, 
cor 

Substituting  for  S  gives: 

w  J 


d 

cor 


_  JL  -  6 

"  uS"  ' 


6  ^  _  3t 

4  +  2c  2t  +  d 


(69) 
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For  a  given  sample,  average  thicknesses  of  particles  retained  on  dif¬ 
ferent  screens  were  obtained  and  were  plotted  against  their  corresponding 
average  screen  diameters  on  a  log-log  scale.  The  points  were  observed  to 
lie  approximately  on  a  straight  line.  Since  it  was  not  possible  to  measure 
the  thicknesses  of  small -sized  particles,  a  logarithmic  least-square  straight 
line  (Ref. 73)  was  passed  through  the  points  for  which  thicknesses  were  mea¬ 
surable. 

The  equation  of  the  least-square  straight  line  is  given  by 

log10  t  =  A  log10  d  +  B.  (70) 

The  constants  A  and  B  are  obtained  from  solution  of  the  normal  equations 

E  1  og i o  t  =  A^logio  d  +  nB 

n  n 


£(1o9io  t-logio  d)  -  A£(log10  d)2  +  B^logio  d 
n  n  n 


and  are  given  by 

^El0 9io  dj^£lo9io  t^-  n  E^ 


^E  logic  dj2  -  n^logl0 


1  og i o  dj  ^logio  t^ 


B  = 


Eiog  1  o  t  -  aE  log  i  o  d 

_n _ n _ 

n 


In  the  above  equations  n  is  the  total  number  of  points. 

The  point  of  intersection  of  the  least-square  straight  line  with  the 
average  screen  diameter  equal  to  the  average  particle  size  line,  t  =  d,  is 
given  by 

d  =  t  =  dA-10B  (71) 
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For  screen  sizes  larger  than  that  given  by  equation  (71)  the  thick¬ 
nesses  are  modified  using  the  least-square  straight  line,  equation  (70). 

The  modified  thicknesses  are  used  in  equation  (69)  to  correct  the  dia¬ 
meters.  For  screen  sizes  smaller  than  that  given  by  equation  (70),  the 
particles  are  assumed  to  be  cubic  and  no  correction  for  the  diameter  is 
requ i red . 

Once  the  Rittinger  constant  Kr  is  determined  for  a  given  rock  the 
value  of  the  breakage  energy  of  the  fragments  from  chipping  and  thermal 
spallation  may  be  approximated  from  screen  analysis  as  well  as  the  energy 
distribution  among  the  particle  size  fractions. 

Usually  the  factor  100/dQ  is  small  in  comparison  to  the  other  terms 
and  may  be  neglected,  in  which  case  equation  (68)  becomes 

Er  =  KrSsT  (72) 

Assuming  that  «r  is  known  and  a  given  specific  energy  is  desired,  an 
equivalent  particle  diameter  of  material  may  be  calculated.  For  example, 

O 

for  «r  =  3.0  x  10  ,  Er  =  10  joules/cc  and  S  =  1.0,  the  equivalent  particle 

-4 

diameter  is  =  3  x  10  cm.  This  calculation  is  in  a  sense  an  estima¬ 
tion  of  the  effective  average  particle  surface  area  which  should  be  sought 
in  a  given  process  if  the  specific  energy  of  10  joules/cc  is  desired. 

The  shape  factors  were  calculated  for  fragmented  rocks  as  a  function 
of  particle  size,  as  previously  described,  larger  corrections  being  required 
for  large  particles,  the  factors  being  employed  for  the  chipping  tests  only. 
Measurements  of  particles  from  the  drop  crushing  tests  indicated  that  they 
were  approximately  equidimensional  and  no  correction  factor  was  applied. 

The  probability  that  incomplete  fractures  are  induced  in  otherwise  intact 
particles  is  hypothetically  a  function  of  particle  size.  The  existence 
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of  partial  fractures  represents  surface  energy,  and  could  be  entered 
into  the  Rittinger  equation  by  decreasing  the  value  of  the  shape  factor. 
However,  for  particle  size  distributions  having  even  small  amounts  of 
fine  material,  the  latter  dominates  the  calculations  and  the  shape  fac¬ 
tors  have  little  effect  on  the  value  of  l<y. 

In  view  of  the  above  it  is  logical  to  assume  that  the  effective  Rit¬ 
tinger  values  are  useful  to  approximate  the  actual  energy  of  breakage  for 
particle  size  distributions  that  are  reasonably  similar.  The  limits  of 
parameters  which  define  such  similarity  have  not  been  investigated  and 
are  to  be  the  subject  for  further  research. 

The  values  in  Table  5  and  the  figures  in  Appendix  E  show  the  Rittinger 
constant  actually  varies  with  a  number  of  test  parameters.  Its  value  de¬ 
pends  on  drop  height,  drop  weight,  specimen  geometry,  the  number  of  drops 
and  other  factors.  This  variation  can  be  minimized  by  standardization 
of  the  test  procedure. 

Drop  crushing  tests  (Table  5)  were  conducted  with  the  specimens 
placed  on  edge  as  in  the  Brazil  tension  test.  Just  as  Protodyakonov  used 
an  impact  test  on  multiple  fragments  as  a  measure  of  hardness  it  is  ex¬ 
pected  that  with  some  modification  of  the  drop  crushing  test  it  will  be 
possible  to  determine  a  controlling  factor  in  the  resistance  of  rocks 
to  thermal -mechanical  fragmentation  by  a  standardized  crushing  test. 

The  values  of  the  Rittinger  constant  for  thermally  treated  and  untreated 
specimens  may  provide  a  measure  of  the  relative  resistance  to  thermal - 
mechanical  fragmentation. 

The  column  titled  "Specific  Surface  Breakage  Energy"  in  Table  9 
shows  how  the  Rittinger  constant  as  determined  from  drop  tests  can  be 
used  in  a  fragmentation  test  to  determine  the  response  of  the  rock  to 
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Table  5.  Rock  Fragmentation  (Crushing)  Constants 


Energy 


Spec¬ 

Drop 

into 

Test 

imen 

No. 

Height 

Rock 

K  2 

K,3 

Kh4 

No. 

Specimen 

Type1 

Drops 

(ft) 

(joules) 

r 

K 

0 

1 

Missouri  Red 

Granite 

1 

1 

5 

626.47 

.00566 

.33858 

.08508 

2 

tl 

1 

1 

10 

544.71 

.00303 

.27299 

.05880 

3 

II 

1 

2 

1 

401.10 

.00485 

.24740 

.07030 

4 

II 

1 

2 

3 

592.37 

.00317 

.26774 

.05909 

5 

II 

1 

2 

5 

521.52 

.00299 

.23146 

.05171 

6 

II 

2 

1 

2 

417.38 

.00491 

.29280 

.07477 

7 

II 

2 

1 

4 

382.39 

.00410 

.25262 

.061 74 

8 

II 

2 

1 

5 

416.97 

.00379 

.27633 

. 06400 

9 

i: 

2 

1 

10 

257.64 

.00188 

.15207 

.03305 

10 

Charcoal  Grey 

Granite 

1 

1 

5 

509.86 

.00449 

.28331 

.07079 

11 

II 

1 

1 

10 

660.85 

.00412 

.33665 

.07677 

12 

II 

1 

2 

1 

463.75 

.00613 

.29134 

.08476 

13 

II 

2 

1 

2 

440.56 

.00558 

.35271 

.09093 

14 

II 

2 

1 

4 

533.31 

.00454 

.36135 

.08248 

15 

li 

2 

1 

5 

449.38 

.00444 

.31293 

.07470 

16 

II 

2 

1 

10 

440.56 

.00281 

.23838 

.051 92 

17 

Kitledge  Pink 

Granite 

1 

1 

5 

672.58 

.00444 

.35375 

.07906 

18 

ll 

1 

1 

10 

568.16 

.00361 

.27944 

.06302 

19 

II 

1 

2 

1 

463.75 

.00429 

.26031 

.06626 

20 

II 

2 

1 

2 

417.38 

.00437 

.27233 

.06775 

21 

II 

2 

1 

4 

371.00 

.00320 

.22442 

.05231 

22 

II 

2 

1 

5 

223.81 

.00205 

.12891 

.03063 

23 

It 

2 

1 

10 

173.91 

.00143 

.09396 

.021 56 

24 

Black  Granite 

(Buena  Gabbro) 

1 

1 

5 

695.63 

.00791 

.43169 

.11540 

25 

II 

1 

1 

10 

463.75 

.00391 

.25957 

.06383 

26 

ll 

1 

2 

1 

463.75 

.00682 

.34256 

.09952 

27 

II 

2 

1 

2 

417.38 

.00690 

.39482 

.10851 

28 

ll 

2 

1 

4 

510.13 

.00645 

.38422 

.09749 

29 

II 

2 

1 

5 

579.69 

.00486 

.41011 

.09170 

30 

II 

2 

1 

10 

243.40 

.00221 

.16490 

.03770 

31 

Sioux  Quartzite  1 

1 

5 

465.11 

.00484 

.29000 

.07630 

32 

II 

1 

1 

10 

676.37 

.00396 

.36273 

.07906 

33 

ll 

1 

2 

1 

417.38 

.00690 

.30176 

.09211 

34 

II 

2 

1 

2 

463.75 

.00758 

.42778 

.11468 

35 

ll 

2 

1 

4 

556.50 

.00608 

.41567 

.10022 

36 

ll 

2 

1 

5 

672.30 

.00620 

.44858 

.10698 

37 

II 

2 

1 

10 

81 1 . 57 

.00751 

.49679 

.11805 
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Table  5.  Rock  Fragmentation  (Crushing)  Constants  -  Continued 


Test 

Specimen 

Spec¬ 

imen 

No. 

Drop 

Height 

Energy 

into 

Rock 

K  2 

K  3 

K  4 

No. 

Type1 

Drops 

(ft) 

(joules) 

r 

> 

Kr 

38 

39 

Dresser  Basalt 

ll 

1 

1 

1 

o 

1 

1 

5 

846.14 

.01811 

.61479 

.20823 

40 

ll 

10 

927.50 

.01186 

.52120 

.15473 

41 

ll 

2 

1 

1 

1 

1 

1 

463.75 

.01891 

.39242 

. 1 5849 

42 

<1 

c 

2 

417.38 

.01 564 

.46303 

.17136 

i  b 

43 

II 

Z 

o 

4 

579.69 

.01446 

.55401 

.18307 

44 

ll 

c 

o 

5 

626.06 

.01594 

.58814 

.19185 

C 

10 

347.81 

.  00858 

.29238 

.09251 

1 

One  specimen  (1), 
Two  specimens  (2) 

2 


2.75"  long  by  1"  diameter  placed  on  edge 
,  each  1"  long  by  1"  diameter  placed  on  edge 


Rittinger  constant 


3«ick  constant 


4 

Bond  constant 


118 


variations  in  different  test  parameters.  Particle  size  distributions 
from  fragmentation  tests  may  be  used  with  the  Rittinger  constant  to 
indicate  the  specific  energy  differences  due  to  bit  size  and  style, 
offset,  kerfing,  simultaneous  versus  independent  operation  of  heat 
source  and  mechanical  device,  traverse  speed,  etc. 

Total  Crushing  Energy 

Many  of  the  empirical  crushing  laws  utilize  the  total  process  input 
energy  to  determine  that  required  for  particle  size  reduction.  However, 
for  thermal,  mechanical,  explosive,  or  other  fragmentation  processes 
there  is  an  inherent  amount  of  waste  energy  which  is  not  utilized  directly 
for  fracturing. 

For  example,  in  the  process  of  blasting  with  a  completely  confined 
explosive  as  much  as  75  percent  of  the  total  available  energy  from  the 
gaseous  products  of  detonation  may  be  absorbed  by  the  rock.  However,  de¬ 
pending  upon  the  geometry  of  the  rock,  i.e.,  number  of  free  faces,  as 
little  as  about  1/6  of  the  energy  may  be  utilized  in  fragmentation,  the 
remainder  being  dissipated  into  the  surrounding  rock  as  stress  waves, 
and  in  kinetic  energy  of  flying  fragments,  etc.  For  mechanical  breakage 
by  a  pneumatic  drill,  only  a  portion  of  the  energy  of  impact  causes  pul¬ 
verization,  fragmentation  and  fracture,  the  remainder  being  dissipated 
into  the  rock,  in  friction  losses,  etc. 

For  thermal  fragmentation  by  sources  applied  to  a  plane  surface 
much  heat  energy  may  be  lost  to  the  surroundings,  and  of  the  heat  trans¬ 
mitted  to  the  rock  only  a  fraction  of  it  is  absorbed  in  fracture  processes, 
the  remainder  showing  up  as  residual  heat  in  the  rock. 

The  calculated  values  of  the  effective  Rittinger  constant  determined 
by  drop  crushing  tests  are  assumed  to  be  representative  of  the  fracture 
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(comminution)  energy  only,  and  hence,  can  be  utilized  as  a  basis  for 
determining  the  approximate  efficiency  of  other  fragmentation  processes 
resulting  in  similar  particle  size  distributions  for  tne  same  material. 

They  also  furnish  approximate  energy  values  for  defining  an  engineering 
index  property  for  resistance  to  fragmentation. 

For  calculation  of  energy  requirements  and  work  indices  for  industrial 
crushing,  the  parameters  are  often  defined  for  a  certain  percentage  of 
crushed  material  less  than  a  certain  size,  and  others  consider  only  particle 
size  distributions  down  to  a  certain  cutoff  size,  neglecting  very  fine 
material.  Hence,  comparisons  of  energy  consumption  recorded  in  the  litera¬ 
ture  for  the  several  fragmentation  processes  in  use  or  under  investigation, 
should  take  into  account  possible  differences  in  the  bases  for  calculations. 

For  many  brittle  materials  the  surface  energy  may  be  taken  as  a  mea¬ 
sure  of  useful  energy  expended,  and  the  Rittinger  law  (Fig.  49)  as  a  good 
first  approximation  for  specific  energy  values.  Factors  discussed  above 
would  cause  apparent  deviations  from  the  Rittinger  law  depending  upon  how 
the  energy  input  is  evaluated.  For  all  empirical  laws,  however,  it  is  ob¬ 
vious  that  energy  consumption  can  be  drastically  reduced  by  lessening  the 
amount  of  fine  material  produced  in  the  fragmentation  process. 
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CHAPTER  VI 


THERMAL-MECHANICAL  FRAGMENTATION  TESTS 

Introduction 

In  order  to  reduce  scaling  (similitude)  problems  to  a  minimum,  to 
utilize  full-scale  equipment,  and  to  provide  large  enough  specimens  to 
furnish  a  reasonable  representation  of  large  rock  masses,  all  thermal 
fragmentation  tests  were  carried  out  on  blocks  of  rock  about  30  inches 
on  each  side.  It  was  desirable  to  vary  parameters  such  as  speed  of 
traverse,  magnitude  of  heat  flux,  intensity  and  frequency  of  bit  impact, 
as  well  as  characteristics  of  the  bits  employed.  Inasmuch  as  an  evalua¬ 
tion  of  specific  energy  of  fragmentation  is  of  vital  interest,  the  ex¬ 
perimental  system  was  designed  so  that  cuttings  from  drilling  and  frag¬ 
ments  from  spallation  could  be  collected  for  screen  analysis.  Although 
some  finer  particles  were  losr.,  it  is  believed  that  the  losses  were  not 
significant  except  on  the  preliminary  series  of  tests  on  Missouri  Red 
Granite,  which  were  used  to  determine  the  operating  characteristics  of 
the  equipment. 

Thermomechanical  Test  Equipment* 

Mounting  Frame.  The  rock  test  specimens  for  thermomechanical  frag¬ 
mentation  as  indicated  above  were  30  in.  (approximately)  cubes  weighing 
about  one  and  one-half  tons.  The  test  apparatus  was  designed  so  that  it 
would  mechanically  move  the  test  blocks  longitudinally  across  the  points 
of  action  of  the  pneumatic  hammer  and  flame  jet.  A  steel  framework  was 

★ 

Note:  Fquipment  designed,  built  and  tested  under  THEMIS  Contract  No. 
DACA-45-69-C-0087. 
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Fiy.  50.  Schematic  of  Experimental  Apparatus  for  Thermo- 
mechanical  Fragmentation 
123 


constructed  for  mounting  the  latter  equipment  and  the  system  (Figs.  50  &  51) 
provided  transverse  as  well  as  vertical  movement  of  the  hammer  and  jet  re¬ 
lative  to  the  rock.  The  cutting  face  of  the  rock  was  placed  in  a  vertical 
position  to  prevent  secondary  breakage  of  particles  by  the  jackhammer  by 
allowing  the  particles  to  fall  from  the  vertical  face  of  the  rock.  Also, 
this  permitted  the  thrust  of  the  jackhammer  against  the  rock  to  be  con¬ 
trolled  by  means  of  an  air  cylinder. 

Flame  Jet.  An  FSJ-6  stone-shaping  torch  was  the  primary  source  of 
heat  for  thermal  weakening/spalling.  The  torch  was  operated  on  oxygen 
and  No.  1  fuel  oil  at  pressures  and  flow  rates  suggested  by  the  manufac¬ 
turer  (Table  6).  Stoichiometric  calculations  give  the  maximum  combustion 
energy  as  1.78  x  105  j/sec  for  the  recommended  operating  conditions.  Com¬ 
plete  combustion  of  hydrocarbons  produces  a  blue  flame,  but  the  flame  jet 
has  an  orange-yellow  flame  indicating  the  presence  of  unburned  carbon 
particles.  Losses  from  incomplete  combustion  and  other  waste  heat  such 
as  the  energy  in  the  cooling  water  reduce  the  maximum  available  energy. 

The  actual  energy  available  to  the  rock  is  reduced  because  of  the  rela¬ 
tively  small  portion  of  the  gas  in  contact  with  the  rock  and  the  short 
residence  time,  resulting  in  large  residual  heat  and  kinetic  energy  losses 
in  the  high  velocity  exhaust  gases.  Losses  due  to  vaporization  of  the 
cooling  water  were  avoided  by  modifying  the  torch  so  that  the  water  was 
not  directed  against  the  rock  surface.  The  standoff  distance  from  the 
torch  to  the  rock  surface  and  the  rate  of  traverse  of  the  rock  are  two 
primary  factors  which  affect  the  heat  coupling  efficiency  and  the  conse¬ 
quent  energy  input  into  the  rock.  Estimates  indicate  that  the  maximum  heat 
energy  available  to  the  rock  varies  from  6  to  8  percent  of  the  total  com¬ 
bustion  energy  (Ref.  19). 
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The  temperature  of  the  flame  is  approximately  2400  deg  C  and  the 
velocity  of  the  exhaust  gas  varies  from  4,000  to  10,000  ft/sec  (Ref.  19). 
The  combustion  chamber  is  one-inch  inside  diameter  and  the  exhaust 
nozzle  is  3/16  in.  DeLaval  type. 

Table  6.  Recommended  Operating  Conditions  for  FSJ-6  Flame  Jet 


Operation  Operation 


Consumable 

Flow  Rate 

Pressure 

(lb/in.2) 

Flow  Rate 

Pressure 

(lb/in.2) 

Oxygen 

100  (ft3/hr) 

150 

500  (ft3/hr) 

135 

No.  1  Fuel  Oil 

15  (lb/hr) 

120 

22  (Ib/hr) 

120 

Water 

85  (gal/hr) 

65 

85  (gal/hr) 

65 

Pneumatic  Hammer.  The  mechanical  removal  of  rock  was  accomplished  by 
a  pneumatic  hammer.  The  work  output  of  a  pneumatic  hammer  is  given  by 
Pfleider  and  Lacabanne  (Ref.  74)  as: 

P  CP3/2A3/2S,/z 

ti  = - 7 - 

wi/2 

where 

Ei  =  pneumatic  hammer  work  output  (ft-lb/min) 

C  -  constant 

P  =  mean  effective  air  pressure  (lb/in.2) 

A  =  area  of  piston  face  (in.2) 

S  =  length  of  working  stroke  (in.) 

W  =  weight  of  piston  and  rifle  nut  (lb) 
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The  above  equation  is  derived  as  follows: 


Factor 

Equation 

Dimensions 

Force  on  working  face  of  piston 

F  =  PA 

(lb) 

Work  on  power  stroke 

_  PAS 
e  -  Tr 

(ft-lb) 

Acceleration  of  piston  on 

a  = 

a  w 

(ft/sec2) 

power  stroke 

Time  of  power  stroke 

from  S  =  1/2  at2 

t  _  WS  % 
z  ~  6PAg 

(sec) 

Time  of  piston  round  trip 

T  =  (HK)  5^ 

(sec/blow) 

Blow  per  minute 

60_  6PAg_ 

n  1+K  WS 

H 

(blow/min) 

Total  work  output  per  minute 
(work  per  blow  x  blows  per  min.) 

/zA: 

1/2  $1/2 

bl  T2TT+K) 

~W 

(ft-lb/min) 


The  constant  "c"  is  determined  from  the  value  of  K  (fraction  of  time  of 
piston  power  stroke  required  for  return  stroke),  since: 


C 


(60)  6*5  qh 

■“i?(t+k  r 
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The  mean  effective  pressure  is  a  function  of  thrust  and  throttle 

setting.  The  cylinder  jacket  of  the  pneumatic  hammer  was  drilled  and 

tapped  and  a  pressure  transducer  was  installed  to  monitor  the  pressure. 

The  average  output  of  the  transducer  was  taken  as  the  mean  effective 

pressure.  The  area  of  the  piston  for  the  pneumatic  hammer  used  in  this 

2 

research  was  5.363  in.  .  The  inlet  operating  pressure  was  maintained  at 

2 

approximately  90  lb/in.  in.  in  order  to  permit  the  use  of  the  manufac¬ 
turer's  value  of  2-9/16  in.  for  the  working  stroke.  The  weight  of  the 
piston  and  rifle  nut  is  4.31  lb. 

Thrust  Cylinder.  Initially  it  was  planned  to  position  the  jacKhammer 
relative  to  the  rock  and  to  clamp  it  during  a  test.  This  proved  unsatis¬ 
factory  because  the  jackhammer  could  not  move  longitudinally  with  the 
contour  of  the  rock  face.  A  screw  carriage  was  then  employed  on  which 
to  mount  the  jackhammer  which  allowed  the  operator  to  move  the  jackhammer 
during  the  test  to  prevent  stalling,  but  it  also  permitted  the  thrust  to 
vary.  Finally  an  air  thrust  cylinder  was  installed  which  permitted  trans¬ 
verse  movement  of  the  jackhammer  while  a  constant  thrust  was  maintained 

against  the  rock.  Calibration  of  the  air  cylinder  indicated  that  a  pres- 

2 

sure  of  13.2  lb/in.  was  necessary  to  produce  a  desired  thrust  of  80  lb. 

Rock  Table.  The  movement  of  the  rock  relative  to  the  flame  jet  and 
jackhammer  was  accomplished  by  placing  the  test  block  on  a  movable  table 
mounted  on  wheels  on  an  inverted  angle  iron  track.  The  speed  of  the  table 
was  calibrated  and  operating  speeds  selected  between  0.51  in. /sec  and 
4.5  in. /sec.  At  speeds  less  than  0.51  in. /sec,  the  jackhammer  began 
drilling  a  hole  rather  than  traversing  the  face. 
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Table  7.  Energy  vs  Table  Speed  for  Table  8 


Output  Energy3 
(joules) 

Table  Speed  0.50 
in/ sec 

Table  Speed  0.95 
in/ sec 

Table  Speed  1 .40 
in/sec 

Efficiency 

Percent 

Jackhammer1 

75,800 

40,350 

27,580 

65 

Flame  Jet2 

10.51  x  106 

5.59  x  106 

3.82  x  106 

2-7 

1  Based  on  77,350  joules/min 

2  Based  on  10.72  x  106  joules/min 

3  Thirty  inches  of  travel 
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Table  8.  Calibration  Data  for  J50A  Jackhammer* 


s-  c 

QJ  *r- 


C  E 

LU 

</) 

P  qj 
3  »— 
Cl  3 
P  O 

OnI2 


O  O  O 
CM  O  CM 
CO  l£>  (V) 

#i  M  • 

cn  o  cn 

r-CVJr- 


o  o  o 

CO  O  CO 
f—  «d“  i — 

r— *  r— 

LT)  Vt*  LO 


O  O  O 
CM  CO  CO 
VO  ID 

*  vs  #s 

cn  m  in 

VDkOVO 


o  o  o 

CM  CTi  CM 
LO  CO  ID 

in  r--  cn 
co  r-s  vo 


o  o  O 
i—  r-  cm 
*cf  in 

#*  *s  *» 

r—  CO  0> 

oo  co 


U  <D 
C  lO 
<D 

cn  ld  <n 

o  o  o 

o  o  o 

Q  O  O 

o  o  o 

d-  i 

CO  CO  CO 

co  vo  m 

Cn  (n  Cn 

o  o  cn 

i—  cn 

QJ  O 

• —  f —  r— 

CM  CM  CM 

CM  CM  CM 

co  ro  cm 

CO  CM  CM 

S-  f— 
Li-  O 

u 

> 

•r- 

P 

QJ  CM 

o 

L-  * 

0> 

3  C 

4- 

CO  •»- 

P 

LU 

0J  -O 
L-  »— 

C 

O- 

<T3 

<U 

s. 

mom  o  m  o  oo 

*—  CM  I—  CM  i—  CM  r^l£> 

r—  CM  CM  CM  CM  CM 


o 

o  o  o 

o  o  o 

• 

•  •  • 

•  •  • 

LO 

i— •  o>  r^. 

O  CO  N 

CM 

CO  CM  CM 

CO  CM  CM 

c 

o 

•r— 

p 

n3 

p 

10  #"*< 

o  O  O 

o  o  o 

o  o  o 

o  o  o 

o  o  o 

P 

i-  1— 

kO  N  CVJ 

ID  CM 

vo  r>.  cm 

vo  cm 

vo  r*.  cm 

•r— 

-C  ' 

CM  LO  CT» 

cm  lo  cn 

CM  LO  O 

CM  U*>  <Ti 

cm  lo  cn 

03 

o 

h~ 

<b 

f—  cn 
P  c 
P  T- 
O  P 
U  P 

-C  Q) 

h-  in 


cm  co  in 


CM 

c 

*r~ 

_Q 


o 

c 

cn 

•F— 

LO 

dJ 

S- 

3 

o 

</> 

1 

<1> 

S- 

CM 

O- 

<U 

a» 

o 

c 

J- 

■r- 

p 

P 

oo 

★ 


130 


i 

F 


Preliminary  Tests 


Initial  tests  were  conducted  with  the  thermal  fragmentation  equip¬ 
ment  to  determine  its  stability,  strength,  operating  characteristics,  and 
to  obtain  initial  test  results  to  establish  optimum  values  for  such  factors 
as  jackhammer  thrust,  table  speed,  and  standoff  for  the  flame  jet.  Some 
welds  were  strengthened  in  addition  to  the  modifications  noted  earlier; 
a  vacuum  system  to  collect  fine  particles  was  installed  after  this  first 
series  of  tests. 

The  variables  utilized  in  the  calibration  tests  were  table  speed, 
fragmentation  without  and  with  heat,  and  flame  jet  standoff  (Table  9). 
Eighteen  tests  were  run,  four  of  which  were  shallow  kerfing  experiments. 
Particulate  material  was  collected,  weighed  and  subjected  to  a  screen 
analysis.  The  specific  breakage  energy  was  then  calculated  using  a  value 
of  Kr  determined  from  drop  crushing  tests.  These  values  of  specific  energy 
are  low  because  some  of  the  very  fine  material  was  not  recovered  before 
the  subsequent  installation  of  the  vacuum  system. 

Some  of  the  usable  ranges  of  operation  of  the  segments  of  the 
apparatus  are  interdependent.  For  example,  if  the  thrust  of  the  jackhammer 
is  too  high  it  will  drill  a  hole  into  the  rock  and  will  jam  the  table. 

The  flame  jet  cannot  be  operated  simultaneously  with  the  jackhammer  because 
the  high  velocity  gases  disperse  the  fine  cuttings. 

Table  Speed.  For  the  initial  tests  (Table  9)  speeds  ranging  from  0.5 
to  1.4  in/sec  were  employed.  Higher  speeds  were  later  needed  and  a 
larger  pulley  was  installed  on  the  motor  drive  for  subsequent  experiments. 

In  the  first  group  of  tests  it  was  found  that  a  table  speed  of  0.95 
in/sec  gave  the  lowest  specific  energy  for  pneumatic  chipping  and  this 
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Table  9.  Thermomechanical  Fragmentation  -  Preliminary  Tests 

Total  Specific 

Surface  Surface  Specific 


Material 

Breakage 

Breakage 

Total** 

Energy 

Description  of  Test 

Removed 

Energy 

Energy 

Energy 

(joules/cc) 

(grams) 

(cc) 

(joules) 

(joules/cc) 

(joules) 

JH 

Total 

la 

TS 

0.50 

B1 

NK  w/o  heat 

413.6 

159 

4554 

28.63 

49,270 

309 

309 

b 

TS 

0.95 

B1 

NK  w/o  heat 

285.2 

no 

3220 

29.37 

26,227 

239 

239 

c 

TS 

1.40  B1 

NK  w/o  heat 

185.3 

71 

2286 

32.08 

17,927 

252 

252 

2a 

TS 

0.50 

B1 

NK  NO  10  in.* 

542.5 

209 

3522 

16.89 

364,570 

235 

1747 

b  TS 

0.95 

B1 

NK  NO  10  in. 

425.1 

164 

2944 

18.02 

193,927 

160 

1186 

c  TS 

1.40 

B1 

NK  NO  10  in. 

198.5 

76 

1686 

22.07 

132,527 

236 

1735 

3a 

TS 

0.95 

B1 

NK  w/o  heat 

285.2 

no 

3220 

29.37 

26,227 

239 

239 

b  TS 

0.95 

B1 

KA  w/o  heat 

491.8 

189 

5001 

26.44 

26,227 

138 

138 

c 

TS 

0.95 

B1 

KB  w/o  heat 

429.1 

165 

4213 

25.53 

26,227 

158 

158 

4a 

TS 

0.95 

B1 

NK  NO  10  in. 

425.1 

164 

2944 

18.02 

193,927 

160 

1186 

b  TS 

0.95 

B1 

KA  NO  10  in. 

869.1 

334 

3008 

8.99 

193,927 

79 

580 

c 

TS 

0.95 

B1 

KB  NO  10  in. 

673.0 

259 

3401 

13.15 

193,927 

101 

749 

5a 

TS 

0.95 

B1 

NK  w/o  heat 

285.2 

no 

3220 

29.37 

26,227 

239 

239 

b  TS 

0.95 

B2 

NK  w/o  heat 

237.1 

91 

2381 

26.12 

26,227 

283 

287 

c 

TS 

0.95 

B3 

NK  w/o  heat 

239.4 

92 

2781 

30.21 

26,227 

285 

284 

d 

TS 

0.95 

B4  NK  w/o  heat 

365.4 

140 

3646 

25.95 

26,227 

187 

186 

6a 

TS 

0.95 

B1 

NK  w/o  heat 

285.2 

no 

3220 

29.37 

26,227 

239 

239 

b  TS 

0.95 

B1 

NK  NO  14  in. 

344.1 

132 

2403 

18.16 

193,927 

199 

1465 

c 

TS 

0.95 

B1 

NK  NO  10  in. 

425.1 

164 

2944 

18.02 

193,927 

160 

1186 

d 

TS 

0.95 

B1 

NK  NO  6  in. 

796.4 

306 

4846 

15.83 

193,927 

86 

633 

7a 

TS 

0.95 

B1 

NK  NO  10  in. 

425.1 

164 

2944 

18.02 

193,927 

160 

1186 

b 

TS 

0.95 

B1 

NK  Off  1  10  in 

.  332.5 

128 

3368 

26.34 

193,927 

204 

1516 

c 

TS 

0.95 

B1 

NK  Off  2  10  in 

.  344.6 

132 

2689 

20.27 

193,927 

199 

1463 

8a 

TS 

0.95 

B1 

NK  NO  10  in.  S 

425.1 

164 

2944 

18.02 

193,927 

160 

1186 

b 

TS 

0.95 

B1 

NK  NO  10  in.  I 

1009.3 

388 

5307 

13.68 

193,927 

68 

499 

*  Number  in  inches  is  standoff  distance 

**  Based  on  jackhammer  efficiency  =  65% 
Flame  jet  efficiency  at  14  in.  =  2 fl¬ 
at  10  in.  =  3% 

IS  _  =  Table  speed  in. /sec 
Bl~  ~  -  1-1/2  in.  cross  bit 

B2  =  2-1/4  in.  cross  bit 

B3  =  1-3/4  in.  button  bit 

B4  =  2  in.  button  bit 


NK  =  No  kerf 

KA  =  Kerf  A  (1.2  in.  wide  x  0.3  in.  deep 

KB  =  Kerf  B  (1.6  in.  wide  x  0.2  in.  deep; 

NO  =  No  offset 

Off  1  =1  in.  offset 

Off  2  =  2  in.  offset 

S  =  Simultaneous 

I  =  Independent 
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speed  was  employed  in  the  remainder  of  the  preliminary  test  series.  It 
was  not  necessarily  the  optimum  speed  for  thermal  weakening  or  fragmentation. 

Flame  Jet  Standoff.  Three  standoff  distances  were  employed,  i.e., 

14,  10,  and  6  inches,  with  the  heat  transfer  efficiency  being  approximately 
2,  3,  and  7  percent,  respectively,  at  those  distances.  The  optimum  speed 
for  chipping  after  flame  jet  heating  at  10  in.  (test  group  3)  standoff 
was  0.95  in/sec.  This  also  gave  the  lowest  total  specific  energy  (1186  j/cc) 
and  the  lowest  chipping  specific  energy  (160  j/cc).  The  best  standoff 
(test  group  6)  was  6  in.  requiring  a  total  specific  energy  of  633  j/cc  and 
a  chipping  energy  of  86  j/cc.  Thus,  the  chipping  energy  is  reduced  by  a 
factor  of  3+  for  flame  jet  heat  applied  at  10  in.  or  6  in.  standoff. 

It  should  be  noted  that  the  flame  jet  was  used  for  heating  without  a 
cooling  water  spray. 

Shallow  Kerfing.  Shallow  kerfing  without  heat  at  a  traverse  speed 
of  0.95  in/sec  showed  the  lowest  specific  energy  for  narrow  deeper  kerfs 
(138  j/cc).  With  heat  the  total  energy  was  580  j/cc  while  the  chipping 
energy  was  79  j/cc. 

Bit  Characteristics.  Only  two  sizes  each  of  two  types  of  bits  were 
available  for  testing,  cross  bits  of  1 h  and  2^  in.  diameter,  and  button 
bits  of  1-3/4  and  2  in.  diameter.  These  were  tested  without  heat  and 
while  the  larger  diameter  button  bit  appeared  to  be  the  most  effective, 
this  type  of  bit  could  not  be  used  effectively  to  traverse  the  rock.  Hence, 
the  majority  of  the  preliminary  and  later  tests  were  conducted  with  the 
smaller  cross  bit,  which  was  most  effective  of  the  other  three  at  the  thrust 
and  traverse  rate  which  proved  operable  for  the  hard  rocks  tested. 
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Thrust ,  The  maximum  thrust  was  employed  at  which  the  drill  bits 


would  traverse  the  rock  without  being  stopped  in  one  position  to  drill 
a  hole.  This  also  required  operating  the  jackhammer  at  No.  3  throttle 
setting. 

Conclusions .  In  general  it  can  be  stated  that  the  values  of  breakage 
energy  calculated  by  means  of  K  which  vary  from  about  9  to  32  j/cc,  are 
probably  too  low,  because  some  of  the  very  fine  material  produced  in  the 
drilling  was  lost. 

The  chipping  energies  both  with  and  without  heat  are  not  optimum 
values.  With  improved  techniques  which  would  permit  operation  of  the 
jackhammer  at  open  throttle  and  optimum  thrust  the  energy  values  could 
be  improved  considerably.  Likewise,  if  a  cooling  spray  of  water  were 
employed  in  conjunction  with  the  flame  jet,  thermal  weakening  and  frag¬ 
mentation  could  well  be  improved. 
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Controlled  Fragmentation  Tests 

Based  upon  observations  made  in  the  preliminary  fragmentation  tests, 
i.e.,  the  large  amount  of  -200  mesh  material  produced  by  the  jackhammer, 
the  table  speed  was  increased  to  4.5  in/sec  for  most  of  the  final  series 
of  tests.  On  Missouri  Red  Granite  the  effect  of  increase  in  table  speed 
(Table  10)  was  to  decrease  the  jackhammer  specific  energy  of  rock  removal 
(w/o  heat) (Test  No.  5)  to  159  j/cc  for  one  pass  and  93  j/cc  for  three 
consecutive  passes  (Test  No.  3). 

Repeated  passes  (5)  of  the  flame  jet  (Test  No.  1)  and  3  for  the  jack¬ 
hammer  reduced  the  jackhanmer  energy  to  35  j/cc,  while  the  values  for  the 
shallow  kerfing  tests  (Tests  No.  2  and  4)  showed  some  advantage  in  the 
deeper  kerf. 

For  Charcoal  Grey  Granite  the  heat  from  both  the  flame  jet  and  a  No. 

5  weld  torch  appeared  to  have  a  detrimental  effect  rather  than  assisting 
in  the  fragmentation.  One  reason  could  be  that  the  heat  caused  micro¬ 
fractures  to  form  in  the  surface  of  the  rock,  thus  decreasing  the  efficiency 
of  transmission  of  energy  from  the  drill  bit  into  the  rock.  Shallow 
kerfing  also  showed  no  improvement  over  single  passes  with  or  without 
heat. 

The  pattern  of  behavior  for  Kitledge  Pink  Granite  was  similar  to  that 
of  the  Charcoal  Grey  Granite,  but  the  amount  of  energy  required  for  frag¬ 
mentation  was  lower.  The  Kitledge  Pink  Granite  appeared  to  be  softer  and 
easier  to  drill  than  the  other  rocks  tested.  For  the  Buena  Gabbro  the  heat 
improved  the  jackhammer  energy  consumption  for  both  six  passes  in  the  same 
groove  and  for  the  shallow  kerfing  test. 
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Table  10.  Thermomechairical  Fragmentation 
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The  pattern  was  somewhat  similar  for  Sioux  Quartzite  although  the  specific 
energy  values  were  about  one-half  of  those  for  the  Buena  Gabbro  for  the 
shallow  kerfing  test. 

The  results  for  Dresser  Basalt  showed  no  improvement  with  heat  from  the 
flame  torch  either  in  multiple  passes  or  in  shallow  kerfing.  A  rather 
marked  improvement  was  made,  however,  with  heat  from  the  No.  5  weld  torch. 

Even  with  a  high  table  speed  of  4.5  in/sec,  the  flame  jet  may  have  caused 
melting  or  softening  of  the  whole  surface  contact  layer  of  minerals  while 
the  heat  from  the  weld  torch  was  concentrated  along  a  narrow  line,  per¬ 
mitting  the  jackhamner  to  chip  more  efficiently. 

The  specific  energy  of  breakage  calculated  from  particle  size  dis¬ 
tribution  analyses  utilizing  Rittinger  numbers  determined  by  drop  crushing 
show  some  interesting  correlations.  Within  the  series  of  tests  for 
Missouri  Red  Granite,  with  the  exception  of  Test  No.  2,  there  is  rough 
correlation  between  the  calculated  and  jackhairmer  energies.  The  low  values 
for  Test  No.  1  indicate  that  the  multiple  passes  of  the  flame  jet  resulted 
in  the  formation  of  more  large  size  particles. 

Because  of  vhe  variation  in  test  conditions  there  is  likewise  only 
a  rough  correlation  between  jackhammer  energy  values  and  calculated  values 
for  the  other  types  of  rocks.  Difference  in  grain  size  (large  for  Missouri 
Red  Granite,  small  for  Buena  Gabbrc)  also  has  an  influence  on  the  energy 
required  to  fragment  the  rock. 

Test  results  indicate,  however,  that  for  given  test  conditions  the 
correlations  would  be  much  better,  particularly  if  adjusted  for  the  grain 
size  and  mineral  content  of  the  rock,  i.e.,  the  hardness  and  strength  of  the 
constituent  minerals.  Thus,  the  drop  crushing  test,  or  a  modification  of  it, 
appears  to  have  promise  as  a  test  to  predict  rock  resistance  to  fragmentation. 
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Deep  Kerfing 

The  minimum  total  specific  energy  required  to  remove  rock  (quart¬ 
zite)  by  shallow  kerfing  was  391  j/cc.  A  series  of  deep  kerfing  tests  was 
conducted  on  Charcoal  Grey  Granite,  Sioux  Quartzite,  and  Dresser  Basalt 
(Tables  11  &  12).  In  the  first  test  (No.  2)  on  quartzite  the  two  slots 
were  only  2  in.  deep  with  a  6  in.  ridge.  The  slots  were  cut  using  a 
jackhammer  assisted  with  No.  5  weld  torches.  The  ridge  was  then  removed 
by  focusing  the  flame  tip  along  one  of  the  slots  causing  a  stress  con¬ 
centration  at  the  bottom  of  the  slot,  which  served  as  the  initiating  point 
for  a  fracture  which  removed  the  ridge  at  its  base. 

For  all  of  the  remaining  tests  the  slots  were  cut  with  the  jackhammer 
alone.  In  Test  No.  3  four  holes  were  drilled  to  serve  as  loci  for  breakage. 

On  this  test,  however,  the  quartzite  spalled  to  a  wide  kerf  which  rapidly 
reached  a  stable  curved  surface  which  required  excessive  heat  for  additional 
spallation.  The  other  configurations  are  as  shown  in  the  footnote  of  Table  11. 
In  test  No.  4,  the  whole  face  of  a  30  inch  cube  was  split  off  by  running  the 
flame  jet  along  the  bottoms  of  the  outside  two  slots.  The  major  part  of 
the  energy  required  was  for  cutting  the  slots  with  the  jackhammer.  Also, 
more  heat  was  used  for  the  basalt  than  necessary  because  melting  and  solid¬ 
ification  held  the  ridge  in  place  after  it  had  fractured.  The  figures  in 
Table  11  were  adjusted  (Table  12)  to  give  the  energy  required  for  simulated 
operational  conditions,  i.e.,  one  slot  cut  per  ridge  removed  with  the  heat 
applied  to  each  slot. 

For  all  of  the  rocks  tested  the  coefficient  of  thermal  expansion  was 
adequate  to  cause  fracture.  None  of  the  tests  were  designed  to  give 
optimum  performance  of  either  the  jackhammer  or  the  flame  jet.  It  is 
estimated  that  the  total  energy  required  could  be  reduced  by  a  factor  of  3 
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Table  11 


Deep  Kerf  Tests 


Rock 

Wt  -  JH 
grams 

Vol  -  JH 
cc 

Wt  -  FJ 
grams 

Vol  -  FJ 
cc 

Specific  Energy  j/cc 
JH  FJ  JH  +  1 

1. 

Grey  Granite 

19,600 

7,180 

43,549 

15,952 

120 

5.6 

41 

2. 

Sioux  Quartzite 

5,695 

2,157 

29,318 

11,105 

387** 

10.0 

71** 

3. 

Sioux  Quartzite 

2,142 

811 

20,828 

7,889 

219 

12.3 

32 

4. 

Sioux  Quartzite 

25,220 

9,553 

162,200 

61 ,439 

210 

2.1 

31 

5. 

Basalt 

28,976 

10,614 

62,902 

21,179 

282 

6.5 

98 

Test  Geometry 

1.  two  slots  5  in.  deep  by  30  in.  long,  6  in.  ridge 

2.  two  slots  2  in.  deep  by  30  in.  long,  6  in.  ridge 

3.  four  drill  holes  \\  in.  diam,  3%  in.  deep  on  6  in.  centers 

4.  three  slots  5  in.  deep,  33  in.  long,  5  in.  ridge 

5.  three  slots  4  in.  deep,  36  in.  long,  4  in.  ridge 

**includes  acetylene  torch  energy  used  in  cutting  grooves 

Table  12 

Deep  Kerf  -  Energy  Per  Kerf 


Rock 

Wt  -  JH 
grams 

Vol  -  JM 
,  grams 

Wt  -  FJ 
grams 

Vol  -  FJ 
cc 

Specific  1 
JH  FJ 

Energy  j/cc 
JH  +  FJ 

*1. 

Grey  Granite 

9,000 

3,590 

43,549 

15,952 

120 

5.2 

31 

2. 

Sioux  Quartzite 

2,848 

1 ,078 

29,318 

11,105 

387 

10.0 

43 

4. 

Sioux  Quartzite 

8,407 

3,184 

40,550 

15,359 

210 

2.1 

38 

5. 

Basalt 

9,658 

3,538 

31  ,451 

10,589 

282 

6.5 

75 

*See  Table  VI-  for  test  geometry 
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for  the  granite  and  quartzite,  and  by  a  factor  of  4  or  5  for  the  basalt 
by  utilization  of  narrower  and  deeper  kerfs,  wider  ridges,  and  by  ap¬ 
plication  of  a  mechanical  pressure  to  wedge  the  ridges  as  thermal  cracks 
begin  to  form. 
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CHAPTER  VII 


CONCLUSIONS  AND  RECOMMENDATIONS 

The  investigations  conducted  on  the  thermal -mechanical  disintegra¬ 
tion  of  hard  rock  dealt  with  the  theory  and  measurement  of  temperature 
distribution  and  heat  coupling  of  flame  sources,  induced  thermal  stresses, 
effects  of  impact  indentors,  fragmentation  energy  determination,  and  labora¬ 
tory  tests  on  fragmentation  of  large  cubes  of  rock  using  a  flame  jet  and 
a  pneumatic  hammer.  Several  important  facts  were  disclosed  by  the  re¬ 
search,  leading  to  a  concept  for  a  thermal  fragmentation  rock  breakage 
round  for  tunnel  excavation.  The  concept  is  based  upon  a  technique  for 
thermal  fracturing  the  rock  similar  to  blasting  with  explosives,  but 
using  a  minimum  amount  of  heat  to  induce  thermal  stresses  for  fracture. 

Thermal  processes  for  excavating  rock  which  require  melting  of  a  signifi¬ 
cant  portion,  or  even  heating  a  large  percentage  of  the  rock  removed,  do 
not  appear  to  be  economically  feasible. 

Heat  Coupling  and  Temperature  Distribution.  The  complexity  of  known 
methods  of  evaluating  mathematical  solutions  for  temperature  distribution 
equations  in  three  dimensions  for  rock  materials  with  moving  heat  sources 
make  computer  calculations  impractical.  An  existing  analytical  solution 
for  a  moving  (point)  source  for  a  one-dimensional  problem  was  applied  to 
basalt  utilizing  known  thermal  properties.  The  results  showed  extremely 
high  temperature  gradients  in  the  rock  directly  beneath  the  heat  source. 

Steady  state  solutions  for  a  hemispherical  source  show  that  the  highest 
temperature  gradients  are  at  the  edge  of  the  source. 

Temperature  measurements  (for  a  flame  jet  source)  at  shallow  depths 
in  the  rocks  showed  that  the  critical  temperature  gradients  occurred  within  1  mm 
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of  the  surface.  Hence,  fractures  would  be  expected  to  originate  within 
this  depth,  and  spallation  phenomena  would  be  similar  in  many  respects  to 
that  observed  for  thermal  shock  of  brittle  ceramic  materials. 

Because  of  the  low  conductivity  of  minerals  and  the  consequent 
necessity  of  placing  thermocouples  very  close  to  the  surface,  surface 
temperatures  and  near  surface  temperature  gradients  of  the  rock  were 
difficult  to  measure.  Approximations  were  made,  however,  indicating  that 
the  heat  coupling  of  the  flame  jet  could  be  within  6  to  18  percent,  the 
lower  figure  being  the  most  reasonable. 

Thermal  Stresses.  The  basic  concepts  of  thermal  stress  analysis  in 
elastic  materials  are  well  defined.  Tractable  mathematical  equations  have 
been  derived  for  steady  state  conditions  for  several  geometries,  heat  sources, 
and  solid  homogeneous  bodies  with  temperature  independent  properties.  The 
process  of  evaluating  solutions  becomes  increasingly  complex  for  two-  and  three- 
dimensional  problems,  for  moving  heat  sources,  for  temperature  dependent 
material  properties,  and  for  most  types  of  heat  sources  which  vary  with 
time  or  distance. 

While  it  was  not  possible  to  evaluate  existing  solutions  of  three- 
dimensional  stress  equations  for  a  moving  point  source,  or  for  temperature 
dependent  properties,  existing  solutions  for  one-  and  two-dimensional 
problems  were  found  to  be  useful  in  furnishing  concepts  of  temperature 
distributions  for  the  more  complex  problems.  Thus,  a  moving  point  source 
on  the  edge  of  a  plate  indicates  that  the  highest  temperature  gradients 
exist  just  in  front  of  and  under  the  leading  edge  of  the  heat  source, 
giving  the  location  of  the  highest  stresses.  Other  configurations  have 
similar  characteri sties . 

In  thermal  shock  processes  the  high  thermal  gradients  usually  exist 
near  the  front  of  the  temperature  change.  These  give  rise  to  high  tensile, 
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compressive,  and  shear  stresses  which  will  cause  brittle  materials  to  fail. 
Some  investigators  have  described  these  types  of  thermal  stress  failure 
mathematically  by  means  of  the  Griffith  theory. 

The  nature  of  thermal  stress  fractures  observed  in  Sioux  quartzite 
indicate  that  most  probably  they  originate  at  the  (plane)  surface  and  form 
spalls  only  when  they  intersect  other  fracture  surfaces  which  allow  a 
particle  to  be  separated  from  the  rock  body.  Otherwise,  the  cracks  die 
out  with  depth  or  lateral  extension. 

In  deep  kerfing  tests  it  was  found  that  heating  at  the  bottom  of 
a  kerf  (slot)  on  one  side  of  a  ridge  caused  high  local  stress  concentra¬ 
tions  at  the  corners  of  the  groove.  When  the  stresses  and  stored  thermal 
strain  energy  reached  sufficiently  high  values  a  crack  propagated  to  the 
bottom  of  the  adjacent  kerf,  severing  the  ridge.  The  thermal  energy  con¬ 
sumption  for  removing  such  ridges  was  2  to  5  joules/cc,  and  may  be  much 
smaller  for  more  carefully  controlled  conditions  and  for  deeper  kerfs. 

Thus,  large  fracture  propagation  is  not  entirely  due  to  very  local  high 
temperature  gradients,  but  is  due  to  other  factors  as  w^ll ,  such  as  the 
overall  stress  field  created,  the  total  strain  energy,  and  heat  dissipation 
in  the  rock.  Therefore,  a  combination  of  thermal  stresses  and  geometry 
for  efficient  fracture  of  large  blocks  of  the  rock  is  necessary.  That  is, 
small  scale  fragmentation  from  a  plane  or  nearly  plane  surface  represents 
an  inefficient  utilization  of  strain  energy  and  requires  heating  of  a  large 
percentage  of  the  rock.  Residual  heat  represents  wasted  energy  and  large 
amounts  of  hot  rock  may  cause  environmental  problems  in  a  thermal  excavation 
system. 

The  condition  that  thermal  stresses  may  be  localized  without  inducing 
excessive  unwanted  stresses,  strain  or  strain  energy  in  large  portions 
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of  an  elastic  body  indicates  a  marked  advantage  may  be  achieved  in  con¬ 
trolling  the  fracture  pattern  and  the  resulting  fragmentation.  By  contrast, 
the  whole  of  a  compressive  specimen  in  a  testing  machine  is  strained  before 
it  ruptures.  Likewise,  in  most  blasting  operations  all  of  the  blasted  rock 
is  subjected  to  large  stresses  and  strains,  resulting  in  a  high  degree 
of  fragmentation. 

The  facts  that  the  location  of  thermal  stresses  can  be  controlled 
with  respect  to  a  favorable  geometry  of  the  rock,  and  that  geometries 
may  be  created  which  permit  fracture  control  with  low  energy  expenditures 
are  the  bases  for  the  design  of  a  thermal  fragmentation  tunnel  round 
described  later. 

Mechanical  Fragmentation.  The  effects  of  shapes  of  indentors  and 
friction  between  the  indentors  and  rock  were  investigated  mathematically. 

It  was  found  that  the  force  required  for  penetration  was  least  for  sharp 
wedge-shaped  indentors  with  no  friction.  Resistance  to  penetration  in¬ 
creases  rapidly  for  increased  friction  and  bluntness  of  the  nose  of  the 
wedge,  as  well  as  with  the  included  angle  of  the  indentor.  Thus,  the 
sharper  the  angle  the  easier  penetration  is  achieved.  However,  sharp 
edges  wear  so  rapidly  that  a  compromise  must  be  made  between  ease  of 
penetration  and  the  angle  which  will  give  optimum  wear.  For  semicyl indrical 
indentors,  the  force  required  for  penetration  increases  rapidly  with  depth, 
but  quickly  levels  off  as  the  depth  approaches  a  value  of  half  of  the 
radius  of  curvature  of  the  indentor. 

The  slip-line  hypothesis  based  upon  the  theory  of  plasticity  gives 
an  indication  of  why  heat  treated  rocks  may  be  more  difficult  to  chip 
than  solid  rocks.  That  is,  if  near  surface  fractures  are  created  by 
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heat,  then  the  friction  along  the  cracks  may  be  less  than  that  within  the 
solid,  thus  restricting  the  area  of  plastic  deformation,  failure  and 
possible  chipping.  Also,  the  existence  of  fractures  in  the  surface 
impedes  the  transmission  of  stress  waves  which  normally  are  responsible 
for  some  of  the  chip  formation.  Thus,  for  some  rocks  which  have  micro¬ 
fractures  near  the  surface  due  to  thermal  stress,  the  minerals  at  the  area 
of  impact  are  pulverized,  and  small  particles  are  produced  rather  than 
large  fragments. 

Further  experimental  studies  of  individual  indentation  processes  will 
be  required  to  check  the  applicability  of  the  theory  which  has  been 
developed. 

Fragmentation  Energy.  Of  the  laws  which  have  been  proposed  by  re¬ 
searchers  to  represent  the  relationship  between  particle  size  distribution 
and  specific  fragmentation  energy,  Rittinger's  hypothesis  and  its  modi¬ 
fications  appear  to  offer  the  best  basis  for  evaluation  of  thermal - 
mechanical  fragmentation  processes.  For  some  types  of  thermal  fragmentation 
which  produce  only  relatively  large  fragments,  the  flaws  in  the  unbroken 
solid,  and  the  manner  in  which  fragmentation  occurs  have  an  effect  upon 
the  apparent  energy  distribution. 

The  technique  of  employing  a  drop  ctuomng  device  to  evaluate 
fragmentation  energy  furnishes  a  means  of  approximating  the  apparent 
particle  size-energy  relationship  for  calculation  of  an  effective  Rittinger 
constant,  Kr>  This  in  turn  may  be  employed  to  calculate  the  approximate 
useful  energy  consumed  in  the  thermal  fragmentation  process  where  the 
particle  size  distribution  is  similar.  Because  producing  very  fine 
material  consumes  so  much  energy,  even  small  amounts  of  finer  particles 
have  a  dominant  effect  in  fixing  the  value  of  ,  and  particle  shape  factors 
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employed  mostly  for  larger  particles  have  a  minor  effect  on  calculations. 

For  the  drop  crushing  and  thermal -mechanical  fragmentation  tests  con¬ 
ducted  in  this  project  is  was  found  that  50  percent  or  more  of  the  crushing 
energy  was  utilized  in  producing  -200  mesh  material.  Thus,  it  will  be 
profitable  to  develop  rock  fragmentation  methods  which  will  produce  a 
minimum  of  fine  material. 

Also,  the  value  of  Kr  gives  an  approximate  measure  of  the  energy  re¬ 
quired  to  break  the  rock.  It  is  planned  to  standardize  and  refine  the 
crushing  test  so  that  it  can  be  used  as  a  reliable  index  for  measuring 
rock  resistance  to  breakage. 

Thermal -Mechanical  Fragmentation.  The  controlled  test  results  show 
that  for  some  rocks  thermal  treatment  with  the  flame  jet  made  them  more 
susceptible  to  surface  mechanical  fragmentation  by  a  pneumatic  hammer, 
while  other  rocks  became  more  resistant.  Both  single  pass  and  shallow 
kerf  methods  of  fragmentation  required  relatively  large  amounts  of  heat 
from  the  flame  .jet,  and  smaller,  but  significant  amounts  from  a  #5  weld 
torch.  Both  types  of  tests  required  that  the  rock  removed  be  heated  to  a 
critical  temperature,  and  only  part  of  the  heat  was  utilized  in  fracture. 

The  remaining  residual  heat  is  believed  to  represent  a  large  percentage 
of  the  total  heat  transmitted  to  the  rock. 

Deep  Kerfing.  For  deep  kerfing  tests  much  less  heat  per  unit  volume 
of  rock  removed  was  required,  because  only  a  small  portion  of  the  rock  was 
heated.  As  previously  described,  for  5  in.  deep  kerfs  and  ridges  6  in.  wide 
the  total  heat  energy  transferred  into  the  rock  was  from  2  to  5  joules/cc. 
Thus,  a  combination  of  utilization  of  thermal  strain  energy  with  proper 
geometry  to  provide  nuclei  for  crack  initiation  at  the  point  of  high  thermal 
stress,  plus  guidance  of  the  fracture  to  create  ‘large  fragments,  permits  an 
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efficient  use  of  the  thermal  energy. 

Based  upon  the  results  of  deep  kerfing  it  is  planned  that  the 
feasibility  of  a  thermal  fragmentation  round  be  investigated.  As  with 
explosive  energy,  the  most  efficient  usage  of  thermal  energy  for 
fracturing  for  excavation  purposes  is  obtained  when  the  energy  is  re¬ 
leased  within  the  rock,  and  fracture  processes  are  enhanced  by  the  presence 
of  free  faces  which  allow  fracture  propagation.  That  is,  some  of  the 
basic  principles  of  rock  breakage  by  explosives  may  be  applied  to  the 
design  of  a  thermal  fragmentation  round. 

Since  the  energy  required  for  pneumatic  drilling  varies  from  75  to 
200  joules/cc  for  medium  hard  to  hard  rocks,  an  optimum  design  of  a  round 
will  require  a  minimum  of  mechanical  removal  of  rock  for  forming  a  kerf 
and  drilling  holes  for  placement  of  heating  elements.  This  may  be 
accomplished  by  means  of  a  shallow  round  utilizing  the  principle  of  a 
burn  cut  round. 

One  concept  provides  that  the  "cut"  portion  of  the  round  may  be  a 
vertical  or  horizontal  slot  near  the  bottom  a  tunnel  face  with  holes 
drilled  adjacent  thereto,  the  breakage  proceeding  outward  therefrom.  The 
depth  of  the  round  will  be  determined  by  the  size  of  blocks  which  can  be 
handled  by  the  muck  removal  system.  The  spacing  of  holes  also  will  be 
governed  by  the  size  of  blocks  desired,  the  amount  of  heat  required  versus 
spacing  and  depth,  and  the  properties  of  the  rock.  The  required  hole 
diameter  is  a  factor  related  to  the  critical  size  of  heating  elements, 
the  heat  flux  and  the  heat  transfer  properties  of  the  rock.  The  frag¬ 
mentation  sequence  contemplated  is  similar  to  an  explosive  round.  If  the 
round  is  relatively  shallow,  most  of  the  broken  rock  will  fall  onto  a 
conveyor  or  within  reach  of  a  simple  mucking  machine. 
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Laboratory  feasibility  studies  are  planned  on  three  to  four  foot 
cubes  of  three  or  four  types  of  hard  rock,  including  basalt,  granite, 
and  quartzite.  Heating  elements  of  the  electrical  resistance  type  have 
furnished  sufficient  heat  to  fracture  small  blocks  in  preliminary  tests. 

The  kerf  can  be  cut  readily  and  holes  drilled  with  a  pneumatic  hammer. 

Problems  expected  are  in  the  design  of  structurally  strong  heating 
elements  which  will  stand  up  under  both  laboratory  and  field  conditions. 

For  initial  field  tests  in  a  subsequent  phase  of  the  research  program 
jackhammers  and  electric  elements  may  be  mounted  on  separate  arms  of 
available  automated  jumbos. 

Recommendations.  In  view  of  the  findings  of  the  investigations 
described  in  this  report  the  following  are  deemed  to  be  areas  worthy  of 
additional  research: 

1.  Continue  theoretical  investigations  of  thermal  shock,  heat 
transfer,  temperature  distribution,  stress  distribution,  and  partitioning 
of  thermal  energy  between  strain  energy,  residual  energy,  and  other  means 
of  absorption. 

2.  Continue  studies  of  evaluation  of  fragmentation  energy  of  rocks 
by  thermal  and  mechanical  means,  including  further  refinements  of  the 
"drop  crushing"  technique  to  provide  a  measure  of  rock  resistance  to 
fragmentation, 

3.  Continue  studies  of  rock  fragmentation  by  thermal,  mechanical 
and  by  a  combination  of  the  two  processes,  utilizing  in  line  acetylene 
heaters,  infrared  heat,  quenching,  and  roller  type  cutters  and  cones. 

4.  Make  a  study  of  the  basic  thermal  and  mechanical  properties 

of  minerals  and  rocks  which  affect  their  resistance  to  different  techniques 
for  fragmentation. 
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5.  Assess  other  means  of  exploiting  the  surface  microfracturing 
of  rock  caused  by  heating. 

6.  Investigate  the  feasibility  of  a  thermal  fragmentation  round  for 
a  tunnel  in  hard  rock  as  outlined  under  Thermal -Mechanical  Fragmentation. 
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APPENDIX  A 


ROCK  PROPERTIES 
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Table  A-2.  Rock  Composition 


Rock 

Missouri  Red  Granite 

Charcoal  Grey  Granite 

Kitledge  Pink  Granite 

Black  Granite 
(Buena  Gabbro) 

Sioux  Quartzite 

Dresser  Basalt 


Composition  (percent  by  volumel 


Quartz 

Potash  feldspar 
Plagioclase  feldspar 
Hornblende 
Hypersthene 
Magnetite 


32.5 

42.5 
25.0 
trace 
trace 
trace 


Feldspar,  plagioclase 
(sodic  end)  and  microcl 
Quartz 
Hornblende 
Biotite 

Clays,  Zircon,  Apatite, 
Sphene 


ine  63.0 
16.0 
17.0 
3.0 


1.0 


Labrador ite 

Augite 

Olivine 

Quartz 

Chert 

Mica 

Feldspar  (Labradorite) , 
microcl ine 
Augite 
Magnetite 


45.0 

40.0 

15.0 

99.0 

0.5 

0.5 


50.0 

40.0 

8.0 
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APPENDIX  B 

TEMPERATURE  DISTRIBUTION  CURVES 
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Fig.  B-l.  Subsurface  Temperature  Measurements  -  Dresser  Basalt 
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Fig.  B-2.  Subsurface  Temperature  Measurements  -  Dresser  Basalt 
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Fig.  8-4.  Subsurface  Temperature  Measurements  -  Dresser  Basalt 
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Fig.  B-10.  Subsurface  Temperature  Measurements  -  (Mo.  Red  Granite) 
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Subsurface  Temperature  Measurements  -  (Mo.  Red  Granite) 
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Fig.  B-12.  Subsurface  Temperature  Measurements  -  (Mo.  Red  Granite) 
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Fig.  B-14.  Subsurface  Temperature  Measurements  -  (Mo.  Red  Granite) 


Fig.  B-15.  Subsurface  Temperature  Measurements  -  (Mo.  Red  Granite) 
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Fig. B-18.  Subsurface  Temperature  Measurements  -  (Mo.  Red  Granite) 
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Fig. 8-19.  Subsurface  Temperature  Measurements  -  Buena  Gabbro 
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Fig.  B-21.  Subsurface  Temperature  Measurements  -  Buena  Gabbro 


stance  of  Torch  from  Thermocouple  (cm) 


Distance  of  Torch  from  Thermocouple  (in.) 
Fig.  B-22.  Subsurface  Temperature  Measurements  -  Buena  Gabbro 
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Fig.  B-24.  Subsurface  Temperature  Measurements  -  Buena  Gabbro 
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Fig.  B-25.  Subsurface  Temperature  Measurements  -  Buena  Gabbro 
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ig.  B-26.  Subsurface  Temperature  Measurements  -  Buena  Gabbro 
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NOWACKI'S  solution  for  half-space 


APPENDIX  C 


Nowacki's  Solution  for  Half-Space 


Suppose  that  on  the  plane  =  0  which  bounds  the  elastic  semi-space, 
a  concentrated  heat  source  of  unit  intensity  is  moving  in  the  origin  of  the 
coordinate  system  in  the  direction  of  the  axis.  It  is  assumed  that 

the  plane  £3  =  0  is  thermally  insulated  and  free  of  traction. 

The  heat  conduction  equation  for  this  case  is 


+  20  fer T(  V  = 
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where 


k  =  —  ,  the  diffusivity,  and 
pc 

Q*Ur)  =  6Ui)6(€2)«(Cj) 

The  solution  of  equation  (C-l)  for  this  problem  is  known 
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T(Cr)  = 
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R  =  (5?  +  +  $)H 
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Solution 
two  solutions, 
heat  source  in 


to  the  stress  problem  is  obtained  in  the  form  of  a  sum  of 
i  and  s.  The  solution  s  is  that  for  the  action  of  a  moving 
the  infinite  space  and  is  given  by 
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Here 

is  the  stress  tensor 

r  =  (£|  +  Kl)h 

k  =  m~ 

4ttk 

u  =  6,  the  shear  modulus 
_  1+v 

m  -  T=v  a 

v  =  Poisson's  ratio 

a  =  coefficient  of  thermal  expansion 

Note  that  on  the  plane  £3  =  0,  |i  =  0  ,  ai3  =  a23  =  0.  The  only 
non-zero  component  of  the  stress  tensor  on  this  plane  is  o33(£  1,  £2,  0) 

=  -p(£i,  £2).  The  solution  5  'is  that  for  the  elastic  semi-space  £3  >  0 
in  which  T  =  0,  the  only  loading  being  the  normal  traction  o33  =-p(£x,  £2) 
on  the  plane  £3  =  0.  Hence,  the  stresses  5-,  are  obtained  from  the 
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integral  expression 
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Iri  this  expression,  J..(£i,  £2,  £3;  m,  02,  0)  denote  the  stresses 

*  v 

at  the  point  (£  )  due  to  the  action  of  a  concentrated  unit  force  situated 
at  the  point  (m ,  02,  0)  on  the  plane  bounding  the  semi-space,  and  are 
given  by 
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The  stresses  in  the  semi-space  due  to  the  action  of  a  source  moving 
on  the  bounding  plane  f,3  =  0  are  obtained  by  adding  the  stresses  given 
oy  equation  (c-3)  and  (C-4)  as 
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APPENDIX  D 


STRESS  AND  FORCE  SOLUTIONS  FOR  FOUR  INDENTORS 


Indentor  No.  1 

A.  Smooth  Surface  (Fig.  III-4(A) ) 

1.  Maximum  shear  theory:  0<  8  <  j  ,  n  =  £ 
c*n  =  2K(  1 

T  -  0 

n 

F  =  2ABonsin6  (force/unit  length) 


2.  Internal  friction  theory 

°n  ■  2K^-^  exp(2»tan») 


T  =  0 

n 


K  cot$ 


F  =  2ABonSine 

B.  Rough  Surface  (Fig.  I II-4(B) ) 

1.  Maximum  shear  theory:  0<  p  <  ^  ,  n  =  0 

°n  = 


F  =  2AB(ansin3+insin0) 


where  0 


B,  if  e  <  4 


7T  ^  TT 

T>  lf  6  4 

AB  =  Edge  of  Indentor,  if  6  <  £ 

Edge  of  False  Nose,  if  B  ?  ^ 

2.  Internal  friction  theory:  0  <  n=0 

on  =  kcos*  [  ( l+sin}>)exp(2^tan4>) -1  ] 
tn  =  K(l+sin<j>)exp(2ytan$) 

F  =  2AB( ansin0+TncosQ) 
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where  t>  =  6,  if  g  £  ^ 
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,  if  6  > 
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Edge  of  Indentor,  if  B  <  ^  - 


2 


Edge  of  False  Nose,  if  8  >  -  |r 

C.  Frictional  Surface  (Fig.  Ill -4(A) ) 


1.  Maximum  shear  theory:  0  <  8  <  ^,  0  <n  <  jj-  . 

=  K(  l+2i|>+sin2n) 

t  =  Kcos2n 
n 

F  =  2AB(onsin8+Tncosg) 

2.  Internal  friction  theory:  0<6<^--|-,0<n<^-+|- 

°n  =  "  s''n(<f>-2n)]exp( 2^tan4> ) -Kcot4> 

Tn  =  ~cbs|nl)  C0S(^-2ri)exP(2’i,tan4) 

F  =  2AB  (nnSin8+xnCOS8) 


1 1 .  Indentor  No.  2 

The  stress  and  force  equations  for  indentor  No.  1  can  be  used  for  this 
indentor  with  the  following  modifications. 

A.  Smooth  Surface  (Fig.  1 1 1 -4{ C) ) 

1.  Along  BC  -  the  same  as  I-A,  except  4<=i^,  and  AB=BC. 

2.  Along  AB  -  the  same  as  I-B,  except  ip=ii^+ip^  . 

B.  Rough  Surface 

1.  6  <  £  -  j  (Fig.  III-4(C))  (Note:  for  maximum  shear  theory  <t»=0.) 

(a)  Along  BC  -  the  same  as  I-B,  except  and  AB=BC. 

(b)  Along  AB  -  the  same  as  I-B,  except  dJ=ip1+^2 . 

2.  b  >  J  -  £  (Fig.  1 1 1-4(D) ) 

The  same  as  I-B 
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III. 


C.  Frictional  Surface  {Fig.  III-4(C):  e  <  j  ) 

1.  Along  BC  -  the  same  as  I-C,  except  <»**  and  AB=£C. 

2.  Along  A8  -  the  same  as  I-B,  except  ^=ih+^2. 
Indentor  No.  3  (r  =  1) 

A.  Smooth  Surface  (Fig.  1 1 1 -5(A) ,  AB=0) 

1.  Maximum  shear  theory:  0  <  ct  <  ir/2,  n  s  J  . 
on  =  2K(l+i|f+0)  where  0  <  o  <  a 

*n  =  0  _ 


F  =  4  j oncos(a-e)d0 
ctior 

n  a^'*linft  ~~  exP[2(4'+©)tan<i>j  -Kcot<j> ,  where  0  <  e  <  a 


2.  Internal  friction  theory:  0<a<ii-»n=-+- 

L  4  2 


T  =  0 

n 


F  =4  / oncos(a-o)d0 
o 

B.  Rough  Surface  (Fig.  III-5(B) 

1.  Maximum  shear  theory:  0  <  a  <  n  =  0. 
(a)  Along  BC  (a  *  J) 

on  =  K[  l+2( ip+G)J  ,  where  0  <  0  <  a  -  ^ 


x_  =  K 


a-i 


F  =  2  /  [ancos(a-o)+Tnsin(a-0)]d0 


(b)  Along  AB 


on  =  Ktl+2('i'+0)]  ,  where  ©  = 


a-jj-  ,  if  a  > 
0  ,  if  a  < 


Tn  =  K 


F  =  /2  AB(0n+Tn) 


2.  Internal  friction  theory:  0  <  a  <  ^,  n  =  0. 

0-3 


(a)  Along  BC  (  a  >  £  +  |) 

on  =  K  ot4>{ 0  +s in<J) )axp [2 (^j-i-©) tan ] -1 } 
tn  =  K(l+sin<j>)exp[2(t|rt-0)tan<|>] 
where  0  <  0  <  a  £  -  | 

JT  "I 

F  =  2  J  [oncos(a-0)+rn  in(a-0)]d0 


(b)  Along  AB 


on  =  Kcot4>  •[(  l+sin4>)expt2(^+Q)  tan<J>] -1 } 


rn  =  K(1+sin<t>)exp[2(ijj+0)tan<}>] 


a 

where  G  =  J 
0 


1  _  £  ifrt>2L  +  i 
^  2*  •  1  cx  -  4^  ^  2 


»  if  a<^-  + 


F  =  2AB  [onsia(|  -  f)  *  Vos(jj-  -  £)] 

C.  Frictional  Surface  (Fig.  I II-5(A) ) 

1.  Maximum  shear  theory:  0<a<^-,0<n<^. 
(a)  Along  BC 


on  =  K[1  +2(^+0)  Hsi  n2n] ,  where  0<G<  a  +  n  - 

t  =  Kcos2n 
n  tt 

otfrr  4 

F  =  2  f  [ancos(a-0)+ipSin  (a-0)]dG 
o 

(b)  Along  AB 

on  =  Kcoti)>{(l+sin^)exp[2(ijj+0)tan4>]-l} 

Tp  =  K(l+sin$)  exp  [2(^+0)  tan<f>] 
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where  0  =  otfr|"  T 
0 


if  a>  I  -n 
1  f  a  <  J-  -n 


F  =  ft  AB  (an+rn) 

2.  Internal  friction  theory:  0<a<|,0<nSjj-+|- 


(a)  Along  BC 

an  *  ^  cbs]'n^  [~W  '  Sl  n(4>-2n)]expt2(^+G)  tanip] -Kcotip 

Tn  cos (4)-2n) exp [2(^+0) tan$J‘ 

where  0  S  Q  S  »n  -  |  . 


a+n  - 


rr 


F  =  2 


I 


[ancos(a-O)+Tnsin(a-0)]dO 


(b)  Along  AB 

on  =  Kcot{  (1  +s  i  n«f>)exp[2(^+0)  tan^] -1  > 

in  =  K0+sin<t)exp[2(ij/+O)tan(}>] 

where  0  =  “  fr<  '  J  ’  f  ’  if  a  -  l  +  I  "  n 
0  ,  if  a  <  }  +  f  -  n 

F  =  2AB[onsin(J  -  |)  +  xn©s(J  -  |)1 

IV.  Indentor  No.  4 

The  stress  and  the  force  equations  for  indentors  no,  1  and  no.  3  can  be  used 
for  this  indentor  with  the  following  modifications: 

A.  Smooth  Surface  (Fig.  III-5(C),  AB=0) 

1.  Along  CD  -  the  same  as  I-A,  except  AB=CD. 

2,  Along  ABC  -  the  same  as  1 1 I-A,  except  B  =  ^  -  a 
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B.  Rough  Surface 

1.  M  l  -  |  (Fig.  111-5(0)) 

(a)  Along  CD  -  the  same  as  I-B,  except  AB=CD. 

(b)  Along  ABC  -  the  same  as  III-B,  except  g  =  -  a. 

2.  3  >  f  -  |  (Fig.  III-5(C),  n  =  0) 

The  same  as  I-B,  excrpt  AB=AD  (AD  is  the  edge  of  the  false  nose, 

not  shown  in  Fig.  III-5(C) ) . 

C.  Frictional  Surfar-  (Fig.  II I-5(C) ,  g  <  J  -  $■  ) 

1.  Along  CD  -  the  same  as  I-C,  except  AB=CD. 

2.  Along  ABC  -  the  same  as  III-C,  except  6  =  j  -  a 


D-6 


APPENDIX  E 

PARTICLE  SIZE  ANALYSIS 
DROP  CRUSHING  TESTS 


Particle  Size  Analysis  -  Drop  Crushing  -  Missouri  Red  Granit 
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Fig.  E-2.  Particle  Size  Analysis  -  Drop  Crushing  -  Missouri  Red  Granite 


Test  No. :  3 

Rock:  Missouri  Red  Granite 


Particle  Diameter  -  in. 

£-3.  Particle  Size  Analysis  -  Drop  Crushing  -  Missouri  Red  Granite 


E-4.  Particle  Size  Analysis  -  Drop  Crushing  -  Missouri  Red  Granite 


Particle  Diameter  -  in. 

Fig.  E-5.  Particle  Size  Analysis  -  Drop  Crushing  -  Missouri  Red  Granite 


Particle  Diameter  -  in. 

.  E-6.  Particle  S1*e  Analysis  -  Drop  Crushing  -  Missouri  Red  Granite 
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Particle  Size  Analysis  -  Drop  Crushing  -  Missouri  Red  Granite 


Fiq.  E -9.  Particle  Size  Analysis  -  Drop  Crushing  -  Missouri  Red  Granite 
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Particle  Diameter  -  in. 

Particle  Size  Analysis  -  Drop  Crushing  -  Charcoal  Grey  Granite 
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Particle  Diameter  -  in. 

Size  Analysis  -  Drop  Crushing  -  Charcoal  Grey  Granite 
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Charcoal  Grey  Granite 


Fig.  E-15.  Particle  Size  Analysis  -  Drop  Crushing  -  Charcoal  Grey  Granite 


Particle  Diameter  -  in. 

-16.  Particle  Size  Analysis  -  Drop  Crushing  -  Charcoal  Grey  Granite 
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Particle  Diameter  -  in. 

Fig.  E-18 .  Particle  Size  Analysis  -  Drop  Crushing  -  Kitledge  Pink  Granite 
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Particle  Size  Analysis  -  Drop  Crushing  -  Kitledge  Pink  Granite 


Size  Analysis  -  Drop  Crushing  -  Kitledge  Pink  Granite 
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Fig.  £-22.  Particle  Size  Analysis  -  Drop  Crushing  -  Kitledge  Pink  Granite 
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Particle  Size  Analysis  -  Drop  Crushing  -  Kitledge  Pink  Granite 


Test 


luacuad 


Particle  Diameter  -  in. 

25.  Particle  Size  Analysis  -  Drop  Crushing  -  Black  Granite  (Buena  Gabbro) 


Particle  Diameter  -  in. 

Fig.  t-26.  Particle  Size  Analysis  -  Drop  Crushing  -  Black  Granite  (Buena  Gabbro) 
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Particle  Size  Analysis  -  Drop  Crushing  -  Black  Granite  (Buena  Gabbro) 
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Particle  Size  Analysis  -  Drop  Crushing  -  Black  Granite  (Buena  Gabbro) 
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Particle  Size  Analysis  -  Drop  Crushing  -  Black  Granite  (Buena  Gabbro) 


Fig.  E-  31 .  Particle  Size  Analysis  -  Drop  Crushing  -  Sioux  Quartzite 
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Fig.  E-33.  Particle  Size  Analysis  -  Drop  Crushing  -  Sioux  Quartzite 
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Fig.  E -34 .  Particle  Size  Analysis  -  Drop  Crushing  -  Sioux  Quartzite 


Particle  Size  Analysis  -  Drop  Crushing  -  Sioux  Quartzite 


Fig.  E-37.  Particle  Size  Analysis  -  Drop  Crushing  -  Sioux  Quartzite 
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Test  No. :  39 

Rock:  Dresser  Basalt 


Fig.  E-41 .  Particle  Size  Analysis  -  Drop  Crushing  -  Dresser  Basalt 
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Fig.  E-43.  Particle  Size  Analysis  -  Drop  Crushing  -  Dresser  Basalt 
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PARTICLE  SIZE  ANALYSIS 
FRAGMENTATION  TESTS 


Particle  Diameter  -  in. 
alysis  -  Fragmentation  -  Missouri  Red  Granite 


Particle  Diameter  -  in. 

Particle  Size  Analysis  -  Fragmentation  -  Missouri  Red  Granite 


Particle  Diameter  -  in. 

Fig.  F-3  Particle  Size  Analysis  -  Fragmentation  -  Missouri  Red  Granite 


Fig.  F-4.  Particle  Size  Analysis  -  Fragmentation  -  Missouri  Red  Granite 


Particle  Diameter  -  in. 

Fig.  F-5 .  Particle  Size  Analysis  -  Fragmentation  -  Missouri  Red  Granite 


Particle  Diameter  -  in. 

Particle  Size  Analysis  -  Fragmentation  -  Missouri  Red  Granite 


Particle  Diameter  -  in. 

F-7.  Particle  Size  Analysis  -  Fragmentation  -  Charcoal  Grey  Granite 


Fragmentation  -  Charcoal 


Fig.  F-9.  Particle  Size  Analysis  -  Fragmentation  -  Charcoal  Grey  Granite 
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Particle  Diameter  -  in. 

Fig.  F-ll .  Particle  Size  Analysis  -  Fragmentation  -  Kitledge  Pink  Granite 


e  Diameter  -  in. 

Fragmentation  -  Kicledge  Pink  Granite 


Test  No. 


Particle  Diameter  -  in. 

Fig.F-15.  Particle  Size  Analysis  -  Fragmentation  -  Black  Granite  (Buena  Gabbro) 


Fig.  F-l 9  -  Particle  Size  Analysis  -  Fragmentation  -  Sioux  Quartzite 


Test  No . :  2C 
Rock:  Dresser  Basalt 
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-20.  Particle  Size  Analysis  -  Fragmentation  -  Dresser  Basalt 


Fig.  F-21.  Particle  Size  Analysis  -  Fragmentation  -  Dresser  Basalt 
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APPENDIX  G 


Energy  Transfer  in  the  Percussive  Drilling  of  Rock 


In  percussive  drilling,  a  piston  is  caused  to  accelerate  along  a 
cylinder  by  the  application  of  compressed  air  to  one  side  of  the  piston. 

The  piston  then  impacts  a  drill  steel,  transferring  all  or  part  of  its 
energy  to  t  le  steel  in  the  form  of  a  strain  wave.  It  is  known  that  the  amount 
of  energy  transferred  to  the  steel  is  proportional  to 


e2 (t)dt 


(G-l) 


in  the  first  incident  wave,  where  e(t)  is  the  strain  at  a  particular  point 
in  the  steel  and  ( t^-t ^ )  is  the  duration  of  the  strain  wave.  The  above  as¬ 
sumes  that  the  reflected  wave  does  not  interfere  with  the  incident  wave  at 
the  measuring  point. 

Rewriting  (H-l)  as  an  equation,  we  have 


(G-2) 


where  K  is  the  appropriate  constant  of  proportionality  and  W.  is  the  amount 
of  energy  transferred.  The  constant  K  can  be  determined  by  impacting  the 
drill  steel  with  a  projectile  of  known  energy  and  dividing  by  the  integral 
of  the  square  of  the  resulting  strain  wave,  under  the  assumption  that  100 
percent  of  the  energy  is  transferred  from  the  projectile  to  the  steel. 

The  amount  of  energy  transferred  to  the  rock  will  be  determined  by 
monitoring  the  longitudinal  strain  waves  in  the  drill  steel.  This  was  ac¬ 
complished  by  Lindquist  and  Anderson  (Ref.  76). 

The  output  of  the  bridge  circuit  is  a  voltage  proportional  to  the 
longitudinal  strain  at  the  bridge  location. 


G-l 


> ,  BBM  ffiyiwii-i 


That  is, 

v ( t )  =  Kl£(t)  16-3) 
where  v(t)  is  the  output  voltage  and  Ki  is  the  constant  relating  strain 
to  voltage.  Equation  (H-2)  can  now  be  written  as 


/tf  t^ 

c2(t)dt  =  K  Kf  f  v2(t)dt 

c.  A. 

1  i 


(,  G-4 ) 


Assuming  that  the  reflected  wave  does  not  interfere  with  the  incident  wave, 
the  energy  in  the  first  incident  wave  is 


W 


/  V. 


(t)dt 


(G-5) 


i  i 


where  T  is  the  time  interval  during  which  the  particular  strain  pulse  is 
being  measured  and  where  the  subscript  v.  denotes  the  fact  that  we  have  an 
incident  wave  and  the  second  subscript  tells  us  we  have  the  first  wave  of 
a  series.  Similarly,  for  the  energy  in  the  first  reflected  wave  we  have 


W 


r  i 


ri 


(t)dt 


(G-6) 


where  the  subscript  r  denotes  a  reflected  wave.  Note  that  the  limits  on 
the  integrals  in  (H-5)  and  (H-6)  indicate  that  integration  is  performed 
only  over  the  duration  of  each  strain  wave.  Further,  since  we  have  assumed 
that  the  incident  and  reflected  waves  do  not  interfere,  the  -  integration 
intervals  are  disjointed  in  time.  This  can  be  accomplished  by  proper  selec¬ 
tion  of  strain  gage  bridge  location  and  length  of  drill  steel  for  a  parti¬ 
cular  hammer.  Expressions  for  the  energy  in  the  second  incident  and  re¬ 
flected  waves,  respectively,  are 


lG-7) 


i 
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G-2 


(6-8) 


The  amount  of  energy  transferred  to  the  rock  on  the  first  incident  wave  is 
merely  the  difference  between  the  energy  in  the  first  incident  wave  and  the 
energy  in  the  first  reflected  wave.*  That  is, 


where  Ej  is  energy  transferred  on  the  first  incident  wave,  Similarly,  the 
energy  transferred  by  the  second  wave,  E2;  the  third  wave  E 3 ,  etc.,  is 


(G-1 0) 


where  n  is  the  number  of  terms  required  to  specify  the  amount  of  energy 
transferred  for  each  piston  blow.  Experimental  evidence  shows  this  number 
to  be  less  than  ten.  Since  at  most  ten  reflections  are  significant,  the 
energy  of  each  piston  blow  is  transferred  before  the  next  piston  blow  occu:  . 

From  the  above  discussion,  it  is  seen  that  energy  transferred  for  each 
piston  blow  i s  (Ey ) , 

EI'  £  E1  =  -  v  -T  “ij  -£  (G-") 

J  =  1  J  =  1  FT  jTH 

For  a  conventional  drill  steel  and  bit,  the  incident  wave  is  compressive 
and  the  reflected  wave  is  tensile  whether  the  bit  is  in  intimate  contact 
with  the  rock  or  whether  it  is  not  in  contact,  the  latter  occurring  at  low 
thrusts.  The  reason  for  this  sign  reversal  is  that  the  effective  impedance 

*It  is  assumed  that  there  are  no  losses  in  the  steel. 


G-3 


of  the  rock  is  less  than  that  for  the  drill  steel  and  that  additional  re¬ 
flection  occurs  because  of  the  poor  coupling  between  the  wedge  bit  and 
the  rock.  Consequently,  the  output  voltage  of  the  bridge  circuit  will 
be  positive  when  the  wave  is  compressive  and  negative  when  the  wave  is 
tensile  (Fig.  G-1 ) 


Fig.  G-1.  Typical  Output  of  Bridge  Circuit  for  Each  Piston  Blow 

From  Figure  G-1  we  see  that  the  bridge  output  is  a  series  of  voltage 
waveforms  which  are  disjointed  in  time.  Next,  define  x(t)  to  be  the 
positive  voltages  and  y(t)  the  negative  voltages  as  shown  in  Figure  G-2. 
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Fig,  G-2.  Bridge  Output,  v(t),  Resolved  into  Positive  Voltage 
Output,  x{t)  and  Negative  Voltage  Output,  y(t) 


A  device  to  separate  v(t)  into  x(t)  is  easy  to  construct  and  will  not  be 
discussed.  For  the  purpose  of  reference  we  will  choose  to  cal!  this  de¬ 
vice  a  Quan  tizer. 

Since  x(t)  and  y(t)  are  disjointed  in  time,  they  are  necessarily 
orthogonal.  From  (H-ll),  we  have 


n  n  n  n 


(t)dt  =  KK i 


X2(t)dt 


(G-12) 


The  total  amount  of  energy  transferred  to  the  rock  for  a  large  number  of 
piston  blows  is  difficult  to  obtain  directly,  but  the  average  power  (P  ) 
transferred  into  the  rock  can  be  obtained  easily  by  dividing  the  expression 
in  (H-12)  by  the  time  interval,  T:  that  is 


Pave  =  KK 


l  j /T*2( t)dt  -  KK?  1/  y2(t)dt  (G-13) 


Since 


i  f  x2(t)dt 
J  n 


equals  the  square  of  the  rms  value  of  x(t) (similarly  for  y(t)),  we  have 
Pave  =  KK?  K  ~  KK*  ^y  =  -  VJ) 

A  A 

where  V^,  and  Vy  denote  the  rms  values  of  x ( t )  and  y(t)  respectively. 
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